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INTRODUCTION 


This  report  contains  a  collection  of  technical  papers  describing 
research  in  adaptive  control  supported  by  AF05R  contracts  R9620-83-C-01 07 
and  FH9620-8U-C-00M . 

The  basic  objective  of  this  research  program  is  to  establish  the 
theoretical  foundations  and  performance  limitations  for  adaptive  control 
applications  to  large  space  structures  (LSS).  An  important  element  of  the 
research  is  to  examine  implementation  concepts  which  can  lead  to  appropriate 
hardware  development. 

The  program  was  originally  formulated  in  late  1982  in  response  to  the 
increasing  concern  that  performance  robustness  of  Air  Force  LSS  type  systems 
would  be  inadequate  to  meet  mission  objectives.  In  particular, 
uncertainties  in  both  disturbance  spectra  and  system  dynamics 
characteristics  (both  time  varying  and  stochastic  uncertainty)  usually 
significantly  limit  the  performance  obtainable  with  fixed  gain,  fixed 
architecture  controls.  The  use  of  adaptive  type  controls,  where 
disturbances  and/or  plant  models  are  identified  prior  to  or  during  control, 
gives  systems  designers  more  options  for  minimizing  the  risk  in  achieving 
performance  benchmarks. 

The  research  was  originally  directed  toward  real-time  adaptation  of 
the  estimator  form  using  variable  order  lattice  filters  to  construct  the 
desired  compensation.  Early  in  the  research,  however,  lack  of  a  well- 
developed  robustness  theory  for  adaptive  mechanizations  required  a 
reexamination  of  the  problem  at  a  more  fundamental  level,  i.e.,  development 
of  model  and  disturbance  uncertainty  bounds  for  which  adaptive  algorithms 
would  exhibit  (stable)  desired  performance.  Toward  this  end  there  have  been 
two  major  accomplishments: 

(1)  Development  of  Theory:  In  examining  the  possible  use  of  "rapid” 
adaptive  control  it  was  necessary  to  generate  new  theory  of  use  on  large 
space  structures.  This  theory  accounts  for  the  effect  of  unmodeled  dynamics 
with  distributed  parameter  systems,  such  as  flexible  space  structures,  and 
extends  current  adaptive  theory  in  several  directions. 

In  the  first  place,  current  adaptive  theory  provides  conditions  for 
"global"  stability,  l.e.,  bounded- input ,  bounded-output  stability  with  no 


limitation  on  the  size  (or  spectrum)  of  the  bounded- 1 nputs  (e.g., 
disturbances  and  references).  Secondly,  the  theory  is  limited  to  finite- 
dimensional  linear  systems.  This  latter  condition  cannot  be  satisfied  by  a 
flexible  space  structure,  which  is  a  distributed  parameter  system.  Also, 
the  disturbance  and  reference  Inputs  effecting  the  spacecraft  have  limited 
magnitudes  and  spectrums  and  these  limits  are  known,  although  not  precisely. 
The  theory  we  have  developed  circumvents  those  difficulties  by  providing 
conditions  for  "local"  stability,  i.e.,  limitations  in  input  size  and 
spectrum  are  accounted.  The  theory  also  allows  for  a  distributed  system  as 
well  as  providing  quantifiable  bounds  on  permissible  model  error.  These 
results  extend  the  state-of-the-art  in  adaptive  theory  beyond  the  current 
limits. 

(2)  Methodology  Development:  The  use  of  "slow"  adaptive  control, 
which  is  more  practical  than  rapid  adaptive  control  in  most  space 
applications,  necessitated  a  new  methodology  development  merging  key  ideas 
in  parameter  estimation,  system  identification,  and  robust  control  design. 

By  "slow"  we  mean  that  there  is  sufficient  time  to  run  batch  identification 
before  the  control  system  is  modified.  The  methodology  we  have  developed 
resolves  a  long  standing  problem  with  adaptive  systems  of  this  type,  namely, 
the  means  to  provide  a  guaranteed  level  of  performance  given  an  "identified" 
model  of  the  system  together  with  the  model  error  between  the  system  and  the 
identified  model.  In  fact,  our  methodology  generates  performance  vs.  model 
error  tables  (to  be  stored  in  the  computer)  from  which  the  control  design  is 
immediately  obtained.  Moreover,  the  order  of  the  control  design  is 
determined  strictly  on  the  basis  of  model  error  and  performance  demand, 
rather  than  trial  and  error  as  has  been  suggested  in  the  past. 
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ISSUES  IN  THE  ADAPTIVE  CONTROL  OF  LARGE  SPACE  STRUCTURES 


by 

Robert  L.  Kosut  and  Michael  G.  Lyons 
Integrated  Systems,  Inc. 

Palo  Alto,  CA  94301 


Abstract 


This  paper  examines  some  basic  theoretical  and  practical  issues  in 
the  adaptive  control  of  large  space  structures  (LSS).  Particular  attention 
is  paid  to  the  practical  Issues  of  model  error,  decentralization,  and 
subsystem  performance  allocation.  It  Is  concluded  that  the  currently 
available  theory  of  adaptive  control,  which  is  based  on  global  stability, 
centralized  information,  and  perfect  modeling,  is  not  well  suited  for  an 
LSS.  A  direction  for  future  research  is  suggested  which  is  based  on  a 
theory  of  local  stability  for  the  adaptive  system. 


•Research  supported  by  AFOSR  under  contract  F49620-84-C-0054 . 
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1 .  INTRODUCTION 

The  high  performance  requirement  a  of  Large  Space  Structnrea  (LSS) 
together  with  potentially  large  nncertaintiea  in  the  aystem  model,  motivate 
the  nee  of  an  adaptive  control  ayatem.  Although  a  great  variety  of  adaptive 
control  achemes  exiat  for  lumped  parameter,  amall  acale  ayatema,  e.g.  [1], 
theae  methodologiea  cannot  be  directly  applied  to  the  LSS  becauae  of  the 
following  iaauea: 

(1)  Model  error  -  The  actual  ayatem  ia  a  diatributed  parameter  ayatem, 
theoretically  of  infinite  dimenaion,  whereaa  the  adaptive  acheme 
muat  be  baaed  on  a  reduced  order  model  (ROM)  of  finite  dimenaion. 
Thia  diacrepancy  introducea  one  Lind  of  model  error,  the  effect  of 
which  ia  often  referred  to  aa  'apillover.'  Another  claaa  of  model 
errora  are  thoae  attributable  to  uncertainty  in  parametera  (e.g., 
mode  ahapes),  neglected  non-linearitiea ,  and  other  uncertain 
unmodeled  phenomenon  (e.g.,  reaidual  modes  included). 

(2)  Decentralized  control  -  In  aome  caaea  the  physical  size  and 
complexity  of  the  LSS  makes  it  impractical  to  use  a  centralized 
control  structure  due  to  considerations  of  actuator/sensor  costs, 
system  reliability,  and  computational  requirements,  as  well  as  the 
step-wise  deployment  (and  removal)  of  sub-sections. 

(3)  Performance  allocation  -  Since  the  performance  requirements  are 
stringent,  it  is  necessary  to  allocate  performance  in  an  efficient 
manner  so  that  sub-systems  can  help  One  another. 

In  this  paper  we  explore  the  above  issues  from  the  point  of  view  than  an 
LSS  can  be  represented  as  a  large-scale  interconnected  system  [2] .  The 
interconnection  model  used  is  composed  of  a  number  of  uncertain  subsystems 
which  are  linked  to  other  subsystems  by  an  interconnection  operator,  which  is 
also  uncertain.  Uncertainty  in  the  subsystems  and  interconnections  is 
expressed  by  using  the  notion  of  a  conic  model  [3]-[4],  i.e.,  representing  a 
complicated  uncertain  dynamic  system  as  belonging  to  a  set  of  systems 
generated  from  simpler  dynamic  systems. 

By  using  this  representation  the  issues  enumerated  above  can  be  brought 
within  a  single  framework  which  facilitates  the  analysis  and  synthesis  of 
adaptive  controllers  as  discussed  in  [5}-[6]. 


IS  2 .  BACKGROUND  DISCUSSION 


The  development  of  a  design  methodology  for  adaptive  control  of  LSS 
involves  many  different  issues.  A  comprehensive  discussion  of  the  theoretical 
and  practical  problems  involved  in  both  LSS  control  and  adaptive  control  is 
well  beyond  the  scope  of  this  paper.  In  this  section  we  present  a  very 
selective  discussion  of  the  issues  that  seem  most  relevant. 
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2.1  Adaptive  Control 

Adaptive  methods  have  achieved  a  great  amount  of  success  in  producing 
stable,  convergent  adaptive  controllers  and  adaptive  observers  for  systems 
whose  structure  is  known  and  whose  parameters  are  constant  but  poorly  known  or 
slowly  time-varying.  Adaptive  schemes  may  be  direct,  as  shown  in  Fig.  1, 
i.e.,  the  available  control  parameters  are  directly  adjusted  (adapted)  to 
improve  the  overall  system  performance  or  indirect,  as  shown  in  Fig.  2,  i.e., 
the  system  parameters  are  identified  (based  on  the  assumed  system  structure) 
and  the  control  commands  are  generated  from  these  parameter  estimates  as 
though  they  were  the  actual  values. 

The  use  of  such  methods  on  distributed  parameter  or  large  scale  systems, 
like  LSS,  is  greatly  limited  by  the  modeling  problem — the  adaptive  scheme  must 
be  based  on  a  ROM  (Reduced  Order  Model)  of  the  actual  system  and,  hence,  the 
order  of  the  model  is,  and  must  remain,  substantially  lower  than  the 
controlled  system  [7]~(8]. 

The  crux  of  the  problem  with  adaptive  control  is  to  guarantee  that  the 
adaptive  controller  that  is  designed  on  the  ROM  will  not,  through  a 
combination  of  spillover  and  model  uncertainty,  diverge  and  ultimately  go 
unstable.  This  problem  is  extremely  difficult  and  has  only  recently  been 
addressed,  even  in  the  general  context  of  adaptive  system,  e.g.,  [5], [6], [9], 


2.2  LSS  Modeling  and  Model  Error 

Traditionally,  control  design  is  based  on  models  of  the  system  which  have 
been  validated  thoroughly  by  extensive  testing.  Since  the  structural 
integrity  of  LSS  does  not  permit  ground  tests,  the  usual  approach  to  modeling 
is  not  feasible  for  LSS.  The  primary  difficulty  is  validating  the  model, 
i.e.,  determining  a  Quantifiable  measure  of  model  uncertainty.  In  this 
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regard,  the  LSS  has  the  interesting  property  of  being  an  infinite  dimensional 
system,  at  least  theoretically  so,  bnt  practically  of  very  high  order.  Thus, 
the  order  of  the  design  model  and  controller  is  not  Lnovn  a  priori.  In 
consequence,  adaptive  controllers  for  LSS  should  be  not  only  parameter 
adaptive,  but  order  adaptive  as  veil.  This  leads  naturally  to  the 
consideration  of  order-recursive  lattice  adaptive  control  [10]. 

A  very  natural  means  to  determine  model  error  (sometimes  referred  to  as 
plant  uncertainty)  is  to  perform  an  experiment  which  compares  the  model  with 
data  from  the  actual  system  (or  plant).  If  there  is  no  error  between  the 
model  and  the  plant,  then  we  have  perfect  knowledge  of  the  plant.  Normally, 
the  situation  is  the  opposite — the  error  is  non-xero  and  represents  how  close 
the  model  is  to  the  plant.  If  we  quantify  this  experiment,  by  defining  a 
specific  measure  of  the  error  size,  then  this  gives  a  sensible  statement  as  to 
model  accuracy.  For  example,  during  experimental  modeling  using  system 
identification  methods  the  model  uncertainty  is  measured  as  the  difference 
between  the  measured  output  and  the  model  output.  Bounding  this  model  error, 
for  all  possible  input/output  pairs,  results  in  a  set  characterization  of 
plant  uncertainty.  For  example,  a  set  description  of  an  uncertain  LTI  plant 
is  to  define  a  ball  in  the  frequency  domain.  The  center  of  the  ball  is  the 

nominal  plant  model,  and  the  radins  defines  the  model  error.  This  set  model 

description  is  one  type  of  a  more  general  set  description,  referred  to  as  a 

conic-sector  [4] .  The  uncertainty  in  the  plant  induces  an  uncertainty  in  the 

input/output  map  of  the  closed-loop  system  which  can  again  be  characterized  by 
a  conic  sector.  Performance  requirements  for  the  control  system  can  be 
translated  into  statements  on  the  conic  sector  which  bounds  the  closed-loop 
systems,  making  it  possible  to  check  whether  a  given  design  meets 
specifications,  and  providing  guidelines  for  robust  controller  design,  e.g.. 

Ill). 


2.3  Decentralized  Control  for  LSS 

In  the  context  of  LSS  control  design,  what  we  mean  by  a  decentralized 
control  is  the  following:  The  control  system  is  made  up  of  a  number  of  sub¬ 
controllers  (local  controllers)  which  have  limited  authority  over  the  LSS  and 
which  use  limited  information  about  the  LSS.  The  limitations  on  control 
authority  and  the  information  pattern  are  the  main  features  of  the 
decentralized  control  problem.  The  general  structure  of  such  a  decentralized 
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control  system  is  illustrated  in  Figure  3.  The  dashed  lines  indicate  a 
partial  information  eicbange.  e.g.,  the  local  controller  receives  reference 
coamands  (or  discretes)  from  a  higher  level  control  (the  coordinator)  and/or 
information  from  other  local  controllers  in  the  form  of  an  'aggregated'  state. 

In  a  decentralized  control  we  also  need  to  determine  the  effect  of 
partial  information  on  closed-loop  performance.  The  two  hinds  of  system 
ignorance,  i.e.,  model  uncertainty  and  partial  information,  can  be  viewed 
nnder  one  framework  by  considering  the  controlled  LSS  as  an  interconnected 
system,  e.g. ,  [2] . 

An  interconnected  system  is  a  system  which  consists  of  several  subsystems 
interacting  through  various  interconnection  operators.  The  key  feature  used 
here  is  that  knowledge  about  the  subsystems  and  interconnection  operators  is 
incomplete. 

Techniques  for  decentralized  synthesis  more  or  less  adhere  to  the 
following  steps  (see,  e.g.,  [12]). 

Step  1:  (Decomposition)  Identify  the  individual  subsystems  and  the 
interconnection  constraints  between  subsystems. 

Step  2:  (Local  analysis)  Design  the  local  controller  so  that  each 

individual  subsystem  satisfies  specified  local  requirements. 

Step  3:  (Global  analysis)  Verify  that  the  interconnection  of  the 

individual  subsystems  satisfies  specified  global  requirements. 

Step  4:  (Robustness)  Verify  that  the  total  system  performance  is  robust 
with  respect  to  failures,  disconnections,  parameter  changes, 
etc. 

In  practice,  these  'steps'  overlap  and  iterations  are  required.  However,  the 
first  step,  decomposition,  is  necessary  to  begin  the  design  process. 

There  are  several  methods  available  for  decomposition.  These  may  be 
broadly  grouped  into  generic  categories  based  on:  time-scale  separation, 
frequency  separation,  and  performance  properties  (e.g.,  observability, 
controllability,  quadratic  cost,  etc.).  All  df  these  can  be  viewed 
graphically  as  well  as  tabularly,  and  many  physical  systems,  LSS  included, 
possess  all  three  types  of  decompositions  (see,  e.g.,  [13],  Chapters  1  -  III 
for  complete  exposition). 

It  is  important  to  emphasize  that  many  decomposition  methods  are  purely 
mathematical  and  may  decompose  the  system  into  simpler  numerical  problems 
convenient  for  parallel  or  distributed  processing.  In  the  LSS  environment. 
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decomposit ions  normally  result  because  a  natural  separation  is  physically  or 
geographically  present  between  functional  components  of  the  system.  For 
example,  consider  the  following  task  oriented  decomposition: 

(1)  High  authority  actuators  and  sensors  with  low  bandwidth  for  rigid 
body  control 

(2)  High  authority  actuators  with  bandwidths  to  20  Hz  for  slewing,  with 
possibly  a  series  connected  low  authority  actuator  for  final  small 
motion  slew  correction 

(3)  Medium  authority  actuators  and  sensors  for  vibration  isolation  of 
disturbances 

(4)  Low  authority  actuators  for  isolation  of  critical  structural 
subsystems  (mirrors,  focal  plane,  etc.) 

(5)  Very  low  authority  actuators/sensors  for  active  damping  or  resonant 
absorption. 

Some  decompositions  result  from  spatial  differences;  weak  dynamic 
interaction  effects  can  be  easily  identified.  A  decomposition  also  occurs 
from  temporal  differences;  phenomena  occurring  at  different  time-scales,  e.g.. 
a  separation  between  fast  and  slow  mode*  or  between  low  frequency  and  high 
frequency  effects.  For  example,  groups  of  the  modes  can  be  separately 
controlled  by  separate  controllers  which  do  not  destabilize  each  other. 
Specific  combinations  of  weak  dynamic  coupling  and  separation  of  slow  and  fast 
modes  can  often  be  identified,  e.g..  Figure  4. 

In  many  cases,  delegation  of  control  authority  is  'politically' 
practical.  It  is  unrealistic  to  assume  that  the  manufacturer  of  one  device 
will  ever  design  another  manufacturer's  device,  or  even  that  both  will 
delegate  complete  authority  to  a  aystems  house.  The  only  communication 
possible  in  this  case  is  to  assess  each  manufacturer  with  specifications  so 
that  the  operating  devices  do  not  compete.  For  example,  simultaneous  on-orbit 
assembly  of  different  parts  of  the  LSS  may  be  accomplished  using  temporary 
vibration  control  systems  built  by  various  manufacturers.  This  motivates  a 
design  specification  which  includes  tolerances  that  allow  for  some  variety, 
auch  that  the  overall  differences  do  not  impair  on-going  missions  or 
constructions  in  other  parts  of  the  LSS. 
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The  inpot /output  view  of  interconnected  system  is  completely  compatible 
with  conic-sector  model  nncertsinty  descriptions.  Representative  theorems 
(see.  e.f.,  [2])  refer  to  interconnected  systems  of  the  form 


*»  '  "i  '  £  "*j  Ti 


i  *  1.  ...» 


where  G.  .  ....  G_  are  the  subsystem  operators,  and  H,,,  ....  H  are  the 
*  ■  11  mm 

interconnections.  The  key  features  of  the  theorems  are: 


(1)  If  the  subsystems  and  interconnections  have  quantifiable 
input/output  properties,  e.g.,  conic-sector  bounds  and/or  passivity, 
then  the  total  system  will  exhibit  performance  properties  directly 
attributable  to  the  subsystems  and  interconnections,  e.g.,  conic- 
sector  bounds  and/or  passivity. 

(2)  The  total  system  properties  can  be  obtained  by  combinations  of 
subsystem  and  interconnection  properties. 

This  latter  point  is  extremely  important  in  the  decentralised  setting.  This 
means  that  one  subsystem  can  'help'  overcome  the  deficiencies  of  another 
system.  Furthermore,  this  also  gives  a  clue  to  the  question  of  allocating 
subsystem  performance  in  an  efficient  way  so  that  a  desired  overall 
performance  is  achieved. 


3  .  BOM'S! .  .(  NON  AJUn  I VI  j  f  ONTBOI 

Before  attempting  to  develop  a  methodology  for  adaptive  control  of 
uncertain  (decentral ized)  systems,  it  is  logical  to  consider  the  non-adaptive 
case  first.  In  fact,  the  adaptive  design  procedure  should  build  on  the  robust 
design  procedure. 

Historically,  research  in  robust  control  theory  has  proceeded  from  an 
input/ontput  view  of  systems,  e.g.,  [3], [14].  The  more  recent  of  these 
results  [4], [15]  are  variations  on  the  Small  Gain  Theorem  [3],  The  theorem 
asserts  that  if  the  'loop  gain'  of  a  feedback  system  is  less  than  unity,  then 
the  closed-loop  system  is  stable.  Bowever,  to  properly  utilize  the  theorem  it 
is  necessary  to  isolate  the  source  of  the  model  error.  This  is  accomplished 
by  what  is  called  a  'loop-transformation.' 

Many  variations  on  loop-transformations  are  available  [11],  [IS],  and 
small  gain  theory  can  be  applied  to  analyze  the  robustness  of  criteria  other 
than  stability,  e.g.,  tracking  response  transients.  Moreover,  the  technique 
can  also  be  used  to  assess  the  impact  of  various  kinds  of  error  sources.  For 
example,  typical  sources  of  model  error  in  spacecraft  systems  include: 

(1)  numerical  errors  due  to  approximate  modeling  techniques,  e.g.,  high 
order  NASTRAN  models. 

(2)  actual  parameter  changes  in  the  LSS  ,e.g.,  thermal  effects,  gravity, 
spacecraft  and  antenna  dimensions,  mass  distributions,  etc. 

(3)  unmodeled  dynamics,  e.g.,  effect  of  reduced  order  modeling, 
neglected  residual  modes  ('spillover'),  neglected  actuator/ sensor 
dynamics,  non-linearities,  etc. 

(4)  incomplete  data  obtained  from  on-earth  testing,  e.g.,  partially 
assembled  structures  in  simulated  zero-g. 

Therefore,  uncertainty  in  the  baseline  model  arises  from  both  actual  causes 
and  intentional  approx isiat ions  of  complicated  phenomena.  In  many  cases,  these 
are  indistinguishable. 

Negative  results  from  the  robustness  analysis  may  warrant  redesign  of  the 
controller,  or  even  upgrading  the  reduced  order  design  model  if  adequate 
robustness  cannot  be  achieved  (e.g.,  [11]). 

3 .1  Application  to  LSS 

Consider  a  controlled  spacecraft,  as  depicted  in  Figure  5,  with  the 
following  model : 
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y,  “  >'  4  ®  « 


n  ■=  sensor  noise 
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Figure  5.  Block  Diagram  of  Controlled  Spacecraft 


Dynamic  Model 

y  -  M(s)n  +  d  ,  d  =  disturbance 
Controller 

n  =  C(s)(r-ys>  ,  r  =  reference 

where  M(s)  is  a  finite  dimensional  transfer  function  matrix  model 
representing  the  dynamics  of  the  actuators  and  spacecraft  and  C(s)  is  the 
transfer  function  matrix  of  the  controller. 

Let  the  actual  dynamics  be  represented  by 
y  *  P(s)u  +  d 

where  P(s)  is  not  necessarily  finite  dimensional.  For  example,  P(s)  can 

be  either  a  high  order  NAS TRAN  model  or  represent  data  from  the  actual  system, 

whereas  M(s)  is  the  reduced  order  control  design  model.  Thus,  the 

controlled  output  is: 

y  =  (I  +  PC)-1  d  +  (I  4  PC)"1  PC  (r-n) 

H  =  I  -  B  . 
yd  yr  yd 

(The  complex  variable  's'  has  been  suppressed  for  brevity  of  notation, 
unless  needed  for  clarification.) 

Suppose  that  M  is  a  reduced  order  model  of  P.  Let 
P  *=  M  +  Ar 

where  represents  the  effect  of  neglected  residual  modes.  Thus,  following 

[11],  the  closed-loop  response  is: 

-  -1  -1  -1 
+  MC)  Ml  ■*  iC(I  +  MC)  1 )  A  C(1  +  MC)  1 
yd  yd  .  t  t 

B  «  (1  +  MC) 


there  U  .  it  the  nominal  transfer  function  with  no  model  error,  i.e., 
yd 

A  *  0. 
r 

Similarly,  we  can  examine  the  way  in  which  other  kinds  of  model  error 
enter  into  the  closed-loop  dynamics.  The  spacecraft  model,  for  example,  may 
not  include  actuator  dynamics.  This  omission  can  be  represented  by  the  model 
error  form, 

P  *  M(I  +  A^) 

where  A#  represents  the  deviation  from  the  dynamics  of  an  actuator  with 

infinite  bandwidth.  In  this  case,  the  closed-loop  response  is: 

H  .  =  S  .  -  (I+MC)"1  M  [I+A  (I  +  CM)"1CM)_1A  Cd+MC)'1 
yd  yd  a  a 

Applying  the  Small  Gain  Theorem  [11].  [15],  the  actual  spacecraft  system  is 
stable  if: 

(1)  The  spacecraft  model  M  is  stabilized  by  C,  i.e.,  the  transfer 


functions  (I  +  MC) 


C(I  +  MC) 


(I  +  MC)  *M,  and 


(I  +  CM)  CM  are  all  stable. 


and  either: 


(2)  Reduced  order  model  errors  are  bounded  by 

I  |A^(  jin)  1 1  <1/  IIC(I+MC)_1(jw)||,  u  l  0 
or 

(3)  Actuator  model  errors  are  bounded  by 

I  |A  (jit)  1 1  <  l/ll<I+CM)_1CM(ja.)l|,  u  l  0 
where  the  norm  II. 1 1  can  be  any  matrix  norm.  Typically,  the  maximum 
singular  value  o(.)  is  used.  However,  this  may  be  unnecessarily 
conservative;  other  measures  are  available,  e.g.,  the  Perron  eigenvalue  [4]. 
The  selection  of  the  appropriate  matrix  norm  will  be  examined. 

In  [7],  there  are  several  examples  of  these  robustness  tests  using  the 
ACOSS  model,  CSDL-I.  Table  1  summarizes  these  robustness  tests  for  generic 
model  errors  bounded  in  singular  value  by 
o[A(  ja»)]  £  6(«)  ,  u  2  0 

where  6(w)  is  determined  from  input/output  tests,  e.g.,  RMS  tests.  Table  1 
shows  the  stabil itv  martins,  denoted  by  &gB»  defined  as  the  maximum  bound  oi 
model  error,  which  (at  the  specified  location,  e.g.,  actuator,  sensor,  etc.) 
ensures  stability.  Thus, 

6(w)  <  6^(10)  •  «  2  0. 

Note  that  the  tests  shown  presume  only  one  location  for  uncertainty. 
Bounds  on  simultaneous  errors  at  different  locations  are  easily  obtsined  [4], 
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STABILITY  AND  PERFORM  AN  CF.  ROBUSTNESS  MARGINS 


17].  Again,  the  upper  singular  value  norm  o(.)  can  be  replaced  by  any  other 
matrix  norm  I  I  .  1 1  . 


Performance  Robustness 

A  similar  procedure  will  be  used  to  determine  an  upper  bound  on  model 
error  to  ensure  a  specified  level  of  performance,  i.e.,  performance 
robustness.  Let  desired  performance  be  defined  by 


o[(By<J  -  B^)(jfc>)]  <  p(«o)c[Byjj( j«o)] 


Thus,  p(u>)  specifies  a  bound 

nominal  B  . .  For  example,  if 
yd 

by 


on  the  relative  deviation  of  B  .  about  the 

yd 

the  effect  of  reduced  order  modeling  is  bounded 


o[Ar<jw)]  <  p(w)tp(u)+l]"1/o[C(I+MC)"1(jw)] ,  w>  0 

then  the  desired  performance  robustness  is  guaranteed. 

Similar  expressions  can  be  obtained  as  a  result  of  other  sources  of  model 

error,  e.g.,  sensor  model  error.  Table  1  summarizes  these  performance 

robustness  tests  for  generic  model  error.  The  table  shows  the  performance 

margins .  denoted  by  &  ,  defined  as  the  maximum  bound  on  model  error  (at  the 
“  ‘  pm 

specified  location)  which  ensures  the  specified  performance  tolerance.  Thus, 

6(w)  6  (u>)  ,  «  >  0  . 

pD 

guarantees  performance  robustness.  This  also  guarantees  stability,  since, 

6  (<■>)  <  6  (<i>)  ,  ui  >  0. 

pm  sm 


4.  ADAPTIVE  CONTROL 

It  is  compelling  to  pass  directly  from  the  preceding  notions  about  robust 
control  to  the  following  indirect  'adaptive'  control  algorithm: 


Identification 

Step  1:  Using  input/output  data  estimates  the  free  model  parameters, 
thereby  obtaining  M. 

Step  2:  Using  the  same  input /output  data  in  Step  1  obtain  the  upper 
bound  on  the  residual  (unmodeled)  dynamics. 
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Design 

Step  3:  Vsing  the  model  M  from  Step  1  along  with  the  model  error  bound 
in  Step  2,  determine  a  control  compensator  C  such  that 
performance  is  bounded  above  by  a  desired  level.  If  no  such 
control  C  can  be  found,  return  to  Step  1  and  upgrade  the  model 
fidelity. 

Reconfiguration 

Step  4:  Reconfigure  the  existing  control  in  accordance  with  Step  3. 

Step  5:  Return  to  Step  1  and  repeat. 

Although  this  process  appears  entirely  reasonable,  there  are  several  open 
questions.  In  particular: 

(1)  What  is  the  best  identification  procedure  for  Step  1?  For  example, 
what  are  the  advantages  of  output  error,  ARMA  models,  lattice  forms, 
etc? 

(2)  Bow  is  the  control  design  in  Step  3  actually  implemented  in  Step  4? 
If,  for  the  example,  the  new  controller  is  put  in  place  instantly, 
then  there  may  very  well  be  a  transient  introduced  such  that  system 
performance,  although  stable,  ia  unacceptable.  If,  on  the  other 
hand,  the  new  controller  is  gradually  phased  in,  e.g., 

u  =  ( l-a)u  +  a  u 

'  '  OLD  NEW 

where  a  varies  slowly  from  ’O'  to  *1*,  then  the  question  is:  how 
slowly? 

(3)  If  in  Step  3  no  control  is  found  to  satisfy  performance,  then  how  is 
the  model  fidelity  upgraded?  Should  we  add  modes  to  the  model?  Or 
perhaps  the  test  for  model  error  is  too  conservative.  If  so,  then 
how  can  we  build  a  hierarchy  of  model  error  tests? 

(4)  In  Step  3,  what  is  the  design  procedure  for  selecting  the  feedback 
given  a  nominal  model  M,  a  bound  on  model  error,  and  a  desired 
performance  level? 

With  the  exception  of  (2),  all  the  above  issues  pertain  to  robust  (non- 
adaptivc)  design.  This  should  not  be  surprising,  since  one  must  assume  the 
existence  of  a  tuned  robust  control,  which  could  be  attained  by  the  adaptive 
system.  Thus,  in  order  to  prove  the  existence  of  a  tuned  control,  it  must  be 
possible  to  design  one. 


App 1 1  cat  Ion 


Graphs  1-5  show  magnitude  of  model  error  vs.  frequency  us;lng  various 
models  compared  to  the  input-output  data  for  the  CSDL  H2  structure.  In  all 
the  figures  we  have  plotted  the  "optimal"  stability  margin  determined  on  the 
basis  of  an  accurate  knowledge  of  the  first  18  modes  (up  to  2  Hz).  In  order 
to  achieve  guaranteed  performance  levels  near  specification,  it  is  necessary 
that  model  error  be  significantly  smaller  than  the  optimal  margin  in  the  2 
Hz  range. 

Graph  1  shows  model  errors  for  an  accurate  10  mode  model  (0.6ii  Hz)  and 
an  accurate  18  mode  model  (1.77  Hz).  To  achieve  performance  it  is  necessary 
to  identify  modes  11  to  modes  18  (0.81  to  1.77  Hz)  in  order  to  reduce  the  10 
mode  model  error.  Graph  2  shows  a  blow-up  of  the  2  Hz  region  of  interest. 

In  Graphs  3  and  ii,  respectively,  we  show  the  result  of  two  different 
13  mode  ID  procedures.  Both  procedures  use  2  mode  models  where  the 
parameter  estimates  are  obtained  from  data  which  is  filtered  over  narrow 
frequency  bands.  The  ID  procedure  in  Graph  5  sweeps  overlapping  modal  bands 

II- lll,  1 3“  16,  and  15-18.  The  models  are  then  added  together  to  form  the  13- 
mode  model.  Although  the  model  error  is  very  close  to  Ideal  (1 8-mode  model) 
below  1  Hz,  there  is  considerable  peaking  near  modes  17  and  18.  The 
procedure  used  to  obtain  Graph  i»,  however,  shows  a  more  uniformly  small 
error.  In  this  case  a  13-mode  model  is  obtained  by  sweeping  through  narrow 
non-overlapping  frequency  bands,  i.e.,  modal  bands  11-111,  12-15,  1 B*4 16,  and 

III- 18.  Using  this  latter  13-mode  model  it  is  possible  to  obtain  guaranteed 
performance  very  close  to  optimal  even  though  modes  111,  15,  and  17  are  not 
completely  known. 


In  summary,  a  fundamental  issue  in  the  design  of  an  adaptive  controller 
for  an  LSS  is  robustness  to  reduced  order  modeling,  parameter  uncert a inty ,  and 
unmodeled  dynamics.  Current  theory,  which  provides  conditions  for  stability 
(or  convergence)  of  adaptive  systems,  is  limited  to  global  stability  and 
relies  on  the  passivity  of  a  particular  subsystem  operator  [1],  [5].  In  this 
regard,  since  the  LSS  has  an  inherent  passivity  with  co-located 
actuators/sensors,  it  ia  natural  to  exploit  this  for  adaptive  control,  e.g., 
[16],  117).  However,  the  presence  of  actuator  dynamics  easily  violates  the 
passivity  requirement. 

For  example,  consider  actuator  error 

P  =  M(I  +  A^) , 

A  unknown  but  stable  where  M,  the  nominal  model,  is  in  model  form,  i.e., 

a 

,  ,  -1 
M  =  sB'<t(s  I  +  2sZfl  +  D  )  <*>'B 

with  Q  and  z  diagonal  matrices  of  modal  frequencies  and  modal  dampings, 
respectively,  <p  is  an  orthonormal  matrix  whose  columns  are  the  approximate 
mode  shapes;  and  B  is  the  actuator  influence  matrix.  Thus,  M  is  passive, 
in  fact,  M  is  positive  real  (PR),  i.e.,  M  is  exponentially  stable,  and 

{j[M(j«o)  :=  J  k  [M  ( j«u)  +  M(-jui) '  J  >0,  «d  e  R 

where  k(.)  is  the  smallest  eigenvalue.  For  scalar  systems,  [i[M(jw)]  - 
Re [M( ju)] .  In  [5]  it  is  shown  that  P  remains  passive  if  M  is  passive,  A# 

is  stable,  and  A  is  bounded  such  that 

a 

o[A>(jw)]  <  g[X(ju)]/o[M(jM)] .  w  e  R 

As  discussed  in  [5],  this  bound  is  very  conservative  and  easily  violated 
by  even  the  most  benign  actuator  dynamics,  e.g.,  a  second  order  actuator 
model.  On  the  other  hand,  the  SPR  condition  is  a  sufficient  condition  and  not 
necessary  for  stability  of  the  adaptive  system.  Further,  practical  evidence 
from  actual  applications  supports  the  fact  that  SPR  is  not  needed  to  provide 
high  performance  adaptive  systems,  e.g.,  [18]. 


The  need  for  the  SPK  condition  can  be  eliminated  by  considering  local 
stability  ratber  than  global  stability  [6].  Local  stability  refers  to 
stability  where  known  restrictions  exist  for  the  system  external  inputs, 
uncertain  parameters,  and  unmodeled  dynamics.  For  example,  persistent 
excitation  induces  exponential  stability  1 2 8] .  Since  an  exponentially  stable 
system  is  inherently  robust,  it  is  logical  to  expect  that  unmodeled  dynamics 
could  be  tolerated.  In  [6]  several  mechanisms — including  persistent 
excitation —  are  examined  which  ensure  stability  of  the  adaptive  system, 
without  SPR,  provided  certain  other  restrictions  are  enforced,  e.g.,  slowly 
varying  signals,  approximate  SPR,  and  restricted  signal  magnitudes. 

Application  to  LSS 

Consider  the  LSS  adaptive  control  system 
y  =  d  +  Pu 

where  d  is  the  disturbance;  P  is  the  mxm  transfer  function  matrix  across 
colocated  actuators/sensors.  Let  each  of  the  m  control  signals  be  given  by 

Ui  “  ”®iyi  ’  i  =  * 

Av 

where  e^(t)  is  the  adaptive  gain  at  each  colocated  station.  The  objective 
of  the  adaptive  controller  is  to  suppress  the  vibrations  due  to  the 
disturbance,  while  achieving  a  specified  damping  in  certain  critical  modes. 
Note  that  the  adaptive  control  structure  is,  in  effect,  decentralized.  A 
typical  parameter  update  law  [1]  is, 

• 

0.  *  g.e.«4  *.  >  o 

i  i  i’o  i 

*i  *  WjfsJyj 

Local  stability  conditions  exist  (see  (€]  for  details)  if  the  system  S: 

I — >  <*i  ,e2'yl  *  *kown  in  *table  where 

L  ■  diag  ^  ,  ....  j 


If  •  ••  9 
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M  K  y*Hi5*  4  e*H2y* 

Hj  «=  W(1  +  PC#)_1P 
H2  -  (I  ♦  PC.)'1 
W  *  diag  (1  ,  ....  W  ) 

X 

C,  =  diag  (0*  .  ....  0*) 

y#  «  Hjd,  e#  =  Wy# 

•  • 

The  constants  0  ,  ...»  0^  are  the  tnned  adaptive  gains;  y#  and  e#  are 

the  responses  of  the  corresponding  tnned  system.  Global  stability  is 
guaranteed  if  e0  and  y0  approach  zero  asymptotically.  and  are 

exponentially  stable,  and  is  SPR.  The  SPR  condition  on  is  not 

needed  for  local  stability.  It  turns  out  that  S  is  exponentially  stable  if 
y*  i«  persistently  exciting  and  e#  is  sufficiently  small  (61 • 


Figure  6.  Feedback  System  S 


5.  CONCLUDING  REMARKS 

In  this  paper  we  have  briefly  discussed  some  of  the  practical  issues 
involved  in  the  adaptive  control  of  an  LSS,  e.g.,  decentralization,  model 
error,  and  performance  allocation.  The  conclusion  is  that  existing  adaptive 
theory  needs  to  be  radically  revised  if  useful  engineering  tools  are  to 
emerge.  A  particular  direction  for  further  research,  as  advocated  here,  is 
the  further  development  of  a  theory  of  local  stability  for  adaptive  systems 
[6],  Ve  have  shown  that  such  a  theory  is  compatible  with  conic-model 
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represent  at i ons  14}  and  interconnected  system  theory  12} ,  thus,  providing  the 
basis  for  resolving  the  issues  enumerated  above. 
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An  Input-Output  View  of  Robustness  in  Adaptive 

Control* 
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An  input-output  theory  of  adaptive  control  provides  a  means  of  determining  the 
robustness  properties  of  adaptive  algorithms. 
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Abstract— The  stability  and  robustness  properties  of  adaptive 
control  systems  are  examined  using  input -output  stability 
theory,  i.e.  passivity  and  small-gain  theory.  A  generic  adaptive 
error  system  is  developed  based  on  the  concept  of  a  tuned 
system— an  ideal  converged  (nonadaptive)  closed-loop  system. 
Using  this  error  system  with  passivity  theory  gives  conditions  for 
global  stability  where  only  boundedness  (in  norm)  is  required  on 
the  external  inputs,  eg.  disturbance,  reference  and  initial 
conditions.  Small  gain  theory  is  used  to  develop  local  stability 
results  where  the  magnitudes  of  the  external  inputs  are  restricted. 
In  the  global  case,  a  particular  system  operator  (not  the  plant)  is 
required  to  be  strictly-passive,  a  condition  which  is  unlikely  to 
hold  in  actual  use  due  to  unmodeled  dynamics.  The  local  results, 
however,  are  not  so  restricted  and  allow  for  unmodeled 
dynamics.  In  this  latter  case  an  estimate  of  the  stability  margin  is 
given  under  a  persistent  excitation  condition. 

I  INTRODUCTION 

At  a  very  basic  level,  the  issues  involved  in  adaptive 
control  design  are  no  different  from  nonadaptive 
(robust)  control  design.  In  either  case  the  goal  is  to 
maintain  specified  performance  properties  despite 
uncertainty  about  the  dynamics  of  the  plant  to  be 
controlled,  as  well  as  uncertainty  about  its 
environment.  In  the  nonadaptive  case  the  problem 
of  robustness  to  unmodeled  dynamics  is  well 
formulated  (e.g.  Doyle  and  Stein,  1981 ;  Zames  and 
Francis.  1983).  However,  research  in  adaptive 
control  theory  has  focused  almost  exclusively  on  the 
case  where  the  plant  can  be  fully  represented  by 
some  member  of  a  family  of  linear  finite¬ 
dimensional  parametric  models  (e  g.  Narendra,  Lin 
and  Valavani,  1980;  Goodwin,  Ramadge  and 
Caines.  1980).  Thus,  the  model  error  due  to 
unmodeled  dynamics  is  presumed  to  be  zero. 
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Unfortunately,  unmodeled  dynamics  can  cause 
adaptive  controllers  to  exhibit  significant  perfor¬ 
mance  degradation  and  instability,  even  with  an 
initial  controller  parameterization  that  closely 
approximates  the  desired  closed-loop  response 
(Rohrs  and  co-workers,  1981,  1982;  Ioannou  and 
Kokotovic,  1 983a, b ).  These  simulated  circumstances 
of  undesirable  behavior  are  in  sharp  contrast  with 
successful  applications  of  adaptive  control  where 
reduced-order  modeling  is  unavoidable  (e.g.  Astrom. 
1 983 ).  This  issue  of  model  error,  then,  is  of  undeniable 
practical  importance,  because  no  actual  plant  is  truly 
linear  and  finite-dimensional. 

Perhaps  the  main  reason  for  the  lack  of  a 
robust/adaptive  control  theory  is  that  the  emphasis 
has  been  on  global  results.  What  we  mean  by  ‘global' 
is  that  the  intent  is  to  require  as  little  a  priori 
information  about  the  plant  parametrization  and 
the  external  inputs  as  possible  to  prove  stable 
behavior.  Because  of  this,  the  resulting  requirements 
(i.e.  assumptions)  are  too  strong,  e.g.  known  plant 
order.  Therefore,  it  is  compelling  to  abandon  the 
requirement  of  global  stability — a  requirement  that 
may  well  be  beyond  the  needs  of  most  actual 
systems — and  develop  conditions  for  local  stability. 
The  term  local'  is  used  in  the  sense  that  the  plant 
uncertainty  and  external  inputs  are  limited  in  a 
defined  way,  e.g.  by  restricting  the  magnitude  and 
spectrum  of  the  reference  commands  and  disturb¬ 
ances,  as  well  as  the  initial  adaptive  parameter  error. 

In  this  paper  we  will  present  an  input-output 
view§  of  robustness  in  adaptive  control.  In 
particular,  we  shall  draw  attention  to  uncertain 
unmodeled  plant  dynamics — often  referred  to  as 
model  error — and  to  uncertain,  but  bounded, 
disturbances.  Based  on  this  view  it  may  be  possible 
to  merge  robust  control  theory  with  adaptive 
control  theory. 

The  next  section  (Section  2 1  formalizes  the 
conversion  of  a  generic  adaptive  controller  to  an 
equivalent  generic  error  system.  The  input  output 
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properties  of  the  error  system  relate  the  perform¬ 
ance  of  the  nominal  control  system  to  that  of  the 
adaptive  control  system.  Section  3  applies  this 
formulation  for  a  specific  continuous-time  adaptive 
model-following  algorithm.  This  permits  the  appli¬ 
cation  in  Section  4  of  stability  results  for  the 
continuous-time  version  of  the  generic  error  system. 
This  section  also  includes  a  discussion  of  the  strictly 
positive  real  (SPR)  condition  imposed  on  an 
operator  within  this  error  system.  Finally,  in  Section 
5,  we  will  examine  the  issues  involved  in  obtaining 
conditions  for  local  stability  and  robustness. 
Though  this  paper  concentrates  on  continuous-time 
systems  (due  to  space  limitations),  this  same 
input  -output  approach  is  applicable  to  robustness 
analaysis  of  discrete-time  adaptive  control  (e.g. 
Kosut.  Johnson  and  Anderson.  1983;  Ortega  and 
Landau.  1983)  as  well  as  time- varying  systems 
(Gomart  and  Caines,  1984). 

2  ERROR  SYSTEM  FORMULATION 
(A )  Tuned  control  concept 

Consider  the  nonadaptive  control  system  of  Fig. 
1,  described  by 

e  =  Hrw[n)w  (1) 

where  e|r)€R"*  is  the  error  output,  w(f)eR"'v  is  the 
external  input,  and  neR""  is  a  set  of  controller 
parameters  to  be  selected  For  our  purposes,  //„.(•) 
represents  a  closed-loop  parametric  feedback  system 
dependent  on  the  adjustable  parameters  in  n.  The 
output  e  of  //,„.(•)  is  the  error  the  control  system 
experiences  in  meeting  its  objective  given  the  external 
input  w.  Portions  of  //,.(•)  are  not  entirely  known, 
e  g.  the  open-loop  plant  imbedded  in  Hew{  ).  The 
input  w(f|  is  also  not  entirely  known  but  can  be 
assumed  to  be  in  a  subset  H  of  bounded  signals.  For 
example,  n  ( r )  can  consist  of  a  set  of  reference 
commands  and  bounded  disturbances.  If  the 
imbedded  controller  were  adaptive,  it  w  ould  adjust  it 
continuously  on-line  soas  to  reduce  the  error ;  but  for 
now  assume  that  n  is  constant  and  will  be  selected  off¬ 
line. 

If  the  control  designer  had  all  the  ‘off-line’  time  in 
the  world  to  ‘fiddle’  with  the  parameters  7t.  then  it  is 
hoped  that  a  satisfactory  adjustment  would  be 
obtained.  Many  strategies  can  be  envisioned  for 
determining  a  satisfactory  n.  In  fact,  such  a 
satisfactory  parameterization  may  not  be  unique 
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but,  rather,  be  any  member  of  a  set  S,  e.g. 

S  =  !7te/?n’,|//FH(7t)  has  desired  properties!  (2) 

Certain  parameter  sets  S  correspond  to  well- 
defined  design  strategies.  Specifically: 

Matched.  Let  S  denote  the  matched  parameter 
set,  i.e. 

5=  {neRn'\HeJn)  =  0\.  (3) 

Robust.  Let  S°  denote  the  robust  parameter  set, 
i.e. 

S°  =  {rt°eR""\\\Hew(n°)w\\/M  <  p°,  Vw()e  W\ 

(4a) 

where  the  norm  ||  j|  is  defined  on  the  underlying 
function  space.  The  finite  constant  p°  represents  the 
robust  performance  specification.  Note  that  S° 
includes  the  ‘optimal’  robust  solution,  i.e.  those 
n  e  R"’'  that  solve 

inf  sup  (||tf,„(n)w||/||w||).  (4b) 

It  H‘ 

Tuned.  Let  SJ, ,  denote  the  tuned  parameter  set 
associated  with  a  particular  w()e  IV,  i.e. 

Sin  =  {^  >€R“m»(*:..>MI/H  <  P* !•  (5a) 

The  finite  constant  p*  represents  a  tuned  perfor¬ 
mance  specification.  In  order  for  (4)  and  (5a)  to  be 
meaningful,  it  is  necessary  that 

P*  <P°  (5b) 

i.e.  the  desired  tuned  performance  is  better  than  the 
desired  robust  performance.  Also,  SJ,.,  will  include 
the  ‘optimal’  tuned  solution,  i.e.  for  each  we  W.  those 
ne  R"n  that  solve 

inf(||WfH.(jr)w||  j|w||).  (5c) 

* 

Ideally,  the  adaptive  control  should  converge  to  the 
optimal  parametrization  of  (5c).  Thus,  the  tuned 
parameter  set,  denoted  by  S*.  is  given  by 

s*=  U  (6) 

W<-»€  H’ 

Note  that  each  element  of  S*  is  satisfactory  for  a 
particular  w(  )e  H7and  that  no  one  element  in  the 
subset  Si,,,  c  S*  need  provide  satisfactory  control 
fora  different  w(  ).  (Although  nj,  ,eS* emphatically 
denotes  the  dependence  of  the  tuned  parameters  on 
w(  ).  we  will  henceforth  denote  membership  in  S*  by 
the  simpler  notation  n*eS*.  where  the  w(  ( 
dependence  is  to  be  implied.) 

The  error  signal  corresponding  to  the  matched 
case  is  identically  zero.  It  is  this  particular  case  that 
has  received  practically  all  the  attention  in  adaptive 
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outiiul  icsc.iich.  unit  about  wliiv.lt  the  sliongcsl 
tllCOICllC.il  I  C  s  1 1 1 1  s  .IK’  available  l  'llllll  t  llll.tlolv  .  Ill 
till’  lllsl  place.  till'  C.I'C  excludes  UlimeUMIIuble 
bounded  dislui  bailees  which  arc  ;i  viitu.il  ceitamty 
in  any  actual  system  Ik  immeastiiable  bounded 
ilisi  in  hanci’s  wc  mean  those  vlistni  bailees  which  can 
not  he  totally  leieeleil  at  the  oulptil  of  the  plant.  In 
the  sccoinl  place,  there  will  always  he  unimuleleil 
dynamics.  i  e.  theie  aie  never  enough  parameters  in 
n  to  solve  //,  „  t -  I  0  in  practice.  These  remarks 
apply  equally  in  a  stochastic  environment.  For 
example,  wheteas  in  the  ileterministic  case  e  =  0.  in 
the  stochastic  ease  /. | »•  |  0.  with  l.\  |  the 
expectation  operator.  Thus,  the  unmeasurable 
houmlcil  ihsturhanees  alluded  to  ahvue  have  their 
stochastic  counterpart  as  processes  which  do  not 
have  zero  mean.  i.e.  /:Jc!  *  0  for  any  n. 

The  more  appealing  of  the  other  two  sets  is  the 
tuned  set  S*.  defined  in  (6).  The  associated  error 
signal 

e*  =  H,  „(7t*)u  (7) 

is  referred  to  as  the  tumil  error  and  W(7r*l  as  the 
tuned  system.  Although  e*(f )  =  0  is  ruled  out  due  to 
the  impracticably  of  n*  e  S.  we  do  not  preclude  the 
case  where  c*(t)  —  0.  This  latter  case  still  presumes  a 
degree  of  idealization.  Consider  the  case  where  the 
external  input  n  consists  of  a  step  reference 
command  with  no  disturbance  and  e*  is  the 
difference  between  the  plant  output  and  the 
reference  command.  Thus.  e*«)  —  0  is  the  ideal 
output  error  for  am  stabilizing  controller  engender¬ 
ing  unit  d.c.  gain  This  class  of  tuned  controllers  can 
be  quite  large  even  if  dim  (rr*)<dim  (7t).  Now 
consider  the  impact  of  a  bounded  disturbance, 
which  is  not  necessarily  of  any  particular  functional 
form,  such  as  a  broadband  bounded  signal.  Clearly, 
with  such  bounded  disturbances  present.  e*(t)  -r*  0, 
and  can  only  be  assumed  to  be  bounded. 

An  important  comparison  for  the  tuned  set  S*  is 
to  the  robust  set  S"  (4).  Let 

e"  =  W.JttV  <«> 

denote  the  robust  error  Recall  from  (5)  and  (6)  that 
the  tuned  parameters  n*  arc  dependent  on  a 
particular  u  |  )  e  It  ,  whereas  the  robust  parameters 
n"  are  not.  Hence,  the  tuned  error  can  never  exceed 
the  robust  error,  i.e  for  a  particular  »»(  )e  W, 

lid  -  l|//,«(rt* |w’||  <  ||i°|l  =  Il//„(n0|w||.  (9) 

Condition  (9) also  follows  from  the  fact  that  p*  <  p° 
(5b|  Note  that  it  is  possible  for  the  robust  set  S°  to 
be  empty  even  though  the  tuned  set  S*  is  not.  If  S"  is 
not  empty,  then  consideration  of  an  adaptive 
controller  is  justified  if  for  some  ’large'  subset  of 


»it-|i  It.  v. moils  t inn’ll  coiitiolleis  exist  such  that 

■  ill  It  CllgCIldcis 

11*  *1!  v  llf"||.  (I0| 

II  this  wete  not  the  case,  then  a  robust  conltollci 
would  suffice  I  Ins  icquiicmcul  ( 10)  is  weaker  than 
the  lequiiemenl  />*  v  /•".  winch  may  not  be 
attainable  lot  all  nth  II .  since  ( 10)  is  required  only 
ovei  a  subset  ol  It  However,  even  if  (10)  holds, 
adaptation  may  cause  the  error  during  adaptation 
to  become  eithei  excessive  or  to  otherwise  exceed 
specifications 

The  usefulness  of detimng  the  tuned  parameter  set 
will  be  borne  out  in  the  next  subsection  The  tuned 
set  is  used  there  to  develop  a  generic  adaptive  error 
system.  At  this  point,  however,  we  rema*1  hat  it  is 
not  necessary  to  solve  the  optimization  problem 
defined  implicitly  in  (5).  rather  wc  only  need  to 
know  that  a  solution  exists  which  is  better  than  the 
robust  solution  (4). 

(B)  Adapt  ire  error  system 

Now  consider  the  adaptive  version  of  ( I ). 
depicted  in  Fig.  2.  and  described  by  the 
input-  output  relations 

iiiai 

n  =  fi[rf<0).  c.  i]  (lib) 

where  n|r  )e  R"~  arc  the  adaptation  parameters  which 
are  generated  from  the  parameter  adaptive  algorithm 
fi,  and  n(0)e  R is  the  initial  parameter  estimate. 
The  adaptive  algorithm  is  driven  by  the  output  or 
adaptationerrore{t)€  R"’  and  the  regressor  £{t)e  R"1. 
The  regressor  is  obtained  from  sensed  signals  w  ithin 
the  feedback  system. 

We  want  to  ultimately  determine  conditions 
under  which  the  adaptive  system  (II)  is  stably 
attracted  to  the  set  of  tuned  systems  (6).  Recall  that 
the  tuned  system  set  is  likely  to  contain  more  than  a 
single  member,  thus  by  stabili.y  wc  mean  stability 
‘about’  a  (possibly  disconnected)  set  rather  than 
about  a  point. 


»<•) 


I  ii.  2.  Adiiptivr  »)Mcm 
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I  Ik-  .iii.iKmv  is  l.iciht.iU'd  In  tiaiisloiming  (lie 
lniDi.il  i>l  I  ig  ?  min  .hi  i  mix  \|  sfciii  lnim.it  I  o  do 
tills  we  must  dell  IK'  ills'  si  OKI  Ills'  «>l  ills'  ildliptlVC 
null) i‘l  (  nnsids'i  .1  single-input  single-output 
(SISO)  plant  imbedded  m //(  )  whose  input  u(f) 
is  given  by  ills'  I’llnnar  s'Spis'ssion 

all)  .,(/)(}(/).  (12) 

Note  dial  this  development  is  nol  limited  to  SISO 
plants  Hie  extension  of  (12)  to  the  multivariable 
ease  involves  a  similar  expression  lor  eaeli  control 
ehannel.  i.e. 

ii.l/l  =  -  »,(/)  r,l/».  i  -  I . »/„  (13) 


where  2,(/|  and  r,(/l  are  vectors  consisting  of 
elements  from  the  regressot  and  parameter  vectors, 
respectively.  However,  only  the  SISO  case  will  be 
considered  here  to  reduce  the  complexity  of  the 
development  and  allow  sharper  focus  on  the 
adaptive  systems  issues.  Normally,  c(/)  consists  of 
the  plant  inputs  and  outputs,  or  filtered  versions 
thereof.  For  example,  in  discrete-time  systems  2(f) 
consists  of  a  finite  record  of  past  plant  inputs  and 
outputs. 

Although  the  bilinear  structure  in  (12)  and  (13) 
remains  the  most  widely  used  and  studied  format, 
nonetheless,  other  structures  (as  yet  underdevel¬ 
oped)  may  be  more  suitable  to  certain  problems 
c.g.  distributed  and  or  nonlinear  structures. 

We  will  now  make  a  strong  assumption  regarding 
the  way  in  which  i<( r )  and  tv(r)  arc  transmitted 
through  //(.)  into  c(f)  and  2(f). 

A  ssumption.  The  map  (it ,  w )  [  — *  (t\  2 )  is  linear  time- 
invarianl  (LTI).  i.e. 

r<M/n  =  r g..j.m  G..„(vf|r H(f)i  =  r h(i)"| 
|_2(f  )J  Lg,  J.M  G;J.s)JLm(/)J  Lu<fU 


where  G(.s)  is  the  open-loop  interconnection  matrix 
whose  elements  are  proper  rational  functions.  (To 
simplify  notation  we  will  use  v  to  denote  either  the 
Laplace  transform  variable  or  the  differential 
operator,  depending  on  the  context.) 

The  adaptive  system  (II)  with  bilinear  control 
(1 3) and  LTI  interconnections  ( 14)  is  shown  in  Fig. 
3.  To  transform  this  system  to  an  error  system, 
define  I  he  parameter  error 

n(f )  =  Mn  -n*  (15a) 


n*  e  S* 


6u)  w  (o) — -{a) 

((n  T 
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and  the  adaptive  control  error 


r(/|  =  c(f)  fi(/|.  (IM 

Thus,  Fig.  3  can  be  redrawn  as  show  n  in  Fig.  4  and 
described  by 

M- -c«4_ 0 ,.l*«4  M,ll.  ,i7„, 

Ls(f)J  L  -H'l  J 


Hence. 


poiy/fM,  «?,(.')] r  w(fn 

Inn J  [Hfjs)  //?.(.v)JL-i’(/iJ 

-"•“'[-rt.'i]  ,,8ai 

where 

HIJs)  =  G„(s)  +  G„(.v  )7r*'(/  +  G.B(.s)7i*)' 1  G.».(.v ) 

(18b) 

W*.(.v)  =  Grv(.v )  +  G,„(v)7r*'(/  +  G.M(v)7T*r  'G<B(s) 

(18c) 

H?M)  =  (f  +  G.,,(x)n*  I  'G.ds)  (I8d| 

/#*,(»')  =  (/  +  G<y(s)Jr*  )-  'G^s).  (I8e> 

The  dashed  box  in  Fig.  4  is  W*(.s).  We  will  refer  to 
H(s)  as  the  tuned  interconnections.  Note  that  the 
tuned  error  (7)  is  identical  to 

e*(f)=  H;js)  w(t).  (19) 

We  also  make  use  of  the  tuned  regressor,  defined  as 

{•(f)=  #/<*„.(%)  w(f).  (20) 

Finally,  the  error  system  (Fig.  4)  can  be  depicted  as 
in  Fig.  5,  where 

«•(»)  =  <•*(»)  —  H*(v)'(H  (21a) 

4(f)  =  {*(0- HJ(M«’(0  (21b) 

i(f)  =  v(M’n(0  (21c) 

n(r)  =  C2 (n<(>),  «•(■).  {(•)].  (2ld) 


Input  output  view  of  adaptive  contiol 


I - 1 


1 i<i  4  Adaptive  error  system 


Figure  5  reveals  that  this  error  system  is  composed 
of  a  nonlinear  system  in  the  forward  loop,  denoted 
by  and  the  LTI  system  H*  (s)  in  the  feedback 
path.  Thus,  the  error  system  is  driven  externally  by 
the  tuned  system  outputs  e*(f)  and  z*U).  and  the 
initial  parameter  error  jt(0)  =  jr(0)  —  7r*. 

(C)  Existence  of  the  timed  controller 

The  designation  of  the  tuned  controller  is  the 
concept  most  important  to  extracting  a  meaningful 
error  system  from  the  description  of  an  adaptively 
controlled  system.  It  might  appear  that  the  ability  to 
specify  this  tuned  controller  presupposes  our 
knowledge  of  an  acceptable  solution  to  the 
underlying  adaptive  control  problem.  This  is  not 
entirely  the  case.  Given  the  parametric  controller 
structure  of  C().  we  need  only  have  an  approximate 
a  priori  knowledge  of  the  system  behavior.  Given  a 
particular  7t*,  we  will  discover  that  the  restrictions 
on  H*.  and  H*  can  be  assessed  from  knowledge  of 
the  tuned  controller  and  bounds  on  the  magnitude 
of  the  plant  modeling  error.  Such  information  is  a 
practical  result  of  a  thorough  plant  modeling  study. 
Thus,  the  study  of  the  stability  of  (21)  will  have 


practical  meaningfulness.  We  will  examine  a 
particular  continuous-time  adaptive  controller  in 
the  following  section  and  derive  the  form  of  //*  and 

H* 

V  CONTI NUOUS-TIMF  ADAPTI  Vfc  MODLL- 
FOLLOWING 

To  characterize  H*  and  //.*,  some  designation  of 
a  tuned  controller  must  be  provided.  We  make  the 
choice  by  assuming  that  no  modeling  error  exists  in 
the  nominal  plant  parametric  model.  We  close  with 
consideration  of  the  degree  of  plant  mismodeling 
allowed  such  that  this  tuned  controller  is  robust,  i.e. 
maintains  stable  control  of  the  actual  system. 
Following  this  discussion,  in  Section  4  we  consider 
the  effect  of  the  adaptive  algorithm. 

(A )  Direct  mode!  reference  adaptive  control 
Consider  the  model  reference  adaptive  control 
(MRAC)  system  shown  in  Fig.  6.  described  by 

r(i)  =  </(f)  +  P{s )u{t )  (plant)  (22a) 

y(r)  =  Hls)r(t)  (reference  model )  (22b) 

e(f)  =  ,v(t)  —  >•(/)  (tracking  error)  (22c) 


I  k.  V  Adaptive  error  system 
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Fig.  6  Model  reference  adaptive  system  (a)  block  diagram,  (b)  controller  detail. 


where  d{t)  consists  of  disturbances  and  plant  initial 
conditions,  and  r(t)  is  the  reference  command.  Let 
{C(),Q[  denote  the  adaptive  controller,  where  (1 
is  the  parameter  adaptive  algorithm  and  C(  )  is  the 
parametric  controller.  Following  Narendra,  Lin, 
and  Valavani  (1980).  let  C(  )  have  the  bilinear  form 

M0  =  -{(f)'A(i) 

=  -£.(0'n.(0  -  C>(0'»U0  -  {,(0'(M0  (23a) 

where  the  regressor  is  given  by  filtered  versions  of  u, 
y  and  r 

{(O'  -  [{.(O'.  {>(0\  {,(0]' 

=  [F(s)M0.  F(s)y(t ),  —  F(s)r(0]  (23b) 


F(s)  = 


J_  til 

Us) . Us) 


Us)  *  s*  +  a,  s*  ‘4...  4  a*.  (23d) 

Thus,  there  arc  3 k  adaptive  parameters.  Using  the 
definition  of  adaptive  control  error  in  (21c),  the 
MRAC  control  signal  (22)  can  be  expressed  as 

MO  =  --J^MO-  ^>0  4  4jy00  -  i>(0 
Us)  Us)  Us) 


where  the  tuned  parametrization  n*(  =  A  —  a)  is 
distributed  among  the  control  elements  as  follows: 

>l*(s)  =  n?s*  ~ 1  4  ...  4  tiJ  (25a) 

^J(s)  =  7tf+,s* (25b) 
A$(s)  =  n^+iS*-1  4  ...  4  rtf*.  (25c) 

Thus,  (24)  becomes 

MO  =  C*(s)  r(f )  -  C*.(s)  y(f)  -  C*  (s)r(r)  (26a) 

where 

*!(0 


C*  (s)  = 


Us)  4  /4*(s)' 


"r  1  us)  +  xrcv)’ 


c:,.(.v)  = 


Us)  4  /l?(.v>- 


We  will  refer  to  C*  =  [C*„  CJr  C*  ]  as  the  tuned 
controller.  The  adaptive  error  system  (21)  cor¬ 
responding  to  the  M  R  AC  of  Fig.  6  is  shown  in  Fig.  5. 
The  tuned  signals  (3.19)  and  (3.20)  are 

e*  m  [<|  +  PC*,)-  'PCI -R)r-t  I  4  PC*)  1 d 

(27a) 

(I  +  PC*  r  ,C*Fr  -  C*  (I  +  PC;,)’  lF  d 
{*  =  (14  PC;,)-  lPC*,Fr  -  (I  4  PC*,)-  xF'i  . 
-Fr 


(24) 


30 


(27b) 
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ii*.  !ii  t  /•<  ;,(  '/•(;,/ 


|/f|  I  until  m  Ait’/n  nmtrol  design 

I  here  arc  any  number  of  ways  to  design  the  tuned 
controllci  ('*.  The  important  point  no  mutter 
what  design  method  is  used  is  that  the  tuned 
design  must  be  robust.  because  the  plant  P(s|  in  (27) 
and  |2S)  is  not  entirely  known.  Recall  from  (5)  that 
the  tuned  controller  is  dependent  on  the  plant.  For 
example,  the  3k  parameters  in  n*  cannot  make 
r  |  — •  e  in  (27a  I  be  identically  zero.  This  can  be  viewed 
as  a  reduced  order  design  problem  or.  as  in  the 
discussion  that  follows,  a  problem  in  robustness  to 
unmodelcd  dynamics. 

Suppose  that  the  actual  plant  can  be  described  by 
Pis  I  =  [I  +  A(.v)]P*(.v)  (29a) 

■4*(s| 


MV +  /».*- 1  +-..  +  HJ 
S "  +  «,s"'  '  +  ...  +  <!„ 


,  m  <  n 


where  P*(x)  is  a  tuned  parametric  model  of  P(s),  i.e. 

the  parameters  (h, . h„.  a, . a„ (  provide  a  good 

fit,  say  at  low  frequencies.  The  transfer  function  A(s) 
represents  unmodeled  dynamics,  i.e.  those  dynamics 
in  P(.s)  not  accounted  for  by  P*(x),  e.g.  high 
frequency  efforts.  Assume  that  A(s)  is  stable  but  is 
otherwise  unknown  except  for  a  bound,  i.e. 

|A(j<u)|  <  £(<»).  Vwc  R.  (29c) 

This  type  of  modeling  uncertainty  is  said  to  be 
unstructured  (Doyle  and  Stein,  1981).  In  more 
general  terms,  (29)  provides  a  set  description  of  the 
plant  rather  than  a  single  parametric  model,  such  as 
P*  (Safonov,  1980). 

We  will  now  examine  the  impact  of  model  error 
on  a  tuned  control  design  based  only  on  the 
parametric  model.  The  model  reference  format 
suggests  that  we  make  e *  as  small  as  possible.  To 
eliminate  the  tracking  error  term  in  (29a)  entirely, 
we  will  use  the  procedure  described  in  Egardt 
(1979),  which  requires  that  the  following  infor- 
mation  is  known: 

(1 )  n  >  m  (P*(s)  is  strictly  proper) 

(2)  n  and  m  arc  known 

(3)  £*(>)  has  all  zeros  strictly  inside  the  left  half 
plane. 


Also,  the  icleicncc  model  ti.mslei  lunclion  j*> 
assumed  to  be 

m  i  y  t, 

•t(s)  s  4  u,s  -4  ...  4  1 1„ 

wheie  />,, ....  <i , . u„  are  pieselected  constants. 

and  wheie  //( s|  is  exponentially  stable,  i.e.  all  zeros 
of  l(s)are  strictly  inside  the  left  half  plane.  It  is  well 
to  point  out  here  that  although  assumption  (2) 
above  can  be  satisfied  by  the  parametric  model 
(29b).  this  is  not  the  case  for  the  actual  plant  (29a) 
due  to  the  presence  iff  the  unstructured  uncertainty 
(29c) 

The  tuned  controller  structure  proposed  in 
Egardt  (1979)  requires  that 

'  Sr  ,3,a) 

RT* 

C-  ‘  KBS- 

where  T*  is  a  stable  monic  poly  nomial  of  degree 
n,  >  n  —  m  -  1,  and  where  the  polynomials  S*  and 
R *  uniquely  solve  the  polynomial  equation 

T*A  =  R*  +  /IS*  (31c) 

with  S*  monic  of  degree  nr.  and  R*  of  degree  n  —  I. 
With  no  model  error  (A  =  0).  this  controller  (31).  in 
addition  to  stabilizing  the  tuned  system,  also  makes 
the  transfer  function  from  r  into  y*  identical  to  the 
reference  model  H(s).  Thus,  the  tuning  of  (31 )  is  for 
the  subset  of  W  composed  of  bounded  reference 
signals  and  zero  disturbance.  The  effect  of  (3 1  )on  C* 
will  shortly  be  made  apparent.  Comparing  (33)  with 
(26)  motivates  solving  for  7t*  from 

L  +  AT  =  B*S* 

A1  =  ~R*  (32) 


A%  =  -  BT*. 
ho 

A  solution  for  n*  exists  provided  that 

k  =  nT  4  m  S:  n.  (33) 

With  this  choice  for  (/4f.  A $,  /4*>.  the  tuned 
controller  is  given  by  (31 )  and  by 

C*  =  L  (34) 

BS* 

(C)  The  effect  of  model  error  on  tuned  system 
performance 

It  is  convenient  to  define  the  transfer  function 
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(I  4  AG*  I 


(.*!//»  ( 


<3M 


(I  .  Ac.*  I 
(I  .  A<>* I 


t*  Ji 

1  *R*  , 

.  1  r 

-  /  d 

/>,.»*  1 

hjl*  7  *  4 

a 

N*  4*  , 

(1  4  A) 

r  -  1  d 

A 

7  *  4  I 

-  I  V 


(36b) 


and  the  tuned  interconnections  (28)  become 


4  (.1  ohm  si  Min  m  < oMinioNs 
I  lie  pm  pose  ol  this  section  is  In  ml  induce  global 
stability  conditions  applicable  to  the  generic  eiror 
system  of  (21 1.  In  the  preceding  section,  vie  specified 
an  adaptive  coni  toller  structure  (  '(.)  from  which  we 
then  developed  the  tuned  system  (r. </||— • 
and  the  interconnection  operators  //*,  and  H*.  We 
now  need  to  characterize  the  adaptive  law  fi  in 
(2 Id).  With  this  connection  we  will  be  able  to 
interpret  some  conditions  under  which  such  a 
continuous-time  adaptive  controller  possesses  a 
(limited)  degree  of  robustness.  Our  interpretive 
remarks  will  address  the  rest rictivc ness  of  the  SPR 
condition  on  II*,  that  arises  in  practically  all  global 
stability  theorems. 


//*  =  (I  +  AG*)  '(I  -t  A) 


h0L 

7*4 


(37a) 


•  A*  l 

(1-4  AG*)  '  --F 

I  B*T*A 

H*  =!  (i  +  AG*r'<i  +  a)~f 
I  T*  A 


•  (37b) 


The  tuned  system  with  no  model  error  (A  =  0)  is 
exponentially  stable,  since,  by  assumption,  the  poles 
of(fl*)~  '.(.4)*  ‘.and  (T*)' 1  are  in  the  open  left  half 
plane.  Hence  e*  and  ci*  are  bounded  if  r  and  d  are 
bounded.  Thus,  the  stability  of  the  actual  tuned 
system  is  guaranteed  if  and  only  if 


(1  +  AGT  1  and  (1  +  AG*)* ‘A 

are  exponentially  stable.  (38) 


Note  that  under  these  conditions,  the  tuned 
interconnections  in  (37b)  remain  exponentially 
stable.  However,  it  is  not  necessary  (nor  possible  by 
assumption)  to  have  a  complete  description  of  A  in 
order  to  satisfy  (38).  For  example,  if  A  is  known  to  be 
exponentially  stable,  then  with  G*  known  to  be 
exponentially  stable,  (38)  holds  if  (e.g.  Doyle  and 
Stein.  1981) 

|A(j<.;)|-|G*(jm)|  <  1.  Vuu.R.  (39) 


Satisfaction  of  (39)  requires  that 


(4)  The  adaptive  algorithm 
We  w  ill  begin  by  specifying  the  adaptive  law(sl  of 
interest.  A  large  class  of  adaptive  algorithms  (2 Id) 
have  the  from 

f.  (r)  =  A  [.;(•).  <u(/)].  n(0|e  Rr  (41a) 
c >(t)  =  2(f)  e(r).  (41b) 

We  will  refer  to  /f  (..  .)as  the  adaptation  gain,  w  hich  is 
a  nonlinear  operator.  In  general  A[\ •]  can  have 
memory,  usually  only  in  ^(j).  The  adaptive 
algorithm  can  also  be  expressed  in  terms  of  the 
parameter  error  n(t)  as 

A(t)~  A  [;(•).  <«(!)].  n(0)  =  rt(0)  -  it*.  (42) 

The  complete  adaptive  error  system  (21 ),  including 
the  adaptive  algorithm  (42),  is  shown  in  Fig.  7. 

The  choice  of  algorithms,  i.e.  the  variety  of 
proposed  adaptation  gains,  is  virtually  unlimited. 
The  following  two  are  our  chosen  representatives: 

Constant  gain  (Narendra.  Lin.  and  Valavani,  1980). 

/4[£(),w(0]  =  /lor->(f) 

(43) 

w here  A„  e  Rr ' r.  A„  =  A'„  >  0. 

Retarded  gain  (Kreisselmeier  and  Narendra.  1982). 

<4  [£(■),  oil)] 

_  f  A{M»,  |(t(f  )l  <  c 

”U,[w(r)  -  (l  -  |n(()|/c)27t(f)].  |7t(f)|  ;>  r 

(44) 


|A(j>’i)|  <  fHo,)  =  l/|G*(jo»)|.  Vwcl|,  (40) 

We  will  show  in  Section  4  that  <>(<")  <  I  is  the  limit 
imposed  on  Mo>)  by  the  usual  global  stability  results 
for  continuous-time  adaptive  systems.  Similar  limits 
arc  also  encountered  with  discrete-time  adaptive 
systems. 


where  A„e  R1"  r.  Av  =  A'{,  >  0,  and  <  £  max  |n*|. 

Wc  will  use  the  concept  of  persistent  excitation 
that  has  proven  important  in  adaptive  control,  as 
well  as  in  adaptive  system  identification. 

Definition  (Anderson.  1977).  A  function 
H  ):R<  —  R"  is  persistently  exciting.  denoted  by 
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I  £  PL.  if  there  exists  positive  constants  a,.  x:  and  a( 
such  that 


3 2I„>  |  i  (0  /(/)</i  >  a,/„.  V.\e  R*.  (45) 

.% 

We  will  discuss  the  implications  of  persistent 
excitation  on  global  stability  below,  as  well  as  in 
Section  5,  in  regard  to  local  stability. 

(B)  A  global  stability  theorem 
Theorem  1.  which  follows,  gives  conditions  for 
global  stability  of  the  adaptive  error  system  of  Fig.  7. 
The  term  ‘global’  refers  to  the  intention  of  seeking 
the  minimal  (reasonable)  restrictions  on  the  tuned 
signals  c*(f)  and  c*(r),  and  the  tuned  interconnec¬ 
tions  H*  (\)  and  H*{s)  resulting  in  the  proof  that  e 
and  4  remain  bounded,  i.e.  (2 1 )  is  stable,  for  any  finite 
tt(0).  (A  detailed  proof  of  Theorem  I  isgivenin  Kosut 
and  Friedlandcr  (1983).)  In  particular,  we  will 
consider  the  following  two  tuned  system  signal  sets  as 
‘inputs’  to  the  error  system: 

m  =  !*’*,  £*.  n(0)|e*,c*  e  L2r\  L  , , 

2*  6//,.  71(0)6/?'’!  f4^) 

M  JS  =  !<-*,  «f*.  n(0)|c’*  e  L,  ,£*  e  Lry , 

n(0)eRr!.  (47) 

Note  that  <■*,  e*€L2n,L,  essentially  implies 
boundedness  and  ultimate  decay  to  zero,  whereas 
inclusion  in  L,  only  implies  boundedness. 

Theorem  I .  For  the  system  of  Fig.  7,  assume  that 

(Al)  //?,(>)  is  strictly  proper  and  exponentially 
stable  (48) 

(A2(  //*(.%)  is  strictly  positive  real  (SPR),  i.e. 
//•  (v )  is  strictly  proper,  exponent  tally  stable. 


and  there  exists  a  finite  constant  /»  >  It  such 
that 

Re  [//*  (jc>)j  >  i> |//*(j(D|’,  Vw6  [0.  j.  ).  (49) 

Under  these  conditions,  algorithms  (43)  or  (44) 
result  in  the  following  properties: 

(i)  If  |e*.  2*.  £(()))€  IF*  then.  n.  e.  c.  and  v 
are  bounded  (in  L,).  zr  ( r )  — ►  0.  and 
<’*(/)  -  t'(r )  — *  0.  In  addition,  if  2*  e  PE.  then 
7r(f)  -» 0  exponentially. 

(ii )  If  /<>*.  2*.  ?r(0))e  H  J  and  2  e  PE.  then  n.  e.  2. 
and  r  are  bounded  (in  L, ). 

(iii)  If  (c*.  c*.  n(0))eH’5  and  2#  PE.  then  the 
results  of  (ii)  still  follow  by  using  the  gain 
algorithm  (44). 

Comments  on  Theorem  I.  Though  theoretically 
significant,  these  results  do  not  offer  the  engineering 
design  guidelines  we  would  like  to  obtain.  The 
major  reason  is  that  W*  (.s)eSPR  (condition  (,42)) 
is  virtually  impossible  to  achieve  for  any  actual 
system.  The  primary  culprit  here  is  the  effect  of 
unmodeled  dynamics.  Details  on  this  issue  may  be 
found  in  Rohrs  and  co-workers  (1982).  Further 
discussion  will  be  provided  in  the  following 
subsection. 

Another  technical  hurdle  is  that  the  only  realistic 
case,  insofar  as  the  tuned  signals  (t'*,2*)  are 
concerned,  is  when  (<**. 2*)e  M'J.  This  is  the 
situation  induced  by  continual  bounded  disturb¬ 
ances.  such  as  would  normally  be  encountered.  But 
in  this  case  the  theory  requires  that  cither  2  e  PE  as  in 
part  (ii)or  that  the  adaptation  gain  is  retarded  as  in 
part  (iii ).  With  bounded  disturbances  present  it  is  not 
known  how  to  guarantee  2  e  PE.  since  2  is  generated 
inside  the  adaptive  loop.  Note  that  part  (i)  only 
requires  that  the  tuned  regressor  2*  e  PL  rather  than 
the  actual  regressor  2  e  PL  as  in  part  (ii).  However, 
this  requires  (<•*.  2  f  IF*,  w  hich  is  only  possible  when 
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tilt1  n-niiol  slimline  pi  ox  ides  ,ixy  mplnlu'  model 
1 1 '  I  It  <  \N  lilt’  .Mill  (llstlll  Kitin'  Il'Il'l'tlOH  lilts  Is  ills' 
ll.lxxli  I'.lxe  studied  III  till-  llU~l.il  Ills  Obviously, 
unmodeled  dy  it. minx  mid  hounded  dixiui  bailees 
slimui.its  this  iils.il  situation  A  iuilhei  dilliculiy 
legaiding  .  t  1*1  isth.it  this  osstti  s  .it  the  expense  of 
;in>  xet-poml  t spill. tilt’ll,  wlusli  delenoialex  in  the 
pi iVssikv  t>l'  1*1  signals  l  sing  gam  letaidatmn  does 
not  isi|ims  persistent  excitation.  hnl  iloss  require 
some  u  prion  inlt'i in. Hum.  is.  as  in  (44 1.  the 
foicknow  ledge  of  an  uppei  hound  on  |n*|.  which  is 
no!  toodillieult  loohtain  Although  retardation  docs 
handle  hounded  ihsiurhanee.  llie  SI’R  eondition  is 
still  requited 

1(1  hi  pursuit  <>l  ilu‘  SI’R  (  o  million 

II  Ill'll  ii  limn  points  lo  llh ■  slurs, 
mill  u  tool  It >o ks  ui  his  finger." 

Anonymous. 

The  mli”':  of  this  aphorism  is  to  divert  any 
lingering  anxieties  about  the  SPR  eondition.  It  —  the 
SPR  condition  simply  will  not  do  as  a  major 
building  block  in  adaptive  control  theory.  But  that 
does  not  mean  a  total  abandonment  of  our  aim;  it 
suggests,  rather,  a  redirection.  We  should  be 
establishing  a  different  path  to  the  ’stars.'  For  now, 
however,  we  will  remain  earthbound  and  address 
the  reslrictncness  of  the  SPR  requirement. 

A  necessary  condition  for  H*  e  SPR  is  that  H*  (s) 
have  a  relative  degree  of  one.  As  pointed  out  by 
Rohrs  and  co-workers  (1982).  this  imposes  the 
requirement  that  the  relative  degree  of  the  plant  is 
known,  eg.  examine  the  effect  on  the  plant  P  in 
(28a  |.  This  knowledge,  however,  is  unavailable  due 
to  the  presence  of  unmodeled  dynamics,  as  assumed 
in  (29). 

The  same  type  of  restriction  can  also  be  seen  as 
follows.  From  (37a) 

H*  =  ( I  +  AG* )'  1  ( I  +  A )/?  *  (50a) 

m.  -  £4  (5oi>i 

If  A  is  exponentially  stable,  but  is  otherwise 
unstructured,  then  conditions  for  //*.<:  SPR  include 
(Kosut  and  Fricdlander.  1983) 

(!)/?*€  SPR  (51a) 

(2)  |A(jm)|  <  I.  (51b) 

Since  //*  is  dependent  only  on  the  parametric 

model  /’*.  it  is  not  difficult  to  find  n*  such  that 
ft*  f  SPR  Unfortunately,  the  drawback  is  that 
(51b)  is  a  condition  that  is  almost  surely  violated, 
due  to  typically  unmodclcd  high  frequency  dy¬ 
namics 


So  il  //*,  -  SPR  i .in  nevei  hold,  i .in  we  I'hmm.ilc 
llie  SPR  i ii 1 1 1 1 1 i  n ii  i 1 1  oi  . 1 1 1 1 1  some  ili'iii  lilii'imp  lo 
dcsuahlv  .illi'i  //,*,  ’  I  oi  ilk'  pcilcii  modeling  cave 
(A  Id  it  is  possible  lo  obi. mi  iMouopoh.  1  *>74 ; 

1  . iiul. III.  I *>7S .  I  p. mil.  I *> “**> | 

//.*,  ii  ,,  //,*,  positive  constant  (52) 

Alihouph  a  positive  constant  is  SPR.  and  lienee, 
satisfies  (51).  condition  (51b)  is  still  icquncil  lor 
//•(A  *  0)  lo  be  SPR 

These  disclaimers  lead  iix  away  horn  the  global 
approach  ty pitied  by  T  heorem  I  lo  the  establish¬ 
ment  of  local  stability  results  which  aie  robust  to 
unmodclcd  dynamics  and  bounded  disturbances. 

S  I  (K\l  SI  Alii  I  Id  (OMM1IONS 

In  this  section  we  indicate  a  means  of  obtaining 
local  stability  conditions  To  clarify  the  distinction 
between  local  and  global,  consider,  for  example, 
result  (ii)  of  Theorem  I.  This  result  holds  if 
H*  e  SPR.  //*  exponentially  stable,  (e*.  2*  )e  ITj{. 

2  e  PE.  and  |n(0)|  <  /. .  Aside  from  the  difficulties  in 
establishing  SPR  and  PE.  all  the  conditions  are 
virtually  free  of  any  magnitude  constraints,  and 
hence,  arc  ‘global-  conditions.  In  every  practical 
case,  it  is  more  than  likely  that  magnitude 
information  is  available,  e  g.  a  priori  bounds  on 
||c* || , .  ||2*|| ,  and  |jf(())|.  as  well  as  a  bound  on  the 
gains  of  H*  and  H «*.  For  example.  Egardt  (1979) 
shows  robustness  properties  for  minimum  phase 
systems  with  bounded  output  disturbances.  Dead- 
zone  and  projection  mechanisms  can  handle  small 
unmodeled  dy  namics  as  shown  by  Praly  (1983)  and 
Samson  (1983).  Ioannou  and  Kokotovic  (I983a.b) 
are  able  to  give  an  estimate  of  the  region  of  attrac¬ 
tion  w  ithout  SPR  or  PE  in  the  case  of  high  frequency 
parasitics.  Persistent  excitation,  and  the  resulting 
exponential  stability  property  (sec  equation  (62)  in 
this  section)  also  leads  to  robustness  (e  g.  Anderson 
and  Johnson,  1982a. h;  Anderson  and  Johnstone, 
1983;  Kosut.  1983).  Various  other  gain  normali¬ 
zations  have  also  been  suggested  (e.g.  Gawthrop 
and  Lim.  1982;  Ortega  and  Landau.  1983).  These 
theoretical  results  remain  incomplete,  because  they 
do  not  as  yet  provide  a  useful  means  of  assessing  the 
impact  of  unmodeled  dynamics,  e  g  a  frequency 
domain  hound  on  model  ci-or.  dependent  on  the 
Tcturn-difTcrcncc  gain'  (e  g  Doyle  and  Stein.  1981 ). 

In  this  section  wc  will  show  in  Theorem  2.  under 
mild  magnitude  hounds,  that  the  adaptive  system  is 
(locally)  I- ,  -stable.  This  result  is  quite  general 
because  the  conditions  arc  independent  of  the 
nature  of  the  adaptive  algorithm,  e  g.  dead-zones, 
normalizations,  or  persistent  excitation. 

To  facilitate  the  analysis  wc  will  only  consider  the 
continuous-time  error  system  (21)  with  constant 
gain  adaptation  alorilhm  (43).  It  is  convenient  to 
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tianslonn  (21 1  to  t lie  lollowmg  vaiiuiioiinl  form, 
wind)  is  more  useful  for  local  analysis: 

'  'i  Vs/  (53a) 

<s/  -  /  /('I  (53b) 

where 

v  =  (n.  e,  el  =  In  -  n*.  «■  -  <•*.  i  -  c*(  (53c) 

v,  =  (Ttj.  i'|.  «i|  I.  /(v)  =  (4'n,  tc).  (53d) 

Details  on  transforming  (21 )  to  (53)  are  in  Kosut 
(19X3).  This  form  of  the  adaptive  error  system  is 
obtained  by  linearization  of  (21 1  about  <•*.  c*  and 
n*.  The  linearized  perturbation  response  is  v,., 
almost  identical  to  the  linearized  system  studied  by 
Rohrs  and  co-workers  ( 19X1 ).  which  was  arrived  at 
by  a  ‘final  approach  analysis.*  The  remaining 
nonlinear  terms  xsl  arc  contained  in  /(x),  a 
mcmorylcss  nonlinearity,  and  in  F,  a  time-varying 
linear  operator.  The  characteristics  of  F,  as  well  as 
those  of  depend  on  the  adaptation  gain  and  the 
behavior  of  the  tuned  signals,  e*  and  c*.  For 
example,  with  the  constant  gain  algorithm  (43),  the 
linearized  perturbation  response  is 

7tt  =  </  +  LM)~'n0  +  Ki*e*  (54a) 

t'i  =  (54b) 

(54c) 

with 

XV  -K 

f=  H*0-c*’KN)  (55) 

H*A\  -  i*'KN)  H&fK 

and  where 

L=  -A0  (56a) 

s 

K  +  LM)~'L  (56) 

M  »  (56c) 

N  =  <!*//*  +  c*H?r.  (56d) 

Since  boundedness  of  (t’*,<J*)  and  stability  of 

(H*  .H*.)  are  established  by  definition  of  the  tuned 
system,  it  is  not  difficult  to  sec  that  conditions  for  the 
stability  of  F  and  the  boundedness  of  xL  are 
identical.  In  fact,  this  follows  if  and  only  if  the  system 
S :  ( v0.  fV)|-»  x,  described  by 


Assuming  that  .Sr  /..-stable  we  obtain  the 
following  local  stability 


'Fheorcm  2 

Suppose  /  t  I. ,  -stable  and  \,  cl.,  Hence,  there 
exists  a  constant  c  such  that 

y ,  (/  )  <  <  <  /  .  (5X) 

Under  these  conditions,  if.  for  some  /;  <  2/c. 

||v,|| ,  S  (I  -  a  2)c  (59) 

then 

llv/.ll,  </:.  (60) 

Proof.  The  proof  is  entirely  analogous  to  the  proof 
of  the  linearization  theorem  on  p.  131  of  Desoer  and 
Vidyasagar  (1975).  Details  for  this  case  may  be 
found  in  Kosut  (19X3). 


Discussion 

Theorem  2  asserts  that  the  adaptive  system  is 
stable,  i.e.  bounded  inside  an  r-region.  provided  that 
FeL,  -stable  and  the  linearized  response  is 
bounded  and  sufficiently  small,  i.e.  condition  (59). 
No  claims  arc  made  about  the  mechanism  that 
provides  Fe  Ly -stable  and  ,x\eL, .  As  mentioned 
earlier,  these  are  insured  if  the  map  S  defined  in  (57) 
is  L ,  -stable.  It  is  possible,  of  course,  that  xLe  Ly  but 
||xt||T  exceeds  the  magnitude  constraint  of  (59). 
Instability,  however,  does  not  follow  because 
Theorem  2  only  provides  sufficient  conditions. 

In  order  for  theorem  2  to  be  of  practical  use.  it  is 
necessary  to  provide  stability  of  S  without  relying  on 
passivity  of  H*..  We  will  illustrate  this  by  using 
persistent  excitation.  Consider  the  system 

x  =  -  A  f  H  f  '.v  +  u,  \(0 leR"  (61) 

It  is  shown  in  Anderson  (1977)  that  if  Ae  R"’", 
/1  =  A'>0,  /ePE  and  WeSPR  then  (61)  is 
exponentially  stable,  i.e.  (here  exists  constants  m. 
A  >  0  such  that 

|x(/)|  <  mii*,|x(0)|  +  f  mt:  "|i»(r)|dr.  (62) 

Jo 


x  =  /4(i(w  -  Afx),  x(0)  =  x0  f‘Rn  (57) 

We  will  apply  (62)  to  provide  stability  of  S  as 
is  stable  (Kosut,  1983).  Note  that  the  system  S  is  follows.  The  system  S  can  be  written  as 
identical  in  form  to  the  linearized  parameter  error 

system  (jr„.4*r*)|-*  Bj.  in  (54).  35  x  =  -  A„lttrrl'x  4  A„*  -  Qx  (63a) 
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"lii'U'  , .  //  ,.  ami  tj  an'  iliTincil  via 


If.’, 

//,,  i  II.. 

(63b) 

* 

s  i  C 

(63c) 

■(,,<  M  J/..C  1 

(63d » 

Computing  K>3.i) 

with  (61).  intuitively. 

if ;  *  im  . 

//, ,  i  SI’R.  ami  (J  Millii'ii'iillv  small'  then  the  system 
(63 I  (equivalently  llii-  map  M  remains  exponentially 
stable.  I  Ini'.  by  T  heorem  2.  an  i  -region  iif  loeal 
stability  evist'  I  lie  precise  conditions  are  stated  as 
follows. 

(  omlUirv  2  I 

l  et  II  ■  SPR  and  e  c  1*1  with  corresponding 
positive  constants  /.  and  in  as  defined  in  (62 1.  Then. 
/  e  I. ,  -stable  and  x,  t  if 

i:si!,  ikli,  (2  •+  m,):  ah*) 

-t  lk*|Uk:*||,;'.lH*l  (64) 

and 

V,  1/7,., c )  <  (/.  m  -  r/l  lkli* .  (65) 

Prmtf.  Follows  directly  by  application  of  Small 
Gain  Theory  iZames,  1966)  to  (63).  Details  may  be 
found  in  Kosut  (193). 

Discussion 

Corollary  2.1  shows  that  persistent  excitation  is 
one  mechanism  which  provides  SeL?  -stable,  and 
hence,  boundedness  of  xL  and  stability  of  F. 
Therefore,  if  in  addition,  .v,.  is  sufficiently  small  (59) 
then  the  adaptive  system  has  a  local  stability. 

Other  mechanisms  to  provide  stability  of  S 
include  dead-zones,  retardation  functions,  and 
signal  normalizations.  Their  effect  on  S  needs  to  be 
determined. 

Corollary  2. 1  also  provides  an  upper  bound  on 
the  effect  of  model  error  via  (65).  This  is  not  yet  in 
the  frequency  domain  form  we  would  like,  but  the 
bound  can  be  quite  large.  Hence.  H*  need  not  be 
SPR.  but  only  approximately  so,  e.g.  H*  e  SPR. 
Think  of  H*  being  SPR  only  at  low  frequencies.  In 
the  same  way,  the  signal  if  can  be  viewed  as  the 
dominant  part  off*  causing  excitation  in  that  part  of 
the  spectrum  where  the  model  error  is  small,  e.g.  also 
at  low  frequencies,  loannou  and  Kokoiovic  (1983b) 
also  discuss  this  type  of  frequency  separation  in  the 
regressor  in  the  presence  of  high-frequency  para- 
si  tics. 

These  results  still  remain  incomplete  because  we 
need  to  know  the  relationships  between  m  and  the 
‘size'  of  f ,  e.g  Theorem  2  requires  a  bound  on  Hnj.il, , 
which  is  a  function  of  /,  in  and  consequently  f.  Of 
further  interest  is  the  effect  of  dead-zones  and  signal 


iimm. ili/.ilioii'  on  the  i. iiialiou.il  loim  (53). 
(  Vi !. only  the  n.iluie  ol  the  iiicmoi  \lc>s  iionlincaiity 
/  ( v  |  change-',  .t'  well  a'  the  system  .V 

(.  (  <>\<  I  l 'Ml  INS 

111  this  papei  we  have  piesented  a  framework  foi 
an  input  output  tlieoiy  of  adaptive  control.  This 
viewpoint  pi  ovules  a  means  to  realistically  de¬ 
termine  the  robustness  properties  of  adaptive 
algorithms  Moreover  input  output  concepts  are 
closely  related  to  measurement  techniques,  and 
hence,  can  lead  to  the  determination  of  usable 
engineering  techniques.  In  control  design  and 
analysis  the  most  notable  example  is  the  use  of  Bode 
plots  for  scalar  systems  (Bode.  1945)  and  singular 
value  plots  for  multivariable  systems  (Doyle  and 
Stein.  I9S1 ).  At  the  present  lime,  no  similar 
‘engineering  theory'  exists  for  adaptive  control 
design.  E/i  route  to  establishing  such  a  theory  it  will 
be  necessary  to  resolve  some  of  the  open  issues 
raised  herein.  The  possible  benefit  to  adaptive 
control  engineering  design  is  substantial. 

Acknowledgements —  R.L.K.  was  partially  supported  by  Air 
Force  Office  of  Scientific  Research  (AFOSR)  Contract  F 49620- 
83-C-0107.  C  R.J.vvas  supported  by  Nalional  Science 
Foundation  Grant  ECS-SI 19312. 

REFERENCES 

Anderson.  B  D.  O  (1977).  Exponential  stability  of  linear 
equations  arising  in  adaptive  identification.  IEEE  Trans.  Am. 
Coni  rid.  AC-22.  S3  SS. 

Anderson.  B  D  O.  and  C.  R.  Johnson.  Jr  1 1982a t  Exponential 
convergence  of  adaptive  identification  and  control  algorithms. 
Aulomatica  IK.  I  13. 

Anderson.  B  D  O.  and  C.  R.  Johnson.  Jr  ( 1 982b I.  On  reduced 
order  adaptive  output  error  identification  and  adaptive  1IR 
filtering  IEEE  Trans.  Am  Control.  AC-27.  927-  933. 
Anderson.  B  D  O.  and  R.  M.  Johnstone  (19831.  Adaptive 
systems  and  time  varying  plants.  Ini.  J.  Control.  37,  367  377. 
Astrbm.  K  tl983|  Theory  and  applications  of  adaptive 
control-  a  survey.  Automalica.  19.  471-486. 

Bode.  H,  W  (1943 1  Set  work  Analysis  and  Eecdhat  k  Amplifier 
Design.  Van  Nostrand.  Princeton.  NJ. 

Desocr.  C  A  and  M.  Vidyasagar  (1975).  Feedback  Systems: 

Inpul  Output  Properties.  Academic  Press.  New  York. 

Doyle.  J.  C.  and  G  Stein  (1981 ).  Multivariable  feedback  design: 
concepts  for  a  modern  classical  synthesis.  IEEE  Trans.  Aul. 
Control.  AC-26.  4  17. 

Egardl.  B.  |I979|  Stability  ol  Adaptive  Controllers.  Springcr- 
Vcrlag.  Berlin 

Gawthrop,  P.  J  and  K.  W.  Lim  (1982).  Robustness  of  self-tuning 
controllers  IEEE  Pro,  ..  129.  21  29 
Gornart.  O.  and  P.  Caines  ( 1984 1  Robust  adaptive  control  of  tune- 
varying  systems  McGill  University  Research  Report.  Depart¬ 
ment  or  Electrical  Engineering 

Goodwin.  G  C  .  IV  J.  Ramadge  and  P.  F.  Caines  (1980).  Discrete 
time  multivariable  adaptive  control.  IEEE  Trans.  Aut.  Control. 
AC-25.449  456 

loannou.  P.  A  and  P.  V.  Kokoiovic  (1983a).  Improvement  of 
robustness  of  adaptive  schemes.  Proceedings  of  the  3rd  Yale 
Workshop  on  Appluations  ol  Adaptive  Systems  Theory,  June 
1983.  pp  148  153 

loannou.  I’.  A  and  P.  V.  Kokoiovic  (1983b)  Adaptive  Systems 
null  Kediuvd  Models.  Spnngcr-Vcrlag,  Berlin 
Kosut.  R.  I.  1 19X3 1  Robust  adaptive  control:  conditions  for  local 
stability  Prot eedingsol the  1983  Automatic  (onlrot Conference, 
San  I  rancisco.  (  A.  June  1983 


36 


Input  output  view  of  adaptive  conliol 


kostii.  K  I  .mil  It  I  i  ifill. mill!  1 1 'IS?  I  Pci  loim. nice  robustness 
pi o|Vi lies  ol  .iil.ipliu'  coni l ol  sy stems  P im  codings  o/ I  In  ?  I  si 
III I  (  on/i  icin  i  on  /Viisioii  miJ  Lonlrol.  Oil.nuto.  II. 
I  Vccmhei  I 'IS?,  pp  IS  ?l 

K  os  in.  K  1  mill  H  I  ncillmulci  ( IIS')  Kohusl  adaptive  conliol ; 
conditions  foi  global  stability.  II. I  I  I  him  Aul  lonlrol 

Kosiil.  K  1CK  Johnson.  Jr  unit  li  l>  O  Anilerson  (1983). 
Conditions  foi  local  siahilit)  amt  robustness  of  adaptive 
control  systems.  Prmccdings  ol  the  ??n./  II  I. I  (  on/erem «•  on 
/lei  is  mu  miJ  (  null  ol.  San  Antonio.  1  X.  December  1983,  pp. 
9 7?  *I7<» 

kreisselmcici.  Li  and  k  S  Narendra  (I'ZS?)  Stable  model 
relercnce  adaptive  control  in  the  presence  of  bounded 
disturbances.  ILLL  '/runs.  Auto.  Control.  AC-27,  lib*/  1175 

Landau.  V.  1).  (1978)  Elimination  of  the  real  positivit>  condition 
in  the  design  of  parallel  MR  AS.  I  ELL  Train.  Am.  Cum  ml.  AC- 
23.  1015  1020. 

Monopoli.  R  V.  ( 19741  Model  reference  adaptive  control  with  an 
augmented  error  signal.  ILLL  Tran*.  Aul.  Control.  AC-19, 
474  4S4. 

Narendra.  k.  S„  Y.  H  Lin  and  L  Yalavam  (1980).  Stable 
adaptive  control  design  ILEL  Irani.  Am  Control.  AC-25. 
440  44S. 

Ortega.  R  and  I  D.  Landau  (1983).  On  the  mode!  mismatch 
tolerance  of  various  parameter  adaptation  algorithms  in  direct 
control  schemes:  A  secloricity  approach.  Proceedings  ol  the 


II  AC  II  ork  s/iop  on  Adapt  it,  (  out  ml.  San  I  lancisco.  <  A.  June 
1981 

Praly .  I  ( 1983 1  Robustness  ol  indirect  adaptive  control  based  on 
pole  placement  design  Prm rei/oigs  ol  tin-  II  Al  Workshop  on 
Adaptin'  (  out  ml.  San  Lraricisco,  CA.  June  1983 
Rohrs.  C.  I-..  L  Valavani.  M  Athans  and  Li  Stein  (1981 1 
Analytical  verification  of  undesirable  properties  of  direct 
model  reference  adaptive  control  algorithms  Proceedings  ol 
r/ic  ?0f/i  ILEL  Conference  on  Decision  and  Control,  San  Diego. 
CA.  December  1981,  pp  1272  1284 
Rohrs.  C  E..  L  Valavani.  M.  Athans  and  Li  Stein  (I982| 
Robustness  of  adaptive  control  algorithms  in  the  presence  of 
unmodeled  d>  nanucs  Proceedings  ol  r/ic  2 1  sr  IELI.  l  on/erem  o 
on  Decision  und  Control.  L>rlando,  EL.  December  1982. 
Safonov.  M  G.  (1980).  Stability  and  Robustness  ol  Miiltirmiuble 
Systems  Mi  l  Press.  Lamhridge.  Mass.  1980. 

Samson.  C.  (1983).  Slabilit)  analysis  of  adaptively  contiolled 
systems  subject  to  bounded  disturbances  Aiitomatica.  19. 
81  86 

Zames.  G  (1966)  On  the  input  output  stability  of  lime-varying 
nonlinear  feedback  systems  ILLL  Irons.  Am.  Control.  Part  I  : 
AC-1 1.  228  238;  Part  II:  AC-II.  465  476 
Zames.  G  and  B.  A.  Francis  (1983).  feedback,  minimax 
sensitivity,  and  optimal  robustness.  IEEE  Trans  Aul.  Control. 
AC-28.  586-601. 


A  LOCAL  STABILITY  ANALYSIS  FOR  A  CLASS  OF  ADAPTIVE  SYSTEMS 


by 


Robert  L.  Kosut* 
Integrated  System,  Inc. 
101  University  Avenne 
Palo  Alto.  CA  94301 


and  Brian  D.O.  Anderson 

Dept,  of  Systems  Engineering 
Australian  National  University 
Canberra,  Australia  ACTA2601 


ABSTRACT 


An  analysis  of  adaptive  systems  is  presented  where  a  local  L^- 
stability  is  insured  under  a  persistent  excitation  condition. 


* Research  supported  by  the  Air  Force  Office  of  Scientific  Research  (AFOSR) 
under  Contracts  FH9620-83-C-0107  and  F49620-8H-C-005‘». 


In  this  note  we  combine  some  earlier  results  [1)  -  (4]  to  provide  a 
framework  for  stability  analysis  of  adaptive  systems.  Ve  consider  here  the 
continuous-time  adaptive  control  of  a  scalar  plant*  with  input  s  and  out¬ 
put  y,  described  by 


Plant: 

y  =  d  +  Pu 

(la) 

Control : 

A 

u  =  —  0  *  X 

(lb) 

Adaptation: 

• 

A  ^ 

8  *  T*«,  0(0)  e  RP 

(lc) 

where  P  is  linear  with  strictly  proper  transfer  function  P(s),  d  is  an 

A 

external  disturbance,  6  is  the  adjustable  parameter  vector,  y>0  is  the  cons¬ 
tant  adaptive  gain,  z  is  the  regressor  (information)  vector  consisting  of 
filtered  measurable  signals,  e.g.,  u,  y,  and  references,  and  e  is  an  error 
signal  which  drives  the  adaptation.  System  (1)  can  also  be  described  in  an 
error  system  form  (e.g.,  [7],  [8])  by  proceeding  as  follows. 

Define  the  parameter  error  by 


e  :=  e-  e# 


(2) 


where  0  ,  i  lf  is  a  constant  vector  of  tuned  parameters  i.e.,  the 
parameters  that  would  be  selected  if  the  plant  P  were  known.  Using  (2)  ve 
can  rewrite  (lb)  as 


u  -  -e;  z  -  v 


v  :«  6'z 


(3) 


*  Extension  to  KINO  plant  is  straightforward,  e.g.  (31. 
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•here  v  it  the  adaptive  control  error  signal.  An  equivalent  representa¬ 
tion  of  (1)  it  given  by  the  adaptive  error  system  depicted  in  Figure  1  and 
described  by: 


e 


H  v 
ev 


(4a) 


r  *  x*  -  H  v  (4b) 

•  iv 

v  «  **0  (4c) 


0  *=  yze,  0(0)  *  0(0)  -  6, 


(4d) 


where  (et,  r#)  are  the  outputs  of  the  tuned  system  which  is  defined  as 
system  (1)  with  control  u  *  “6**f  The  operators  and  Hjv  are  linear 
with  strictly  proper  transfer  functions  H  (s)  and  E  (s),  respectively, 

XV 

which  are  dependent  on  the  tuned  parameter  0#.  From  the  definition  of  the 
tuned  system  [3J#  [4],  it  follows  that  H  (s)  and  H  (s)  are  exponentially 

CV  XV 

stable.  By  the  same  reasoning  the  toned  signals  e,(t)  and  x#(t)  are 
bounded. 


Fignre  1.  Adaptive  Error  System 


I 


One  of  the  very  useful  features  of  this  error  system  is  that  the  non¬ 
linear  effect  of  the  adaptive  algorithm  can  be  analysed  separately  from  the 
analysia  of  the  tuned  system.  The  tnned  system  represents  an  ideal  which 
would  be  achieved  with  the  given  structure  of  the  adaptive  control.  Bence, 
the  algebraic  design  procedure  is  separated  from  the  nonlinear  stability 
analysis.  It  is  convenient,  therefore,  to  view  e(>  and  6^  as  'inputs' 

to  the  error  system.  The  assumption,  naturally,  is  that  e,  and  t(  are  well 
behaved  with  e#  small.  Note  that  6^  need  not  be  small.  In  the  ideal  case, 
assuming  perfect  model  following  and  no  disturbances  in  the  tuned  system, 
e#(t)  =  0.  If  the  disturbances  are  of  a  special  kind  then  e$(t)  ->  0, 
i.e.,  the  tuned  system  exhibits  servo  action.  The  more  realistic  case,  how¬ 
ever,  is  when  e#  e  due  to  bounded  disturbances  which  cannot  be  asympto¬ 
tically  rejected. 


Global  Stability  Conditions 

By  global  stability  of  (2)  we  mean  that  all  bounded  inputs  e#,  x#, 
and  8_  produce  bounded  outputs  e,  6,  and  x.  In  general,  no  restrictions 
are  placed  on  the  initial  parameter  error  other  than  boundedness. 
Sufficient  conditions  for  global  stability  can  be  obtained  for  (2)  using 
passivity  theory  (e.g.,  [5],  p.  182).  A  detailed  analysis  can  be  found  in 
[3]-[4l.  One  of  the  conditions  is  that  B^ti)  is  strictly  positive  real 
(SPR),  i.e.,  B^t*)  is  strictly  proper*,  exponentially  stable  and  there 
exist  a  positive  constant  p  such  that 


Re  Hev(ja.)  2.  p|Bev<j«>|2.  V  a>  s  R 


(5) 


•When  B  (s)  is  proper  but  not  strictly  proper,  then  SPR  is  de 
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Re  »  <J»>  2i>0,  Vail. 


f 


Unfortunately,  H  (s)  t  SPK  la  not  robust  with  respect  to  even  slid  modeling 

error,  particularly  high  frequency  unnodeled  dynasties  [6).  For  example, 

H  (s)  e  SPR  implies  that  the  relative  degree  of  B  (s)  cannot  exceed  one, 
ev  ev 

from  which  it  follows  that  applying  this  restriction  to  (1)  imposes  the  same 
relative  degree  restriction  on  P(s)  as  well.  This  is  unrealistic,  even  in 
this  simple  example. 

Local  Stability  Conditions 

Conditions  for  local  stability  require  not  only  that  the  inputs  e#, 
tt,  and  6q  are  bounded,  but  that  these  hounds  are  not  arbitrary.  The  local 
analysis  is  facilitated  by  transforming  the  error  system  (4)  to  the 
variational  form 


i  B  -  Gf(x)  (6a) 

where  x,  x^,  G,  and  fix)  are  defined  by 


(6  f ) 


M  Ni; 

K  :=  (1  +  LM)_1L  (6g) 


and  where  L  has  the  transfer  function, 

L(s)  =  *  r  (6h) 

with  y  from  the  adaptive  algorithm  (1c).  This  error  system  (6)  is  arrived 
at  by  separating  the  nonlinear  cross  product  terms  in  f(x)  from  the  linear 
We  shall  refer  to  as  the  response  of  the  1 inearized 
system.  This  is  almost  identical  to  the  linearised  system  studied  by  Bohrs, 
et  al.  [6a],  which  was  arrived  at  by  a  ’final  approach  analysis.'  Note  that 
in  this  case  the  linearized  system  is  the  input  to  the  nonlinear  system. 

The  operators  I  and  G  are  linear  and  time-varying  due  to  their  dependence 
on  the  tuned  signals.  If  the  linearized  response  z^  in  (6c)  is  small,  and 
if  the  nonlinear  term  f(z)  is  suitably  restricted,  then  intuitively,  x 
would  be  attracted  to  some  neighborhood  of  x^.  The  following  theorem  makes 
this  notion  precise.  We  use  the  notation  y  (.)  and  II  II  to  denote  L  - 
gain  and  L— -norm,  respectively. 

Theorem  1 :  Suppose  there  exist  finite  positive  constants  g,  e,  and  6(e) 
such  that 


yJC)  ±  g  <  1/a 


(7a) 


|x|  <  6(c)  •>  lf(x)  |  <  till 


(7b) 


I  l*Ll I.  1  (1-g  t)  6(c) 


(70 


impl iei 

I  lx  I  I  16(c).  (7d) 

m> 

Theorem  1  follows  directly  from  the  1 inear ization  theorem  of  [5,  p.  131). 
Theorem  1  asserts  that  the  error  outputs  z  of  the  adaptive  error  system 
are  L^-bounded  in  an  e-neighborhood  of  the  linearized  response,  provided 
that  the  linearized  response  is  small  enough  and  that  G  e  L^- stable. 
Condition  (7d)  shows  that  the  actual  response  can  be  arbitrarily  close  to 
the  linearized  response.  Since  Theorem  1  provides  sufficient  conditions, 
instability  does  not  follow  if  z^  e  L*  but  ezceeds  the  magnitude 
constraint  of  (7c). 

The  function  6(c)  in  (7b)  can  be  determined  from  the  definition  of 
f(z)  in  (6b)  and  the  norm  selected.  For  ezample,  if  the  norm  on  Rn  is  de¬ 
fined  as  I z  I  =  Maz  |z.|  and  l|z||  *=  snp|z(t)l,  then 

1  • 
i  t 


6(c)  «  e 

and  using  the  corresponding  induced  matriz  norm,  we  obtain 


(8a) 


g  -  maz  (f1#  $2). 

*1  "  *0(1  +  ,,*JL,k(1  +  n))  (8b) 

S2  *  k(l  +  n) 


where 


t0  2  T.<»gv>J 

*  2  tJN),  k  2  T.<I>  <«•> 
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Although  Theorem  1  provides  conditions  for  local  L^-stabil  ity,  these 

do  not  immediately  provide  a  region  of  attraction,  i.e.,  bounds  on  e#,  x#, 

and  o  •  These  bounds  in  tnrn  are  determined  from  the  aet  of  allowable 
o 

reference  commands,  plant  initial  conditions,  and  disturbances.  Since  e, 
and  x,  are  bounded  by  predetermined  performance  goals  of  the  tuned  system, 
it  follows  that  is  the  unknown  driving  factor  governing  the  sixe  of 
I  lx,  II  .  That  the  initial  parameter  error  vector  occupies  this  position  of 
villainy  should  come  as  no  surprise.  One  way  to  offset  large  initial  para¬ 
meter  errors  is  to  keep  the  adaptation  gain  y  small.  This  has  the  effect  of 
reducing  large  system  transients,  however,  this  may  be  less  than  prudent  if 
the  system  is  initially  unstable  or  lightly  damped. 

No  claims  are  made  in  Theorem  1  about  the  mechanism  that  provides 
x^e  L”  and  G  e  L^-stable.  However,  it  follows  from  the  definition  of  the 
tuned  system  that  e„  e  L  ,  r .  e  Ln  and  H  ,  H  e  L  -stable,  thus,  M  in 
(6f)  is  L  -stable.  Hence,  a  term  by  term  inspection  of  G  ( 6 d )  and  x. 

®  a,  L 

(6c)  reveals  that  x,  e  L**  and  G  e  L  -  stable,  if  and  only  if  0.  t  L^. 

L  •  ff  L 

Looking  at  (6c)  we  can  also  describe  e^(t)  •*  the  solution  to  the 
differential  equation, 

t(t)  *=  -y(M?)(t)  +  Y  w(t)  (9) 


with  w  «  x#e#  and  ((0)  *  6Q.  Referring  to  (6)  and  (9),  the  operator  K  is 
equivalent  to  the  mapping  from  w  into  (.  Hence,  the  stability  analysis 
of  (9)  is  of  fundamental  importance. 

Persistent  Excitation  and  Exponential  Stability 

Equations  similar  to  (9)  have  been  studied  by  invoking  a  persistent 
excitation  condition  on  x0(t).  The  following  definition  and  lemma  from  [1] 
provides  the  basic  result. 

Definition;  A  regulated  function  f(.):R+  ->  R*  is  persistently  exciting, 
denoted  f  a  PE,  if  there  exist  positive  constants  e. ,  o_,  and  a,  such  that 


(10) 


•♦a. 


s 


f(t)f(t)*dt  i  a,  1  . 

i  D 


V*  t  R_, 


Lemma  1 :  Consider  the  differential  equation: 


((t)  =  -rf(t)  (Bf'tHt)  +  y*(  t).  t2  0 


(11a) 


If  f  e  PE  and  H(s)  t  SPR  then  the  map  (((0),*)  l->  (  is  exponentially 
stable,  i.e.,  there  exist  positive  constants  m  and  X  soch  that. 


U(t)l  <me“Xt|t(0)l  +/  me-X(t_T) |w(T) IdT 


(lib) 


The  nsefnlness  of  applying  Lemma  1  to  determine  stability  conditions 
of  (9)  is  made  apparent  by  writing  B ^  as. 


B  =  B  +  B  (12) 

ev  ev  ev 


where  B  is  the  nominal  representation  of  B  and  B  is  the  deviation 
ev  ev  ev 

indneed,  for  example,  by  modeling  error.  Combining  (12)  with  (9),  and  using 
the  definitions  in  (6)  gives. 


%  *  -r*.5ev*;?  ♦  t<U  ♦  r* 


(13a) 


where 


Q  M  -  xj) 


a"  *i 
•  ev  • 


-  *  B  %• 
•  ev  • 


•  IT  < 


(13b) 
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If  llfv  (s)  e  S1K  and  e  PK,  then  using  Lemma  1  gives. 


Uc t >  I  i 


■e"XtU(0)  I 


t 

/ 


yme 


-X(t-x) 


l(QO(T)  *  w(t)ldt  (Id) 


Hence,  k  from  (8)  is. 


k-f  i  *-(I) 


( ISa) 


and  from  (14)  vith  ( 


replaced  by  6  ve  get. 


I  l©Ll I.  i  (l  -  r»q)  1  [ I (^1  ♦  r*l l./x] 


(15b) 


provided  ymq  <  1  where 


q=  H*Jlf  T-<®ev)  +  ^  ^«(Q)  *  <15®> 

Combining  (8),  (15),  and  Theorem  1  gives  the  following  result. 


Lemma  2:  The  adaptive  system  (1)  or  (2)  is  locally  L^-stable  if  for  some 
a  <  1/g. 


o*l 


Iepl  +  T»l !»«•,! 1,/X 

(1  -  g  «)  • 


>  0 


and 


ml* 


i 


(Ida) 


I 


(1(b) 
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I.ctnmt  2  together  with  (15)  and  (8)  provides  sn  explicit  upper  bound  on 

1 1().  1 1  ,  It'  I,  and  the  amount  by  which  H  can  deviate  from  a  nominal  H 
1.  •>  G  ev  ev 

which  is  SPK.  If  the  bounds  are  satisfied  then  Theorem  1  asserts  that  the 
signals  in  the  adaptive  system  (1)  are  all  bounded. 

Unlike  the  global  stability  case  where  the  bound  on  the  deviation  B 

cv 

is  severely  restricted,  the  bound  here  can  be  large. 

Concluding  Remarks 

The  stability  analysis  provided  here  involves  establishing  the 
exponential  stability  of  a  differential  equation  (9)  which  arises  in  the 
study  of  most  adaptive  systems.  Although  the  connection  between  exponential 
stability  of  (a)  and  persistent  excitation  is  known  [1],  it  is  important 
here  to  obtain  specific  formulae  for  the  rates  and  gains  involved,  e.g.  (8), 
(15),  (16).  Other  methods  to  obtain  these  values  can  be  found  in  [9]  and 
[10].  Note  also  that  Theorem  1  only  requires  L^-stabil ity  which  is 
certainly  provided  when  (9)  is  exponentially  stability.  However,  In¬ 
stability  can  be  obtained  by  using  a  nonlinear  adaptation  gain  in  (lc), 

A 

i.e.,  0  •=  yh(x.e).  For  example,  h(z,e)  can  arise  from  using  a  dead-zone, 

leakage,  or  normalization  [11].  Such  schemes  can  be  incorporated  in  the 
general  framework  presented  here  but  require  further  analysis  in  order  to 
obtain  explicit  signal  bounds. 
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ABSTRACT 


2.  NOTATION 


The  question  is  examined  of  uhes  an  adaptive 
control  systees  it  robust  to  unaiodrled  dynamics  and 
unknonn  bounded  disturbances.  Conditions  arc 
presented  tbat  ensure  tbe  eaistcnce  of  such 
robustness  properties,  but  only  locally;  i.e., 
restrictions  are  placed  on  tbe  behavior  of  signals  in 
tbe  ideal,  perfectly  tuned  adaptive  ayttca.  Local- 
L^-ttability  it  investigated  when  certain  tuned 
aianals  are  attuned  to  be  persistently  exciting. 


Let  L  denote  tbe  act  of  Lebetgue  integrable 
functions  E(.):  ->  S*  with  finite  nora  llsll  : 

-  <  /  llx(t)llpdt)1/p  for  p  Cl,-)  and  llxll_ 
anp  Tx(t)l.  where  I. I  denotes  a  none  on  l*. 

Hmilsrly,  l*t  i.n  •  denote  the  ei|cajjOB  of  L*. 
consisting  of  fnnPfiona  x(.)  aneb  that  x_  e  L  .  VI 
2  0,  where  x_(t)  denotes  tbe  truncation  of  x?t) 
at  T.  i.e..  x  <t)  -  x(t>  for  t  i  T,  and  x_(t) 
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Theoretical  inveatiga tiona  on  tbe  stability  of 
adaptive  control  ayateaa  have  focnaed  alnost  entirely 
on  developing  conditions  tbat  guarantee  alobal 
stability,  e.g..  Ill— 131 .  These  rcanlta  are  global 
in  tbe  rente  that  initial  conditions  and  external 
signal  magnitudes  need  only  be  bounded.  Specific 
bounds  are  not  reqoired.  In  addition,  the  results 
provide  sufficient  conditions.  One  of  tbe  conditions 
ia  tbat  a  particular  subsystem  operator  be  strictly 
passive  with  finite  gain  or,  in  tbe  case  of  linear- 
time- invariant  systems,  tbe  operator  ia  strictly 
positive  real  (SPR).  This  condition  results  from 
application  of  tbe  Passivity  Theorem;  specifically, 
tbe  adaptive  system  can  be  reconfigured  into  two 
subsystems:  a  'feedback*  subsystem  (the  adaptation 
law)  tbat  is  passive,  and  a  'feedforward'  subsystem 
which  is  required  to  be  SPR.  This  condition  turns 
out  to  be  qoite  restrictive.  In  the  first  place,  tbe 
SPR  condition  nocessitetes  tbat  tbe  system  transfer 
function  (in  tbe  scaler  cate)  have  a  unitary  relative 
degree.  At  pointed  out' by  rfohrt,  at  al.  (4),  it  in 
virtually  impossible  to  guarantee  unitary  relative 
degree  for  an  actual  system.  Secondly,  tbe  SR 
condition'  bat  extremely  limited  robustness  to 
nnmodcled  dynamics  (3] . 

In  this  paper,  conditions  ere  developed  that 
gnarantee  tbe  existence  of  )ocal  stability  and 
robustness  properties  of  tbe  adaptive  system,  i.e., 
conditions  which  take  into  acount  the  slxe  of  initial 
parameter  error  and  external  signal  magaitodea. 

These  conditions  are  imposed  on  certain  subsystem 
operators,  which  have  a  time-varying  dapeadeace  on 
signals  that  arise  from  an  ideal  fictitious  system 
where  tie  adaptive  gains  are  perfectly  tuned  to  the 
unknovi  t  to  be  controlled.  The  mechanism  for 

local  ati.  tty  which  ia  examined  here  ia  that  of 
persistent  excitation  (7),  (Jl.  Under  these 
conditions,  we  develop  a  specific  hound  on  model 
error  which  ensures  conditions  for  local  stability. 


la  this  section  we  present  an  adaptive  error 
ayatem  which  is  representative  of  a  large  class  of 
adaptive  control  systems.  Tbe  error  system  will  be 
presented  in  two  forms:  a  parameter  variational  form 
and  a  full  variational  form.  Tbe  parameter 
variational  form  was  developed  ia  detail  in  (Sb)  and 
is  used  for  global  stability  analysis.  Tbe  fall 
variational  fora,  to  be  developed  here,  la  nsed  for 
local  stability  analysis. 


To  facilitate  tbe  develojaent  of  tbe  error 
*T*ten.  consider  tbe  simple  model  reference  adaptive 
controller  (KRAC)  depicted  in  Fig.  3-1  with: 


y  •  d  ♦  Pu 


<3. la) 


d  ;•  external  disturbance  ♦  plant  initial 
conditions 


y.  ‘  Br 

*  reference  coauaaad 


(3.1b) 


.  - 1>, 


(3.1c> 


•  :•  adaptive  gains,  a  regressor 
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(3.2) 


•  here  4*(  •  ((*  ,  #,,)  1  it  a  cunataat  vector  o( 

tuned  game;  i.e.,  the  valuea  that  could  ba  aelectad 
il  tbe  plant  P  cere  knovu.  lit  log  (3.2)  ce  can 
recritc  (3.1c)  at 
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(3.3) 


cherc  v  ia  tbe  adapt ive  control  error  aignal.  An 
equivalent  repreaentation  of  (3.1)  ia  given  by  tba 
adaptive  error  ayatea  depicted  in  Figorc  3-2  and 
*!t  imbi  J  b>  : 


e  -  e  -  U 


i  -  B  v 

•  IV 


z‘6 
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%  +  Ue 


(3.4a) 

(3.4b) 

(3.4c) 

<3.4d) 


chere  (e  ,  x()  are  tbe  ootputa  of  tbe  tnned  avataa. 
at  sbocn  in  Figure  3-3;  it  the  initial  value  of 

tbe  adaptive  gain  error,  and  H  8  ,  and  L  are 

tbe  interconnection  operatora.  For  the  aiaple  MRAC 
cate  conaidered  berc  (Fig.  3-1),  tbe  toned  aignala 
are : 


(1  ♦  P»#1)  *d 


c-1. 


♦  Kl  ♦  M#1)  "Pe.j  -  Br)r 


«,  *  («.  ♦  yt.  -*)* 


(3.3a) 

(3.3b) 


and  tbe  intcrconacctiona  have  tranafer  fnnetiona. 


"1. 


R  (a)  -  (1  ♦  P(a)».,)  *P(e> 
ev  *1 


Hiv(t)  -  1(1  4  P(a)«.1)"1P(a).  0)* 


L(a)  •  (l/a)B 


(3.6a) 

(3.6b) 

(3.6c) 


Altboogb  tbe  error  ayatea  (3.4)  baa  been 
developed  here  for  a  very  aiaple  MRAC  ayatea,  tbe 
fora  of  (3.4)  ia  generic  and  appliea  to  practically 
all  aingle-iapot-aiagle-ootpot  adaptive  controllara 
and  filtcra  (3).  Moreover,  tbe  exteaeioa  of  (3.4)  to 
tbe  aoltivariable  caae  reqoiree  only  that  v  end  e 
are  vectora  and  that  B  and  B  are 
aultivariable  of  coapetIKle  diaenalona. 

Specifically.  (3.4c)  and  (3.4d>  are  replaced  by 


v  -  Z’6 


•  -  6#  ♦  LZe 


(3.4c)' 
(3.4d)  * 


ch era  Z  ia  a  block  diagonal  aatrla  of  appropriate 
diaeaaioaa  aneb  that 


Z  -  diag  (tj. 


. .  a_) 


(3.4e)‘ 


ia 


a'  -  ( ij . a;> 
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ilia 


ill  be 
of  (3.4), 
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r  air  ion  to  tl, »  au  1 1 1  v  1 1  1 1  b  1  r  cate  folloct 
inne d 1 1 1 e  I  y  . 

One  of  tbe  very  uteful  ieaturaa  of  tbit  error 
tytlea  it  that  tbe  nonlinear  effect  of  tbe  adaptive 
algoritba  can  be  antlyted  teparately  froa  tba 
analytic  of  the  tuned  eyttea.  Tba  tuned  ayatea 
repraaentt  an  ideal  ebieb  could  be  achieved  citb  the 
given  atrnctura  of  tba  adaptive  control.  Bence,  tba 
algebraic  detign  procedure  it  teparated  froa  tba 
nonlinear  liability  analytia.  It  it  convenient, 
therefore,  to  vlev  e,.  i,.  and  8.  aa  'inputt*  to 

the  error  tytlea.  Tbe  eaiuaption,  naturally,  ia  that 


and 


are  veil  behaved  vitb  a. 


aaall.  Note 


that  8„  need  not  be  aaall.  In  tbe  clattic  cate. 


attuning  perfect  node!  following  and  no  ditturbancea 
in  the  tuned  ayttea.  e,(t)  0.  If  tbe  diatorbancet 

are  of  a  tpecial  kind  then  e,(t)  ->  0,  i.e.,  tba 
tuned  ayatea  exhibit!  tervo  action.  The  aore 
realittic  cate,  bovcver,  ia  vben  e_  e  L  doe  to 
bounded  ditturbancea  which  cannot  be  aayapt oti cal ly 
rejected. 


3 .2  Global  Stability  Conditiont 


Condition!  for  global  atability  reqnira  that 
B  (a)  a  SPk.  Tbia  ariaaa  became  prooft  of  global 
afability  utilize  paaaivity  theory  (c.g.,  (91,  p. 
182).  A  detailed  analytia  can  be  fonnd  in  tl)~(3]. 
Unfortunately,  though  of  theoretical  aignificance, 
theae  type  of  reaultt  do  not  offer  any  practical 
enginaaring  guideline!.  Tbe  aajor  reaton  ia  that 
B  (a)  e  SPK  ia  not  robuat  vitb  retpect  to  even  aild 
modeling  error,  particularly  high  frequency  uaaodelcd 
dynaaict  [4].  Since  B  (a)  a  SPK  inpliea  that  tbe 
relative  degree  of  B  fH)  cannot  exceed  one,  it 
followt  that  apply iag*Ihia  reatrictioa  to  (3.6a) 
iapoaea  a  unitary  relative  degree  reatrictioa  on 
P( a)  a*  well.  Tbia  ia  uarealiatic,  even  la  tbia 
aiaple  exaaple. 


Another  view  of  tbe  reatrictiveneaa  of  B^^  a 
SPK  it  froa  robuatneaa  theory,  e.g.,  [IS).  Suppoac 


that  P( a )  ia  (3.1a)  can  be  exprcaacd  aa  beloagiag 
to  tbe  aet  of  tranafer  fuactioaa 


P(e)  «  (1  4  A( a ) )P(a ) 
|A(  jw)  I  S.  6(a).  (bit 


(3.7a) 

(3.7b) 


Bence,  P(a)  ia  a  aoaiaal  aodel  of  P(a)  and  A(a) 
repreaenta  nodding  error,  e.g.,  bigb  frequency 
uaaodelcd  dynaaica.  We  can  aov  write 


B  -  E  4  1 

av  ev  ev 


(3.8a) 


where  tba  aoaiaal  ia. 


-1 


B  -  (1  4  P»  )  P  e  SPK  (3.8b) 

ev  •! 

and  tba  deviation  induced  by  aodel iag  error  ia. 
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ev 


i  li  ♦  (i  ♦  »..)"*  (j  ♦  pe.. )_1a 

ev 


•3 


•3 


(3.8c) 


It  ia  abowa  ia  (3)  that  the  largeat  tolerable  8(a) 
ia  (3.7b)  to  eaaare  a  SPK  ia  bounded  by 

6(a)  <  1  .  (3.8) 


Again,  tbia  ia  marealiatie  and  ia  violated  even  by 
tbe  aoat  aild  fora  of  aaaodeled  dynaaica.  Note  that 
(3.8)  and  tba  aaitary  relative  degree  reatrictioa 
both  aaceaaarily  ariae  froa  tbe  tn  coaditioa. 
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than 


11>«  nror  i>tl*>  0.4)  can  be  tranafixeied  «o  tha 
td  loving  variational  lor*  abltb  ia  aort  uaeful  for 
local  liability  anelyaia,  I. a., 

a  *  *L  -  Cf (a)  O .10a) 

•bare  (ha  quantiliea  above  are  defined  below  by 


"•II.  i  ■. 

(ii)  In  addition,  if.  foe  tone  p  •  ||,*),j 

conatant  a  aucb  that, 

P 

T  <C)  i  g  <  2/a  (4.4) 
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(3.10c) 


(3 . lOd) 


Ni* 


(3 .10f) 


i  <1  -  I ll»,  II  (4.J) 

P  P  *  L  P 

lto.pl : 

Theoren  4.1  ia  baaed  on  tha  1 inaariaation 
tbeoren  of  (9,  p.  131).  The  proof,  aa  apacialixed 
here,  ia  in  Appendix  A. 


Benefit 

(1)  Tbeoren  4.1.  part  (i).  aaaerta  that  tha 
error  ontputa  a  of  the  adaptive  error  ayatan  are 
L  -bounded  in  an  e-neighborhood  of  the  linearized 
reaponae.  provided  that  the  linearized  reaponae  ia  in 
L_  and  ia  mail  enough  (4.2).  and  that  C  a  Le¬ 
ttable  (4.1).  Condition  (4.3)  ahowa  that  the  actual 
reaponae  can  be  arbitrarily  cloae  to  the  linearized 
reaponae,  if  llx  II  >■  anall  enough  (4.2).  Since 
Theoren  4.1.  part  ( iT,  provider  aofficient 
condition!,  inatability  doea  not  follow  if  x  t  L_ 
but  azcecda  the  nagnitnde  conatraint  (4.2). 


I  :«  (I  ♦  LM)_1L  (3 . 1 0 g ) 

and  where  L  haa  the  tranafer  function. 

Lit)  *  J-  B  (3 .10h) 

with  B  fron  the  adaptive  algorithm  (3. Id). 

The  nodcl  (3.7)  ia  arrived  at  by  aaparating  tha 
nonlinear  croaa  product  terne  in  f(z)  fron  the 
linear  tenax  ia  x.  .  Be  ahall  refer  to  x,  aa  the 
reaponae  of  the  linearized  avaten.  Thie  ia  alnoat 
identical  to  the  linearized  ayaten  atudied  by  Kohre, 
et  al.  [4a],  which  wax  arrived  at  by  a  'final 
approach  aaalyxia.'  Note  that  in  thia  caae  tha 
linearized  ayaten  ia  the  input  to  the  nonlinear 
ayaten  (3.10a).  The  operator!  K  and  G  arc  linear 
and  tine-varying  doe  to  their  dependence  on  tha  twnad 
aignala  (a,,  a,).  Thia  Mdcl  (3.10)  will  now  be 
ntilited  to  develop  local  atability  condition!. 


4.  CONDITIONS  FOB  Mm  STABtUH 


If  the  linearized  reaponae  x.  ia  (3.10)  ia 
anall,  and  if  the  nonlinear  tarn  G  f(x)  ia  anitably 
reatricted,  then  intuitively,  x  would  ba  attracted 
to  tone  neighborhood  of  x,  .  Tha  following  theoren 
aakee  thia  notion  preciee.L 

PW«F  4.1 

(i)  If  3  conatant  g.  each  that. 

T.(C)  1  |.  <  -  (4.1) 

aad  if 

llz,  II.  i  a.  (1  -  I.I./1).  i.  a  (0.  2/g.) 

L  74.2) 


(2)  The  reanlta  in  part  (ii)  are  atronger  than 
in  part  (i)  aince  they  can  only  be  applied  when  x^  e 
L°  for  aone  p  a  [1,*>.  Looking  at  (3.10).  thia  can 
oBly  occur  if  a.e.  e  Ln  which,  ia  practical 
aituationa.  alnoat  aevcf  occur!  dne  to  the  preacnce 
of  diatnrbancec  la  L..  Bence,  part  (ii)  of  Theoren 
4.1  doea  not  offer  any  practical  advice  aad  wc  will 
focua  only  on  part  (i). 

(3)  Fron  (4.2),  the  largeat  upper  bound  on 
||xLl 1^  la  l/2g.  which  occnra  when  a.  “  1/g.. 

(4)  Although  Theoren  4.1  part  (i)  providec 

condition,  for  local  L.-atabil ity ,  theac  do  not 
iaacdiatcly  provide  a  rezion  of  attraction,  i.a.. 
bouada  on  a,,  a,,  aad  0  .  Theac  bonada  in  turn 

are  deteraiaed  fron  the  aet°of  allowable  reference 
connanda,  plant  Initial  condition.,  and  dlctorbaacea. 
Since  e.  and  i,  are  bounded  by  predetermined 
performance  goala  of  the  tuned  ayaten,  it  followc 
that  •  ia  tha  unknown  driving  factor  governing  the 
aize  of0  I  lx  I  I..  That  the  initial  paraaetcr  error 
vector  occnpaca  thia  poaition  of  villainy  abonld  cone 
aa  no  aurpriae.  For  example,  if  ia  anall  (order 

a  )  than  the  adaptive  ayaten  atayc  near  the  tuned 
ayaten  for  anall  (order  a.)  inputa  e.. 

(5)  No  claina  arc  nade  in  Theoren  4^.1  part  (i) 
about  the  nechanien  that  providea  x,  a  L.  end  G  a 
L  -atable.  Bowevar.  it  follow,  fron  the  definition 
of  the  tuned  ayaten  (3.2)  that  e.  e  L  .  a.  e  L. 
and  B  .  B  a  L.-atable.  thwa.  M  In  (3.10f)  ia 
L.-etaSfe.  ieaca,  a  tern  by  tern  iacpection.of  G 
(3.104)  and  x  (3.10c)  reveal,  that  x^  a  L.  and 

G  •  L.-atable,  if  and  only  if: 

(I  ♦  LB)”1*  •  L*  4W  a  B*  (4.4a) 

o  m  o 

•■4 
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K  a  L  -.table 


(4.4b) 


Ul»  I'ltuttl)  .  *«•  hiw  l>i*  lot  low  iii|  frtull. 


l4<ol  i  me  Utoil  ♦  /  ns  ***  'NtiilUt 


1»»|  4jl 

Surpou  that  I  tuned  solution  6(  t  R°  dials. 
Let  ((t)  denote  the  eolution  nt  t i me  t  o(  the 
dilfeteutiel  equation. 

(<t>  -  -H(MlMt)  «  eft)  t  2  0  (4.7) 

Then.  i  e  L*  and  (i  e  l.^-atabl*  if  and  only  If  ( 
e  L*  for  all  t(0)  ■  R°  and  ■  e  L*. 


Using  the  dcfinitiona  of  K  and  M  In  (3.10). 
if  4(0)  •  0  then  (4.7)  represents  I  :  a  |->  4. 
Thus.  K  e  L^- st able  if  «  |->  (  >  L^- at able.  Alao. 
if  4(0)  «  8*  and  «  *  *,*,  then  4  *  6^  Iron 
O.10c£.  Bence,  uitboot  any  further  restrlctiona  on 
0  a  R  or  te  a  L°,  the  reaolt  of  Least  4.1  la 
establ ished. 


Condition!  for  Robustness 

The  usefulness  of  applying  leant  3.1  to 
detsrnine  robust  stability  conditions  of  (4.7)  is 
nads  apparent  by  proceeding  at  in  Section  3.  1.#., 

U  -  H  sfi  (3.4) 

av  ev  av 

•here_  B  la  the  noainal  representation  of  B 
and  B  it  the  deviation  induced  by  nodeliog  error. 
Coabinfog  (3.4)  uith  (4.7),  and  osing  the  definitions 
in  (3.10)  gives. 


-Bi.fi  rU  ♦  Q4  •»  • 

•  ev  • 


Q  B(M  -  x  H  »!) 

•  av  • 


If  H  (a)  a  SPR  and  a,  a  SR.  than  (3.3)  of 
Lena  (3.11  applied  to  (3.3)  gives. 


Leans  4.1  identifies  the  L^-stability  of  the 
tystec  (4.7)  at  being  crucial  to  obtaining  local 
stability  conditions  froa  Theorem  4.1.  Thia 
condition  it  not  sufficient.  Even  if  the  conditions 
of  Lcauna  4.1  are  satisfied  the  adaptive  aystea  ia 
locally  stable  provided  that  I  lx  1 1^  ia  aaall 
enough,  i.e.,  (4.2)  Bust  hold.  Nonetheless, 
establishing  (4.9)  ia  a  first  step. 

Recall  that  1 1 a  1 ia  aaall  if  1 1 ®L 1 1 —  is 
aaall,  hence  it  is  necessary  to  control  the  sixe  of 
these  signals.  Coaparing  8.  in  (3.10c)  to  (4.7). 
soac  of  these  aagnitude  conditions  can  be  secured  by 
lowering  the  adaptation  gain,  i.e.,  the  norm  of  the 
astrix  B  in  (3. Id)  or  (3.10h). 


14(0 1  i  aiXtl4(0)l  ♦  /  aeXa“*,|(Q4>(t)  ♦  w(0|dt 


Therefore,  if  Q  has  a  sufficiently  aaall  gain  then 
(4(0).  v)  |->  4  a  L_- stable,  and  hence,  the  adaptive 
systea  is  locally  L^-atable.  Specific  conditions  ere 
given  ns  follows. 


Suppose  x.  s  PE,  B(s)  s  SPR,  and  hence,  froa 


Suppose  i,  s  PE,  H(a) 
5.1  4(0  in  {5.5*via 


(5.7).  Then,  the  adaptive  systea  (3.10)  is  locslly 
L  -stable  if,  for  soae 
o  s  (0,1), 


Tj(F>  <  1  -  q 


In  thia  section  we  exaaine  persistent  excitation 
as  s  aechaaisa  to  provide  L  -stability  of  (4.7),  and 

hence,  local  L  -stability  of  the  adaptive  aystea 
(3.10). 

Def  tuition  ( 8] 

A  regalsted  function  f(.):R+  ->  Rn  ia 
nersistentlv  excitins.  denoted  f  s  PE,  if 
positive  constants  o^,  Oj.  and  Oj  such  that 

s+o, 

a,  l§  j  /  *  f(t)  f(t)*dt  i  oa  IB.  Vs  a  R+ 

*  (5.1) 

The  relationship  between  persistent  excitation  and 
stability  of  (4.7)  ia  given  as  follows. 


4(0  -  -f(t)  <804X0  ♦  w(t).  »20  (5.3) 

If  f  s  PE  and  B(s)  a  SIR  then  the  aap  (4(0), w) 
I - >  4  1*  exponentially  atabla.  i.e..  there  exiat 

positive  consteats  a  and  A  each  that, 


«  "  ,lee|l.ll,*.M.T2(H“T)  (5.9) 

and  the  operator  F  has  the  integral  fora, 

(Fn) (t)  -  /  ,#<,)<i^*;,)(T)dT 

*  (5.10) 


Proof :  See  Appendix  B. 


(1)  The  a-snperscript  notation  ■*  aeans  that 
if  I  has  transfer  function  1(a),  then  I*  haa 
transfer  function  l(s-n).  Thus,  n  a  (0,1)  ia 
fnrthar  Halted  ao  that  I  (a-o)  and  I  (a-o) 
reaain  exponentially  stable;  otherwiaa  the  L  -gains 
in  (5.0) .  (5.10)  are  infinite.  1 

(3)  that  Theorem  5.1  asserts  is  that  if  a,  a 
PE,  and  If  B  Is  aloes  eaongh  to  Being  BOB,  than 
under  suitsbls'lasll  gain  conditions  (5.S),  loaal 
at  ability  can  he  gnaranteed  via  Theorem  4.1.  The 
arnr  of  the  natter  ia  te  satehlieh  that  y,(F)  n la 
anfficieatly  eaall  despite  a  reasonably  large  5^. 
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fjJi  u.ltjli.ji  i'J  .  JjIU 

Intuitively,  if  the  itnge  of  doalnmt 
fitgukiiik  >f  t(  it  aufficitnlly  teptrated  fro* 

tht  I rcqur tc it t  vhcre  tl  (jw)  It  large,  then  tj(F) 

•  ould  be  imII.  t.|..  i*v  It  pertitlanlly  exciting 

•  I  *  low  *  frtqurncict  where  H  it  approximately 
SPK.  be  ail  I  formal  tic  tbit  SXtion  in  Tbaorca  5  .2 
btloa.  lint,  bortu  i ,  at  need  tbt  following  reaultt 
fro*  | JO)  for  deteraming  a  large  clttt  of 
pertiatenlly  exciting  tignalt. 


(3)  Mine  a  t  (0,1)  In  (3.13)  ctn  bt 
arbitrarily  teiall,  tnd  tirrr 

q  <<  1  it  likely  due  to  ll»fllt  bring  tatll.  It 
folloat  that  t  rcaionablr  approx ime t ion  to  tbt 
robutlnett  teat  (3.13)  tt 


lb  <  jw>  I  < - j - : — 

tv  Bp  tup  la  (w  >1 

J  ,  %2.1/t  *• 


inf 

k 


l(«-«k)‘  ♦  1*1 
1*^  (w^) I 


(full 


(3.14) 


Definition  3.1  1 10] 


A  function  f(.):  B  ->  b"  hat  a  tgectral  lint 
at  frequency  u  of  atiplltude  tj(u)  *C  if 

lin  1  /  f(t)eJ  dt  -  *  (up)  ,  Vt  t  B 
.  t  f  ♦ 

x~>m  1  (5.11) 

ahen  a^(ui)  I  0,  f  bat  a  toectrtl  line  * t  a. 

Leant  3,2  [101 

Suppote  f  e  L  bat  tpectral  linet  at 

ti  of  aaplitudet  a.(w  ),  .... 


frcquencict  n  ,  ....  _ 
a  (w  ).  Then  t  PE.  9 it 
t  P 

rank  |a,(ui.) . e,(ui  )]  »  n 

11  1  P 


w 

(3.12) 


JbeoreB  }_.2 

Suppote  that: 


6.  CONQ  l'niNC  REKAKBS 


6.1  Tett  Procedure* 

The  retulti  of  Tbeoren  S.2  can  be  of  practical 
ute  tince  they  provide  the  batit  for  developing 
robutlnett  teat  proceduret.  The  aott  obviont  nay  to 
apply  Theorca  3.2  it  to  determine  the  model  error 
bound — the  right  hand  tide  of  (3.13) — by  either 
analytical  or  empirical  meant.  Once  the  bound  it 
found,  it  remtint  to  generate  reatonable  ettimatea  of 
the  model  error  and  compare  tbit  to  the  bound. 

An  alternative  procedure  it  to  verify  Leaaa  4.1 
by  direct  empirical  meant.  In  other  uordt,  we  can 
utilize  Theorem  3.2  to  give  qnalitative  gnidelinet  on 
the  required  tpectral  character iitict  of  i,  and 
then  timulate  (4.7)  for  a  variaty  of  initial  ttataa 
((0)  and  input*  w  e  L*.  Tbit  latter  approach  it 
not  theoretically  perfect,  hnt  it  a  practical  meant 
to  gain  under ttanding  of  tha  adaptive  ayttem 
behavior. 


(Al)  x#  t  PE  <“  *,  hat  apectral  line*  at 
f requenciet  e>  ,  ....  *  with  rank 
la  (*r  ).  ....a  (ei  )JP«  n. 

*e  1  *t  » 

(A2)  i,T(»)  *  SPR 


Then,  the  adaptive  tyttem_(3.10)  it  locally  L^-etabie 
if.  for  tome  o  t(O.X),  B^ft-o)  ia  atabla  and 
bounded  by. 


|i  (jui-o)l  <  - 

•v  mp  tnp  la 

><  .2  „  .2.1/2 

♦  (X-a)  ] 

lv7v1  • 


(-k)l 


Va  a  B. 


inf 


(3.13) 


with  m.  X  from  (3.7)  and  q  from  (3.9). 


Proof 

See  Appendix  C. 


6.2  Other  Mechanitn*  for  Local  Stability 

Although  wa  have  focated  on  peraiateat 
excitation  at  a  meant  to  nature  local  atahility.  tbit 
it  by  no  meant  the  only  way.  For  example,  if  the 
adaptation  algorithm  (3. Id)  it  modified  to  include  a 
retardation  (tee,  a.g..  (11).  (12))  then  L(t)  in 
(3.10h)  will  have  tha  form  (13): 

L( t)  B  -  B’  >  0  (6.1) 

where  (a.b)  are  poaitivc  coattanta.  Tbit  maaaa 
that  if  a,  it  a  conttant  vector,  than  the 
liaearixed  ayttem  (3.10c)  ia  L^-atable  by  pataivity 
argument*  (13).  (Note  that  it  it  aot  poatible  to 
prove  (4.7)  atable  for  i,  conttant  with  L(t)  ■ 

( 1 / a ) B  a*  in  (3.10h).)  Bence,  ating  (6.1)  together 
with  thaoremt  on  tlowly  varying  ayttamt  (x,  etaya 
near  conttant  long  aaongh),  we  can  arrive  at 
condition*  for  local  Ls-atability  which  are 
independent  of  peraittent  excitation  (tee  (13)  for 
preliminary  reaulte). 


(1)  Theorem  3.2  providea  an  axel  icit  nnner 
hound  on  the  amount  by  which  H  can  deviate  from  a 
nominal  B  which  it  SPB.  Thai,  if  (3.13)  held* 
and  if  II*.  i*  aafficiently  email  (4.2),  then 
aigaelt  ia  the  adaptive  ayetem  are  guaranteed  to  be 
bounded. 

(2)  Unlike  the  global  atahility  aaae  where  the 
bound  on  the  deviation  I  ia  aevarely  reatriatad 
(3.9),  the  bonad  hare  cae*Se  quite  large.  Moreover, 
th*  bound  can  be  determined  from  the  epeetral 
propertiea  of  af.  Recall  from  Lemma  3.1  that  m 
and  X  ia  (3.13)  are  faactioa*  of  the  apectral 
proptriitt  of  t.. 


APPENDIX  A 

ngfiE  at  msaai 

Wa  firtt  thow  that  f(x)  in  (3.10b)  haa  the 
property  that  Va  >  0, 

lx  I  <  a  ■>  I  f  f  a)  I  <  |  Ixl  (A.l) 

From  (3.10b). 

if (x> i  -  ( ir-51*  ♦  i;;i*),/l 

i  i;i  (i5is  ♦  i;i*),/1 
-  m  ( i.is  -  i;i,),/*. 


by  (3.7b) 


i 


|  I  » 1 2 .  by  bolding  |>  I 


by  I* I  i  «  . 

Non.  iiitint  |<nporiil)y  Ikil  (A.I)  bold*  for  *1) 

i  •  B°.  i.*.  . 


I  f ( a  >  |  i  ( a/2)  I *1.  V|a|  4  * 


(A. 2*) 


I  f ( a  )  |  l  <  c/2  >  1*1.  Vial  >  t 


(A. 2b) 


1  i >«  0.10*) 


ll.ll^  i  <l»Llllfc  ♦  IlGfUJll^ 

4  IliJI-p.  ♦  t„  I  lf(i)  1 1^.  by  (4.1) 
4  111,  ll_  ♦  (1  t/2)  I  la )  I  (A. 3) 


itting  the  temporary  utuptioo  (A. 2).  Since  x.  «  L 
(4.2)  *ad  ll*Lll.p.i  llxjJI^, 


Il.ll^  4  ll*Lllm  *  <gm  «/2)  11*11^ 

<  t  (l~tm  t/2)  +  ( g_  t/2)  ll.ll^  (A. 4) 
by  (4.2).  Silica  g_  t/2  <  1  by  attnaption  (4.2), 


IWI*.  i.  « 


and  hence.  Hill  It.  Looking  back  over  the  proof 
•e  tee  that  the  temporary  aaanaption  (A. 2)  it  never 
violated,  i.e.,  the  behavior  of  f(z)  for  1st  >  • 
(A. 2b)  i*  never  needed  under  the  aaaaaptiona  of 
Theorea  4.1.  Thi*  provea  part  (i). 

Part  (ii)  proceed*  analogontly  ezeept  non  we  nae 
the  L  -not*.  Note  that  part  (ii)  aaca  11*11.  <  • 


aa  an  aaanaption. 


APPENDIX  B 


Be  nae  the  exponential  weighting  technique*  froa 
.  (14).  Let  y  denote  the  exponential  weighting 


(9).  (14) 
operation 


(y“)(t)  :«  y*(t)  :«  eat  y(t) 


Jf  y  «  Bu  then  let  B*  denote  the  nap  n  i->  y  . 
For  eaaaple,  if  B  he*  tranafer  function  H(t), 
then  H°  ha*  trantfer  function  B(a-a), 


83B1I  If 


An  operator  1:1,  ->  Ljf  ha* 

v  if  a  nonnegatfve,  nonincrei 


>(.)  *  Lj  each  that 


,  t  ,  _ 

!(la)(t)r  i  J  >(t-t)ln(f)rdt.  Vt  2  0,  Vu  •  L"#, 
•  (1.2) 


Seppote  I  :  if  ->  if  he*  decaying  L.-aeuory. 
If.  for  eoae  a  >  0,  r  L,-*t»ble,  then  *■  c  l,- 

_ . .ii .  A 


•table. 

Apply  Lena*  B.l  a*  follow*:  Th*  aipoaaatial ly 
weighted  vertioa  of  (3.7)  i*. 


U"(t)l  4  ceU  u,‘l(0)  ♦  /  o'»U'oM,',)|.°(t)  « 


(Oct“HT))dr 


Tbit,  together  with  tbe  definition*  of  F  (3.10)  and 
q  (3.9)  give*. 


lll*lln  i  (■/  2(k-o)  >  U(0)  I  ♦  llw-!!, 


♦  (q  ♦  Yj(F)>  lll*lln 


llting  q  *  tj (F)  <  1  froa  (3.3)  give*. 


IU°II~  i  (i*q_y.(F))~1  ( (a/  2(X-o)  U(0)l 

£  ii.‘iin)  ..... 


Bence,  w®  |->  {“  e  1.  -(table.  Moreover,  nting  (B.3) 
w  |->  l  i*  exponentially  (table  and  bence.  froa 
definition  B.l,  w  I - >  (  hat  decaying  L  -aeaory. 
Therefore,  tbe  condition*  of  Leaaa  B.l  are  *ati*ficd; 
contequcnlly  w  |->  (  a  L^-atablc. 


APPENDIX  C 


Be  need  to  calculate  th*  L  -gain  of  F  (3.10) 
under  the  attnaption*  of  th*  thlorea.  Froa  (3.10). 
F  :  x  l->  y  ha*  th*  fora. 


y  «  Fx  »  Gi(  Hi jx 


where  G  and  B  have  trantfer  function*. 


H(c)  -  Hev(»-«> 


Sine*  F  ia  cental,  th*  L^gain  of  F  i*. 


Ilyll, 

VF)  *  tutt: 

xeLVO  * 


•up  - ■ 


xtL./O 

Theorem) 


wbara  y(j«)  end  x(Jw)  are  tha  Fourier  tranafora* 
of  y  and  x,  re»pectively .  Balng  aaanaption  (Al) 
of  Thaoraa  3.2, 


P  P 

y( ja>  -  C  €  CC3ae>»C Jae- 

Ja<k)ei(*ijT*#(r),*(j**-Jak-Jaj) 


l*nca. 


P  P  ,  2/2 

ly(Ja>|  i  e(a)  1C  C  |a( Ja-Ja.-Ja  I*  1 

*■*  -*  <c.»> 


•  (a)  "  (anp|a|  (ak)|)(*ap|G(ja)l(J«- 
Jak)l.|a#k(ak)t)  k  (C.S) 


llyllj  ■  (  /  ly(Ja)lJda) 


.-,4- 


: :  .  •• '  / . : : 


1/2 


•  '  i  1 

l>L,  >  1  -  >11  /  >J  V_  l»(jw  J**,  >  d*l 

u  *  11  I1! 

‘  pl»u|>  c<«)|  ll»llj  <C.9> 

Using  the  definition  of  |iin  (C.5)  together  vitb 
(5.0  from  Theorem  5.1  gives, 

y^tF)  i  P  sup  c(*>)  <  1  *  <1  (C.10) 

v 

Condition  (5.13)  folios*  by  lubttituting  (C.2)  and 
(C.3)  into  (C.10)  nnd  renrianging  terms. 
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This  paper  reports  some  preliminary  results  concerning  robustness 
properties  of  adaptive  control  systems  to  unmodeled  dynamics  and  bounded 
disturbances.  The  analysis  is  conducted  from  the  viewpoint  of  input/output 
stability  theory.  Generic  representations  are  proposed  for  both  continuous¬ 
time  and  discrete-time  adaptive  systems  and  conditions  for  stability  and 
robustness  are  developed  for  each  case.  These  conditions  require  varying 
degrees  of  a  priori  knowledge  about  the  plant,  e.g..  global  conditions  in¬ 
volving  minimal  knowledge  and  local  conditions  involving  more  restrictive 
assumptions . 


1 .  Continuous-Time  Cate 

A.  Globe!  Analysis 

A  large  class  of  continuous- t ime  adaptive  systems  can  be  represented  by 
the  nonlinear  system  (Fig.  1): 

et  =  c*  -  ll(p)vt.  vt  =  t'  nt 

"  G(p)vt  (Sc) 

*t  -  A(4t,<ot),  »t  =  4tct 

where  et<  e*.  v  c  R,  and  n^,  e  R°.  The  operators  H(p)  and 

G(p)  are  proper  rational  functions  with  real  coefficients  in  the  differen¬ 
tial  operator  p,  i.e.,  (px>t  :*  it*  We  will  refer  to  e^  as  the  output 
error .  as  the  pr rameter  error,  as  the  control  error,  as  the 

regressor  and  A(.,.)  as  the  adaptation  gain.  In  general,  only  e^  and 
are  available  as  measurements.  The  parameter  error  nt  :=  nt  -  n* , 
where  is  the  adaptive  estimate  of  the  true,  but  unknown,  parameter  n* . 

The  signals  e*  and  are  referred  to  as  the  tuned  output  error  and 

regressor,  respectively,  meaning  that  these  signals  are  generated  from  an 
'ideal*  system  with  the  desired  parameters  n  *,  i.e.,  =  n*.  Details  on 

the  relation  between  (S  )  and  the  actual  system  (unknown  plant  +  adaptive 

c 

controller)  can  be  found  elsewhere,  e.g..  Ill —131 •  In  general,  the  unknown 
plant  is  imbedded  in  G(p)  and  H(p),  which,  incidentally,  are  also  func¬ 
tions  of  the  true  parameter  n*. 

Since  yt*  as  well  as  the  plant  are  unknown,  it  follows  that  H(p)  and 

G(p)  are  unknown.  Bowever,  in  order  to  establish  conditions  for  stability 

of  (S  ),  it  it  necessary  to  know  something  about  H(p)  and  G(p).  The 
c 

same  remark  holds  for  knowledge  about  the  tuned  signals  e£  and  (*.  The 
following  theorem  gives  conditions  for  global  stability  of  (Sc>.  The  term 
'global*  refers  to  the  intention  of  requiring  minimal,  but  reasonable, 
restrictions  on  H(p),  G(p),  e*  and  (*.  Proof  of  Theorem  1  is  given  in 

[  21. 
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1  . 


'V< 
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Theorem  1:  Global  Stability 

For  the  system  (S  )  assume  that: 

c 

(Al)  The  elements  of  G(p)  are  strictly  proper  and  exponentially 
stable  (all  poles  strictly  inside  the  left  half  plane) 

(A2)  H(p)  is  strictly  positive  real  (SPR),  i.e.,  the  elements  of 

H(p)  are  strictly  proper,  exponentially  stable,  and  Re  [ U ( j to ) ] 
is  positive  for  all  to  e  [0,“). 

(i)  Suppose  that  the  adaptation  gain  is  constant,  i.e.. 


AUt.  to.  )  “Aw.,  A  =  A'  >  0 
t  t  o  t  o  o 


Under  these  conditions,  if  e* ,  e*  t  L„  L  (=>  e*  ->  0),  and 

2  •  t 


«•.  «•  *  ^ 

then: 

(i-a) 

.  n  .  ,n  ^  ,i 

it  t  Lb  ,  ti  ,  it  c  L  2  f\  L 

(i-b) 

e ,  e  e  UH  L  ,  and  e 
i  • 

( i-c) 

v  t  Ljf)  L^.  v  e  L^,  i 

(i-d) 

t.  i  « i£.  *  -  i*.  i 

t  t 


*2  a  c 


Ct  -  0  exp. 

(ii)  Suppose  that  is  persistently  exciting  [4],  i.e.,  9  constants  o^. 


o^,  a ^  >  0  such  that 


•  +  o. 


J  Vi  "  <*2 


I  ,  V  s  e  R 
n 


(ii-a)  Result  (i)  holds,  and  in  addition,  n^,  v.?  ->  0  exp. 
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( 1  i  —  b )  If  the  elements  of  e*.  #•,  4* ,  and  4*  are  all  in  L^, 

then  the  elements  of  n,  n .  n,  e,  i,  v,  v,  4»  and  4 
are  all  in  L  . 

O) 

(iii)  Suppose  that  the  adaptation  gain  is  retarded .  [5],  i.e., 

I n ^ I  <  c .  c  2  In* I 

In'j  I  /c )  2nt )  .  In"t  I  >.  c 

(iii-a)  Result  (i)  holds 

(iii-b)  Result  (ii-b)  holds. 

Remarks 

The  major  difficulties  in  applying  Theorem  1  are  that,  in  the  first 
place.  H(p)  e  SPR  (condition  (A2))  is  an  unlikely  event  in  actual  systems, 
due  to  the  effect  of  unmodeled  dynamics  [6,7].  Secondly,  the  conditions  on 
e*  as  given  in  (i)  are  also  unlikely,  namely  e*  ->  0.  This  condition 
rules  out  the  presence  of  unmeasurable  bounded  disturbances.  Thus,  the  con¬ 
ditions  on  e*  in  (ii-b)  remain  the  only  realistic  case  insofar  as  the 
tuned  signals  are  concerned.  But,  this  raises  another  problem:  ensuring 
that  either  4t  is  persistently  exciting  (2)  or  that  the  adaptation  gain  is 
retarded  (3).  Notwithstanding  the  difficulty  with  the  SPR  condition  on 
H(p) ,  there  are  specific  problems  related  to  (2)  and  (3).  For  example: 

Persistent  Excitation  (PE):  With  bounded  disturbances  (conditions  (ii-b)) 
it  is  not  known  how  to  guarantee  that  4(  *  PE.  Recall  that  4t  i*  gener¬ 
ated  Inside  the  adaptive  loop,  and  thus,  can  only  be  controlled  from  the 
input,  i.e.,  from  either  e*  or,  more  likely,  from  Since  we  do  not 

know  G(p)  and  H(p)  it  is  not  possible  to  conclnde  beforehand  if  4t  s  PE 
even  if  4*  s  PE.  In  the  special  (unrealistic)  case  of  no  nnmodeled 
dynamics  and  no  bounded  disturbances,  e*  *  0,  (ii-a)  holds  and  4*t  a  PE 
*>  4t  s  PE.  Even  though  this  latter  situation  is  easily  ruled  out,  it  cer¬ 
tainly  makes  sense  that  4*t  e  PE  implies  a  'local'  result.  That  is,  with 
certain  suitable  restrictions  on  signal  size  and  so  forth,  the  system  is 


A(4t.  wt)  = 


A  ii> 
o  t 


A  -  (1  - 
o  t 
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robust.  These  arguments  were  formalized  to  some  extent  in  [3]  and  will  be 
slightly  extended  here. 


Retarded  Update:  In  [5]  it  is  suggested  that  the  update  algorithm  be 
retarded  as  given  by  (3).  Likewise,  in  [8],  a  alow  or  'leaky*  integrator  is 
added.  Although  both  these  schemes  (as  well  as  similar  ones)  do  give 
D.l.B.O.  results,  they  both  require  additional  information  abut  the  plant, 
e.g.,  as  in  (3)  an  upper  bound  on  ln*|.  These  results  can  also  be  con¬ 
sidered  as  'local*  results. 

Slow  Variations:  Together  with  a  retarded  update,  another  mechanism  for  en¬ 
suring  B.I.B.O.  stability  is  to  slow  the  variations  in  ^  (see  [3]).  The 
idea  follows  by  examining  the  simple  constant  gain  retarded  algorithm, 

rr  =  A  to  -  an.  (a  >  o) 

t  o  t  t 

*  -Ao$t  H(p)^«t  +  AQ$tet»  -  a  nt 

If  is  constant  then  exponential  stability  can  be  assured  by  direct  LTI 

techniques.  Thus,  if  ((  varies  slowly  enough  with  respect  to  the 
dynamics  of  H(p)  it  is  reasonable  to  expect  a  similar  result.  We  will  ex¬ 
amine  this  more  closely.  However,  control  of  introduces  the  same 

difficulties  as  in  requiring  c  PE,  i.e.,  only  'local*  results  can  be 

obtained. 


B.  Local  Analysis 

The  system  (S  )  can  be  transformed  to  a  more  useful  form  for  local 
c 

analysis: 


«w  <w  A* 

*  *  *L  ”  *NL 
JNL  *  F  f(?) 


<V 


where : 


63 


*  :*  (n,  e,  4)  :*■  (ji-ji*.  e-e*.  4-4*) 


j 

\ 

i 


*L:~  <nI  '  *L*  V  '  ^*n' 

The  system  (S  )  is  obtained  from  (S  )  by  linearization  of  (S  ) 
c  c  c 

e • .  4* ,  and  n*.  resulting  in  the  linearized  perturbation  response 

The  remaining  nonlinear  terois  i.„  are  contained  in  f(z)  where  F 
time-varying  linear  operator.  The  characteristics  of  F.  as  well  as 
of  x^»  depend  on  the  adaptation  gain  and  the  behavior  of  the  tuned 
nals.  e*  and  4*  (see  13]). 

Consider  the  constant  gain  algorithm  (1)  with  a  retarded  update. 

«w* 

2  depicts  the  resulting  system  (S^)  where: 


about 

V 

is  a 
those 
aig- 

Fig. 


« 


L 


1  =  T  >0.  with  (3) 

p+a  o'  \  0.  otherwise 


M  :=  4*  H(p)  4*'  +  e*G(p)4*/ 


(4) 


N  :=  4*  B(p)  +  e*G(p) 

Thus.  Fig.  2  reveals  that  »L  is  the  response  to  p^  :=  (4*e*,  0,  0) 

with  n  t  0,  whereas  x._  is  the  response  to  :*  (0.  4e.  ('n)  with 

n  =  0.  Clearly,  boundedness  of  the  linearized  response  x.  and  stability 
O  L 

of  the  operator  F  reqnire  stability  of  the  map  q,  into  n,  indicated 

in  Fig.  2  by  K(nQ).  It  is  shown  in  [3]  that  stability  of  I(no>:q  |->  n  en¬ 
sures  the  existence  of  conditions  for  local  stability  of  the  adaptive  system 

(S  )  or  (S  ). 
c  c 

Of  particular  interest  is  the  degree  to  which  it  is  possible  to  main¬ 
tain  stability  despite  arbitrary  dynamics  B(p)  and  G(p),  i.e., 

robustness  to  model  error.  Primary  consideration  is  given  to  unmodeled 
dynamics  in  H(p).  Let, 

H(p)  : «=  B(p)  +  AH(p)  (5) 

where  B  denote  the  nominal  dynamics  obtained  under  ideal  conditions, 

*  n* ;  consequently,  we  may  consider  H  to  be  a  fixed  transfer  function 
which  is  independent  of  n* .  All  errors  will  be  lnmped  into  Ay.  The 
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desired  result  of  the  local  analysis  is  to  obtain  a  quantitative  bound  on 
the  worst  case  model  error  for  which  stability  of  (Sc>  **  guaranteed.  We 
will  do  this  by  analyzing  the  stability  robustness  properties  of  the  map 
,  thus,  the  results  obtained  will  only  verify  that  local  conditions 
exist . 


B.l  Local  Stability  bv  Persistent  Excitation 
Assume  that 


4*  e  PE 


H(p)  c  SPR 


(6) 


Under  these  conditions,  it  follows  from  [4]  that  the  system 


4t  = -*?*<*> 


(7) 


is  exponentially  stable,  i.e..  there  exists  constants  m.  A  >  0  such  that 


l*tl  < 


-At 

m  e  |xol 


(8) 


The  following  result  gives  a  coarse  bound  on  the  model  error 


Theorem  2 


The  system  Kn^)  in  Fig.  2  is  L^-stable  if: 


A/m  >  o  :=  lle*|IJIfll.  T.(G) 


(9) 


and 


T.(Ab)  <  (A/m  -  o)/||{*|l*  (10) 

Proof :  Follows  directly  from  small  gain  theory  (see  e.g.  [9]);  details  are 
in  [3]. 
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Remarks :  Although  sharper  bounds  can  be  obtained  [10],  the  significance  of 

Theorem  2  is  that  !!(p)  need  not  be  SPR  if  e  PE.  The  conditions  of 

Theorem  2  can  be  determined  experimentally  by  simulating  K(nQ)  for  a 

variety  of  n  c  R°,  m  c  Ln ,  and  e  PE.  This  procedure  can  only  yield 
o  ®  t 

an  estimate. 


Local  Stabiliti 


Slow  Variations  with  Retarded  Update 


In  this  case  wc  will  assume  that  is  not  PE,  but  varies  slowly, 

in  a  defined  way,  in  relation  to  the  known  dynamics  of  H(p).  Let 
denote  (*  frozen  at  time  t  =  x.  Let  K^(p)  denote  the  linear  time- 
invariant  operator  given  by. 


K  (p)  :=  [I  +  L(p)S  (p)]_1L<p) 
x  n  x 


(id 


where 


Mx(p)  :=  H(p)  (12) 

Let  R^  denote  the  linear  time-varying  operator 

Rx  t*H(pK*'  +  e*G(p)?*'  -  ^(p)  (13) 

The  operator  is  simply  Kn^)  with  M  *  (p ) ,  i.e.,  fixed  at 

(*.  Thus,  represents  the  effect  of  how  far  {•  is  from  other  values 

(*  under  the  dynamics  of  H(p)  and  G(p). 

Suppose  that  l^(p)  is  exponentially  stable,  i.e.,  there  exists  con¬ 
stants  m,  X  >  0  such  that, 

‘  -X(t-s)  . 

|(Ix(p)u)tl  i  »  J  *  Ivjdi.  n  e  La(0,t),  Vt  a  R  (14) 

o 

The  following  result  is  analogons  to  Theorem  2. 


Zfefi£ZJ£-l: 

The  system  Kt>0)  i»  L^- stable  if: 


67 


X/m  >  o  lle*|latIU*|la>  yJO 


■»  sup  yj  4*114*  ‘-tJHt;*)  (IS) 

t>0 


and 


Y„(Ah)  <  (X/m-  o)/ll(*ll*  (16) 

Proof:  Follows  directly  from  small  gain  theory;  details  in  [3], [10]. 

Remarks:  As  in  Theorem  2,  the  conditions  here  for  local  stability  do  not 
depend  on  B  e  SPR.  and  in  this  case  do  not  depend  on  e  PE.  Thus. 
Theorem  3  is  weaker  than  Theorem  2.  The  key  is  to  establish  (14),  i.e.,  ex¬ 
ponential  stability  of  K^(p).  Note  that  with  U  e  SPR.  and  L(p)  given  by 
(4)  with  a  >  0.  (14)  is  established  by  passivity  arguments.  Tighter  bounds 
on  the  norm  operations  can  be  obtained  [10] .  Also,  the  norms  themselves  can 
be  estimated  by  simulating  candidate  actual  systems  (Fig.  2). 

2.  Discrete-Time  Case 


A.  Global  Analysis 

The  discrete-time  version  of  S  is  somewhat  different,  due  to  the  in- 

c 

herent  system  delay  1^1,  The  following  discrete-time  nonlinear  system  is 
representative  of  most  discrete-time  adaptive  systems: 

«t  B  '  Hl*q_1*vl  t  “  B2*q_1*V2  t 
*  *t  "  Gl(ql)vl.t  "  G2(ql>V2.t 
”t  m  <sd  k* 


The  signal  and  operator  distensions  are  the  sane  as  those  defined  in  (S  ), 

-1  -1  c 
with  q  the  backward  shift  operator  (q  x)^  *  xt-l*  *n  lestr>^*  with  • 

unit  delay  (k  =  1),  (S.)  collapses  to  the  form  of  (S  ).  Specifically, 

o  c 

the  unit  delay  adaptive  systen  is: 


H(q_1)vt.  vt 


G(q_1)v 


(Sd,l> 


”t  *  ”t-l  +  A^t'Wt^'  Wt  ^t*t 

In  this  paper  we  will  only  exanine  ^ .  Details  on  ^  can  be  found  in 

110]. 

A.l .Adaptation  Algorithms 

It  is  an  understatement  to  say  that  the  choice  of  discrete-time  algo¬ 
rithms  is  overwhelming.  However,  following  [11], [12]  they  more  or  less 
belong  to  the  following  almost  generic  types: 

Projection 

2  _1 

A(4t.«t)  “  (1  +  **t*  )  “t  (P> 

Recursive  Least  Squares 


A(W  -  Vt 


s;1  -  s;^  +  ,  Se  «  s;  >  0 


(RLS) 


A(5t,«t>  ■  *t«t 


■I1  -  *1-1  4  l«,l2  .  •„  >  0 
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Available  stability  results  have  dealt  almost  exclusively  with  (S.  .  ) 

-1  d,k 

where,  in  the  deterministic  case,  e»  *=  0,  with  either  H  (q  )  «  1  (or 

-1  1  1 

positive  constant)  and  ll^q  )  =  0,  or  vice  versa,  e.g.,  (111.  (The 
stochastic  version  assumes  e*  has  zero  mean  with  bounded  variance,  e.g., 
(121.)  The  following  theorem  extends  the  deterministic  results  to  the  case 
where  e*  e52*  Thus*  **  approaches  zero  asymptotically,  but  is  not  identi¬ 
cal  to  zero.  Proof  of  theorem  2  is  in  (10). 

Theorem  4:  Global  Stability 

For  the  system  (S.  . )  assume  that: 
d,l 

(Al)  The  elements  of  G(q  1 )  are  proper  and  exponentially  stable  (all 
poles  strictly  inside  the  nnit  disc) 

(A2)  H(q  ‘S  is  proper,  exponentially  stable,  and  for  some  constant 

6  >  0. 

I  H(q_1)  -  1  I  1  6  .  V|q|  =  1  (17) 

Under  these  conditions,  if  e*  c  _  ( *=>  e  — >  0)  and  (*  e  J  °  then  using 
adaptation  algorithm  (P) ,  (RLS),  or  (SA)  results  in  e,  v  t  J[  ^  and  (,  n  t  ^ 
provided  that 

b  <  1  (18) 


Remarks 

(1)  Theorem  4  offers  no  more  than  part  (i)  of  Theorem  1  for 
continuous-time  systems,  in  that  it  is  not  possible  to  insure  an  arbitrarily 
large  model  error.  The  bound  (18)  of  6  <  1  is  as  unrealistic  as  the  re¬ 
quirement  that  H(p)  a  SPR  in  Theorem  1  (  in  fact,  H(p)  e  SPR  implies  that  6 

<  1;  see  (21).  Similar  restrictive  results  for  discrete-time  adaptive  sys¬ 
tems  have  been  reported  in  (13)  and  (14). 

(2)  It  can  be  shown  (10)  that  Theorem  2  is  valid  if,  in  (A2) ,  H  is 
either  an  LTI  operator  in  the  sector: 

I  H( q”1)  -  i( q”1)  |  i  *,  Vlql  -  1 
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where 


(19) 


H(q_1)  e  SPK 

or  if  II  -  Fl(q  *)  is  a  slope-restricted  memoryless  nonl  inesr  i  t  y ,  i.e., 

I  (Hv)  t  -  fKq_1)vt  I  <  6  I  vt  I  (20) 

(3)  With  arbitrary  sector  conditions  on  B,  Theorem  2  holds  for  6  <  1 
if  the  adaptation  gain  is  modified,  e.g., 

AUt.  «t)  =  mtu>t.  Bt  -  <1  +  II  (  W\j~in  (2D 

Other  modifications  like  this  can  be  constructed,  provided  m^  satisfies  cer¬ 
tain  conditions,  e.g.,  if  m  is  a  positive  nonincreasing  function,  then 

1  -1 

sector  properties  on  H  apply  to  the  operator  m^Bm^  .  The  required 
properties  of  mt  relate  to  the  noncausal  multiplier  theory  described  in  [9]. 
Picking  the  right  multiplier  -  which  is  only  needed  in  the  proof  of 
stability  -  is  an  artform  akin  to  selecting  a  suitable  Lyapunov  function  for 
a  nonlinear  system.  The  multiplier  requirements  do,  however,  motivate  a 
myriad  of  modifications  to  adaptation  gains  (as  proposed  in  (21)),  for  which 
multiplier  selection  is  more  easily  facilitated,  see  e.g.  [14].  It  is  un¬ 
clear  at  this  time  whether  these  modifications  can  acheive  practical  sector 
conditions  on  B  for  global  stability,  i.e.,  where  6  >>  1. 

B.  Local  Analysis 

Stability  results  dependent  on  persistent  excitation  or  retarded  update 
have  a  more  'local*  character  than  their  discontinuous-time  counterparts, 
and  thus,  have  been  left  out  of  the  global  analysis.  As  remarked  before 
after  Theorem  1,  these  are  the  known  means  to  insure  J[  ^-stability ,  which  we 
have  argued  is  the  case  most  related  to  the  actual  system  environment. 

The  local  stability  analysis  for  continuous-time  systems  can  be 
developed  analogously  for  the  discrete-time  case,  with  only  minor 
modifications.  Thus,  Theorems  2-3  have  their  discrete-time  counterparts. 

One  major  difference,  however,  is  that  the  nonlinear  term  in  (8#)  is  more 
complicated  doe  to  the  complexity  of  the  adaptation  gain  algorithms,  e.g.. 
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ABSTRACT 


An  input-output  approach  Is  presented  for  analyzing  the  global  stability 
and  robustness  properties  of  adaptive  controllers  to  unmodeled  dynamics.  The 
concept  of  a  tuned  system  is  introduced,  i.e.,  the  control  system  that  could 
be  obtained  if  the  plant  were  known.  Comparing  the  adaptive  system  with  the 
tuned  system  results  in  the  development  of  a  generic  adaptive  error  system. 
Passivity  theory  Is  used  to  derive  conditions  which  guarantee  global  stability 
of  the  error  system  associated  with  the  adaptive  controller,  and  ensure 
boundedness  of  the  adaptive  gains.  Specific  bounds  are  presented  for  certain 
significant  signals  in  the  control  systems.  Limitations  of  these  global 
results  are  discussed,  particularly  the  requirement  that  a  certain  operator  be 
strictly  positive  real  (SPR)  —  a  condition  that  Is  unlikely  to  hold  due  to 
unmodeled  dynamics. 
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1.  INTRODUCTION 


1.1  Background 

The  analysis  and  design  of  adaptive  control  systems  has  been  the  subject 
of  extensive  research  in  the  past  two  decades  [1]-[10],  Adaptive  techniques 
provide  a  way  of  handling  plant  uncertainty  by  adjusting  the  cont. oiler 
parameters  on-line  to  optimize  system  performance.  An  alternative  method  for 
handling  uncertainty  is  to  use  a  fixed  structure  controller  designed  to 
provide  acceptable  performance  for  a  specified  range  of  plant  behavior.  In 
principle,  adaptive  controllers  can  provide  improved  performance  compared  to 
fixed  robust  controllers,  since  they  are  tuned  to  the  uncertain  plant. 

However,  adaptive  controllers  sometimes  exhibit  undesirable  behavior  during 
the  tuning  or  adaptation  process.  For  example,  unmodeled  dynamics  can  cause  a 
rapid  deterioration  in  performance  and  even  instability  [11], [12].  This 
problem  is  not  resolved  by  increasing  the  order  or  complexity  of  the  model. 
Since  the  model  of  any  dynamic  system,  by  definition,  is  not  the  actual 
system,  it  can  therefore  be  argued  that  unmodeled  dynamics  are  always  present, 
ad  infinitum. 

The  main  reason  for  these  difficulties  with  adaptive  controllers  seems  to 
be  that  robustness  to  unmodeled  dynamics  was  not  considered  as  a  design 
criterion  in  the  development  of  the  adaptive  control  algorithm.  The  design 
objective  is  global  stability  of  the  closed-loop  system,  e.g.,  [7],  [9]  and 
various  assumptions  on  the  structure  of  the  plant  are  required  to  achieve  that 
objective.  In  particular.  It  Is  necessary  to  assume  that  the  plant  Is  linear 
and  time  Invariant  (LTI),  that  the  relative  degree  of  the  transfer  function  is 
known  as  well  as  the  sign  of  the  high  frequency  gain.  Such  requirements  are 
not  practical  since  real  plants  are  often  nonlinear  and  time-varying  and  can 
be  accurately  represented  only  by  high  order  (sometimes  Infinite  order  [13  j) 
complicated  models. 

The  need  for  robustness  to  plant  uncertainty  Is  not  unique  to  adaptive 
control.  The  problem  of  robustness  Is  ubiquitous  In  control  theory  and  has 
been  studied  In  the  context  of  fixed  (nonadaptlve)  control  [14]-[17],  These 
studies  rely  on  the  Input/output  properties  of  systems,  e.g.,  [18], [19].  The 
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predominant  reason  to  examine  robustness  Issues  In  this  way  Is  that  the 
characterises  of  unmodeled  dynamics,  such  as  uncertain  model  order,  are 
easily  represented.  Lyapunov  theory,  on  the  other  hand,  is  not  well  suited 
for  this  type  of  uncertainty.  Typically,  plant  uncertainty  Is  characterized 
by  assuming  that  the  plant  belongs  to  a  well  defined  set.  For  example,  a  set 
description  of  an  uncertain  LTI  plant  Is  to  define  a  "ball"  In  the  frequency 
domain.  The  center  of  the  ball  is  the  nominal  plant  model,  and  the  radius 
defines  the  model  error.  This  set  model  description  Is  one  type  of  a  more 
general  set  description,  referred  to  as  a  conic-sector  [15].  The  uncertainty 
In  the  plant  induces  an  uncertainty  in  the  input/output  map  of  the  closed-loop 
system  which  can,  again  be  characterized  by  a  conic  sector.  Performance 
requirements  for  the  control  system  can  be  translated  into  statements  on  the 
conic  sector  which  bounds  the  closed-loop  systems,  making  it  possible  to  check 
whether  a  given  design  meets  specifications,  and  providing  guidelines  for 
robust  controller  design. 

In  this  paper  we  use  the  input/output  approach  to  analyze  the  global 
stability  and  robustness  properties  of  continuous- time  adaptive  controllers 
with  respect  to  unmodeled  dynamics  (although  we  consider  only  conti nuous-tlme 
algorithms,  the  input-output  formalism  can  be  readily  extended  to  the 
discrete-time  case).  By  global  we  mean  that  no  specific  magnitude  constraint 
(other  than  boundedness)  is  placed  on  any  of  the  external  inputs  or  Initial 
condi tons.  We  develop  an  adaptive  error  system  of  a  general  form,  by 
comparing  the  actual  adaptive  system  with  a  tuned  system,  l.e.,  the  control 
system  that  could  be  obtained  If  the  plant  were  known.  This  error  system  Is 
similar  to  the  type  used  In  [7], [8]  where  the  tuned  system  error  output  Is 
zero,  due  to  the  assumption  of  perfect  modeling.  By  relaxing  this  assumption 
we  show  that  the  non-zero  outputs  of  the  error  system  are  the  Inputs  to  a 
nonlinear  feedback  error  system  consisting  of  the  adaptive  algorithm  and  two 
feedback  (interconnection)  opera tors, denoted  by  Hev  and  Hzy  . 

An  Important  consequence  of  this  structure  is  that  the  existence  of 
solutions  (e.g.,  tuned  system  performance)  Is  separated  from  the  stabilty 
analysis  (e.g.,  stability  of  the  nonlinear  error  system).  In  general,  the 
adaptation  law  Is  passive;  consequently.  If  Hev  Is  strictly  positive  real 
(SPR),  then  application  of  passivity  theory  [19]-[21],  provides  global 
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L^-stablllty  of  the  map  from  the  tuned  system  output  to  the  actual  adaptive 
system  output,  even  though  the  adaptive  parameters  may  grow  beyond  all 
bounds.  We  provide  other  conditions  (e.g.,  Hzy  stable)  to  Insure  the 

L  boundedness  of  the  adaptive  gains.  Similar  results  are  developed  to 
« 

Insure  L  -stability  of  the  error  system  by  using  an  exponentially  weighted 

OD 

passivity  theory  [19].  These  results  are  summarized  In  Theorems  1A  and  IB. 

As  a  by  product  of  the  Input/output  view  we  also  obtain  specific  bounds 
on  the  l0  and  L  norms  of  significant  signals  in  the  adaptive  system.  The 
results  are  summarized  in  Corollary  1. 

The  results  In  Theorem  1  and  Corollary  1  are  not  essentially  new  (see 
e.g.,  [7], [8]),  although  they  do  provide  some  extentions  to  previous 
results.  The  main  contribution,  however,  is  the  fact  that  all  the  results  can 
be  obtained  from  a  generic  error  system  and  from  the  application  of  nonlinear 
stablity  theorems  based  on  input-output  properties.  As  a  consequence  of  this 
approach,  it  is  to  be  expected  that  conditions  for  robustness  will  arise  in  a 
natural  way.  Such  robustness  results  are  obtained,  but  unfortunately .they 
have  a  limited  practical  use.  The  main  limitation  is  that  the  global  theory 
(Theorem  1)  requires  that  Hgv  e  SPR  ,  which  in  turn  places  an  upper  bound  on 
the  size  of  the  unmodeled  dynamics  in  the  plant.  The  details  are  contained  in 
Lemmas  4.1  and  5.2.  This  bound  is  quite  restrictive  and  is  easily  violated  by 
even  the  most  benign  model  errors,  thus,  verifying  the  results  obtained  in 
[11],  [12].  To  over  come  this  limitation,  we  construct  an  SPR  compensator, 
based  on  the  scheme  proposed  in  [22]  in  the  context  of  robust  (non-adaptive) 
control.  Although  in  the  adaptive  case  the  supporing  arguments  are  heuristic, 
an  example  simulation  shows  a  positive  result. 

The  input/output  analysis  presented  here  provides  a  generic  framework 
within  which  it  is  possible  to  analyze  the  robustness  of  adaptive  robust 
controllers.  We  believe  that  this  framework  can  be  used  to  develop  practical 
adaptive  control  algorithms  that  can  be  more  readily  applied  to  real  systems, 
than  the  class  of  algorithms  currently  In  use. 

Since  this  paper  merges  Ideas  from  several  areas.  It  is  necessary  to 
Introduce  a  number  of  definitions  and  concepts. 


Since  this  paper  merges  ideas  from  several  areas.  It  Is  necessary  to 
Introduce  a  number  of  definitions  and  concepts. 


2.  SOME  PRELIMINARIES 


2.1  Notation 

The  Input/output  formulation  of  multivariable  systems  Is  the  principal 
view  taken  throughout  this  paper  and  the  notation  and  terminology  used  Is 
standard  (see  e.g.  [18], [19]).  The  Input  and  output  signals  are  assumed  to  be 
Imbedded  In  either  the  normed  function  space 

Lp  s  {x  :  [0,-)  ♦  Rn  |  IUII  <  «}  (2.1a) 

or  Its  extentlon 

Lpe  *  {x  :  [0,T]  ♦  Rn|  ||x||Tp  <  T  <  -}  (2.1b) 

The  respective  norms  ||*|lp  and  11*11 jp  are  defined  as  follows: 

| |x| |  *  lim  1 1 x| I-  (2.2a) 

p  T—  p 

with 

/(/  |x(t)|pdt)1/p  ,  pe[l,-) 

1 1 x| |Jp  *  |  °  (2.2b) 

(  sup  I  x(  t)  j ,  p  -  - 
V  te[0,T] 

where  |.|  Is  the  Euclidean  norm  on  Rn.  Hence,  Is  an  Inner  product 
space,  with  Inner  product  <x,y>y  of  elements  x,  y  e  L2e  defined  by 

T 

<x,y>,  s  /  x(t)'y(t)dt  (2.3) 

o 

and  so  ||x|jT2  *  (<x,x>T)^2  •  If  T  ♦  •  then  Is  an  Inner-product  space 
with  Inner  product  <x,y>  *  11m<x,y>  . 


2.2  Stability 


Systems  considered  In  this  paper  are  described  by  Input/output  equations 
of  the  form  y  =  Gu  where  G:L™  ♦  l"  Is  a  causal  map  from  u  Into  y,  also 

pG  p6 

denoted  u  ♦  y  .  The  system  G  Is  said  to  be  Lp-stable  (or  simply  stable)  If  G 

maps  u  e  Lm  Into  y  c  Ln  and  If  there  exists  finite  constants  k  and  b  such 

that  IlGufjjp  <  k  llulljp  +  b  ,  for  all  T  >  0  and  all  UeLpe  •  The  smallest 

k  that  can  be  found  is  referred  to  as  the  Lp-ga1n  (or  simply  gain)  of  G, 

denoted  y  (G)  . 

P 

Because  we  often  encounter  LTI  systems  It  Is  convenient  to  introduce  the 

following  notation.  Let  R(s)  and  RQ( s)  denote  the  proper  and  strictly  proper 

rational  functions,  respectively.  Let  S  and  Sq  denote  functions  in  R(s)  and 

R  (s)  ,  respectively,  whose  poles  all  have  negative  real  parts.  Thus, 
o 

S  and  Sq  are  the  stable,  lumped,  LTI  systems.  Denote  multivariable  systems 
with  transfer  function  matrices,  by  R(s)nxm,  Snxm  .  etc.  For  example, 

G  e  Sjjxm  means  that  all  elements  of  G  belong  to  SQ  ,  and  so  on. 

If  G  €  Snxm  then  the  following  Lp-ga1ns  are  obtained, 


y.(G)  <  y  (G)  *  /  ^lG(t)]dt 

1  *  0 

(2.4) 

y  (G)  =  SUP  olG(ju)] 

^  -  D 

(2.5) 

where  F(A)  denotes  the  maximum  singular  value  of  the  matrix  A,  defined  as  the 
positive  square  root  of  the  maximum  eigenvalue  of  A*A,  where  *  Is  the 
conjugate  transpose  of  A.  In  (2.4),  (2.5)  G  Is  the  operator,  G(ju)  the 
transfer  function  matrix,  and  G(t)  Is  the  Impulse  response  matrix. 

2.3  Passivity 

The  following  definitions  follow  those  In  [19], [21].  Let 
G:Lje  Lje  and  let  p  be  constants  with  u  >  0  .  Then,  V  u  e  Lge  : 


G 


Is  passive  If, 


<  u,  G  u>j  >  p  (2.6) 

G  Is  Input  strictly  passive  If, 

<  u,  Gu  >T  >  p  +  u»u«T2  (2.7a) 

G  Is  output  strictly  passive  If, 

<  u,  Gu  >T  >  p  +  uiGu«y2  (2.7b) 

(y  and  p  are  not  the  same  throughout).  When  G  e  Smxm  satisfies  (2.7),  G  Is 
said  to  be  strictly  positive  real  (SPR),  denoted  G  e  SPRm  .  Because  SPR 
systems  play  a  crucial  role  In  the  proof  of  stability  of  adaptive  systems,  we 
introduce  the  following  subsets: 

SPR™  *  {G  e  S1"*™!*^  [G(  ju)  +  G(-ju)1]  -  yl)  >  0,  VueR}  (2.8a) 

SPR™  =  {G  e  s“|x{£  [G(  ju)  +  G(-ju)']  -  y  G( -jo,)  ’G(  ju))  >  0,  VueR}  (2.8b) 

where  xjA)  denotes  the  smallest  eigenvalue  of  A.  Thus,  whenever  G  e  Smxm  , 
conditions  (2.7)  can  be  tested  In  the  frequency  domain.  Moreover,  SPR^  and 
SPR™  »  respectively,  separate  the  strictly  proper  SPR  functions  from  the 
proper,  but  not  strictly  proper,  SPR  functions.  In  the  scalar  case,  the 
frequency  domain  conditions  simplify  because  x[G(j<*>)  +  G(-ju>)']» 

2  Re[G( ju)]. 

Certain  unstable  systems  In  R(s)mxm  can  be  passive  by  virtue  of  (2.6). 

In  particular,  GeR(s)mxm  Is  passive  If  G(s)  Is  positive  real.  The  transfer 
function  matrix  G(s)  Is  positive  real  If:  (1)  It  has  no  poles  In  Rets)  >  0, 
(11)  poles  on  the  ju  axis  are  simple  with  a  non-negative  residue,  and  (111) 
for  any  u  e  R  not  a  pole  of  G(ju)  +  G(-ju)'  >  . 


2.4  Model  Error 


The  cornerstone  of  robust  control  design  Is  a  quantifiable  bound  on  the 
error  between  the  model  used  for  control  design  and  the  actual  plant  to  be 
controlled.  In  the  adaptive  control  case  considered  here  the  model  Is  a 
parametric  model,  where  the  parameters  are  not  known  exactly.  The  structure 
of  the  parametric  model  can  be  obtained  analytically  from  physical  laws,  but 
this  invariably  results  In  a  complicated  model.  Often  a  simple  structure  Is 
selected  because  it  is  more  convenient  for  analysis  and  synthesis. 


Let  P  denote  the  plant  to  be  controlled.  In  the  broadest  sense  P  Is  a 


relation  in 
(u,y)eLje  x 


Lle  x  Lle  ’  *'e‘*  the  set  possible  ordered  pairs 

L of  Inputs  u  e  and  outputs  ycL^  that  could  be  generated 


by  the  plant  [18].  The  uncertainty  In  the  plant  Is  denoted  by  (u,y)  e  P  . 


Let  p  :L™  ♦  l"  denote  a  parametric  model  of  the  plant  P  with 
a  pe  .  pe 

parameters  a  e  R  •  The  parameters  can  be  selected  so  as  to  minimize  any 
discrepancies  between  the  model  and  the  plant,  l.e., 

,nf  'y-P«u,Tp  '  ,*-P*“'Tp  (2.9 

aeR* 


L 

We  will  refer  to  a+eR  as  the  tuned  model  parameters  and  to  P  *  p,.  as  the 

a* 

tuned  parametric  model  of  the  plant.  In  general,  P*  Is  dependent  on  the 
Input/output  sequence. 


Most  of  the  previous  work  on  adaptive  control  deals  with  the  case  where 
for  every  (u,y)  e  P  there  exists  a  tuned  parametric  model  P*,  such  that 
P**P.  In  this  paper  we  consider  the  presence  of  unmodeled  dynamics,  thus, 
the  uncertain  plant  P  cannot  be  perfectly  modeled  by  any  parametric  model 
P^  .  Since  we  will  deal  exclusively  with  LTI  plants  P  e  R(s)nxm  ,  It  Is 
convenient  to  describe  this  model  error  In  the  frequency-domain.  Let 
Bj(r)  denote  a  "ball"  In  S  of  radius  r,  defined  by 

B$(r)  {6  e  Snxn\  ^G(ju)]  <  r(tt),  H  e  R}  (2.10) 
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Let  the  plant  to  be  controlled  be  described  by 


P  *  (1  +  a)P* 


(2.11a) 


where  p  e  R(s)nxm  Is  the  plant,  p*  e  R(s)nxm  is  the  tuned  parametric  model, 
and  a  e  Snxn  denotes  the  unmodeled  dynamics.  Further,  the  only  knowledge 
available  about  a  is  that  It  Is  bounded  such  that 

a  c  B  («)  (2.11b) 

where  a(u>)  is  known  for  all  frequencies.  In  other  words,  while  the  operator 
a  is  not  precisely  known,  we  do  know  a  bound  on  its  effect.  This  model 
description  (2.2)  is  used  throughout  the  paper  to  precisely  define  the  plant 
to  be  controlled  in  an  adaptive  system.  Following  Doyle  and  Stein  [16]  we 
will  refer  to  (2.11b)  as  an  unstructred  uncertainty.  Note  that  although  a  is 
stable,  P  and  P*  need  not  be  stable.  Hence,  the  parametric  model  is 
implicitly  required  to  capture  all  unstable  poles  of  the  plant.  Although  this 
is  not  severly  restrictive  -  at  least  on  practical  grounds  -  nonetheless,  it 
can  be  eliminated  by  definng  model  error  as  (stable)  deviations  in  (stable) 
coprime  factors  of  the  plant  [23].  As  the  subsequent  analysis  is  not 
substantially  effected  by  this  choice,  we  will  remain  with  (2.11)  for  purposes 
of  illustration. 

2.5  Persistent  Excitation 

From  [31],  a  regulated  function  F(.)  *  R+  ♦  Rnxm  is  persistently 
exciting,  denoted  F  e  PE  ,  if  there  exists  finite  positive  constants  , 

,  02*  and  such  that 
s+o3 

<*„  I  >  f  F(t)F(t)'dt  >  a,  I  ,  V  s  e  Ra  (2.12) 

2  n  In  + 

The  usefulness  of  a  persistently  exciting  signal  is  in  establishing  the 
exponential  stability  of  the  following  differential  equation  which  arises  in 
many  adaptive  and  identification  schemes,  l.e., 


x  *  -BFHF'x  ♦  w  , 


(2.13) 


* 


1 


It  Is  shown  In  [31]  that  If  Be  Rnxm,  B  *  B*  >  0  ,  He  SPR™  or  SPR+,  and 
F  e  PE  ,  then  (w,  x( 0 >  | — ►x  Is  exponentially  stable,  1.e.,3m,  x  >  0  such 
that 

I x( t) |  <  me'Xt  | x ( 0 )  |  ♦  J  me‘x(t"r)  |w(T)|dT  •  (2.14) 

0 

We  will  utilize  this  latter  result  In  section  IV  in  our  proof  of  stablity  of 
the  adaptive  system. 


■ 

0 


k 
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3.  ADAPTIVE  ERROR  MODEL 

In  this  section  we  develop  a  generic  adaptive  error  model  which  will  be 
used  In  the  subsequent  analysis.  This  requires  defining  the  notions  of  robust 
control  and  tuned  control . 

Robust  and  Tuned  Control 

Consider,  for  example,  the  model  reference  adaptive  control  (MRAC) 
depicted  In  Figure  3.1,  consisting  of  the  uncertain  plant  P,  a  reference  model 

A  A 

Hr,  and  an  adaptive  controller  C(e)  ,  where  e  Is  the  adaptive  gain  vector,  r 
Is  a  reference  input,  d  is  a  disturbance  process,  and  n  is  sensor  noise. 

Denote  by  H(e)  the  closed-loop  system  relating  the  external  Inputs  w  =  ( r * , 
d‘,  n')'  to  the  output  error  e,  as  depicted  in  Figure  3.2..  Also,  let  w  e  W 
denote  the  admissable  class  of  Input  signals. 

The  objective  of  the  adaptive  controller  is  twofold:  (1)  adjust  e  to  a 
constant  9*  e  Rk  such  that  H(e*)  has  desireable  properties;  and  (2)  during 
adaptation,  as  0  is  adjusted,  the  error  is  well  behaved.  In  the  usual 
formulations  [7]  only  (1)  is  considered  and  further  it  is  assumed  that  there 
exists  a  matched  gain,  denoted  by  F  e  Rk  ,  such  that 

H(f)  =  0  (3.1) 

The  presence  of  uncertain  unmodeled  dynamics  in  the  plant  eliminate  the  chance 
of  satisfying  the  matching  condition.  Thus,  it  Is  more  appropriate  to  define 

L- 

a  tuned  gain,  denoted  by  0*  c  R  ,  corresponding  to  each  (u,y,w)  e  P  x  W  , 
such  that 

H(q*)w  <  H(e)w  ,  V  0  c  Rk  (3.2) 

The  error  signal  e*  :*  H(e*)w  is  referred  to  as  the  tuned  error.  Note  that 

each  (u,y,w)  e  P  x  W  engenders  a  possibly  different  0*  .  Also,  It  Is 

Important  to  distinguish  the  tuned  gain  e.,  ,  from  the  robust  gain  e.  c  Rk  . 

*  0 

where 

sup  H(e.)w  <  sup  H(e)w  ,  V  e  c  Rk  (3.3) 

P  x  W  0  P  x  W 
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The  error  signal  e  :=  H(e  )w  Is  referred  to  as  the  robust  error.  It  follows 

0  0  " 

from  these  definitions  that  the  tuned  error  Is  always  smaller  In  norm  than  the 
robust  error,  thus  V  w  e  W  , 

e*  =  H(e*)w  <  en  -  H(e  )w  ,  (3.4) 

The  tuned  controller  Is,  unfortunately,  unrealizable  since  it  requires  prior 
knowledge  of  the  actual  system  H(e)  (or  equivalently,  the  plant  P)  and  the 
input  w.  A  practical  adaptive  controller  is  likely  to  have  a  larger  error 
norm. 

Structure  of  the  Adaptive  Control 


In  summary,  we  consider  the  multivariable  adaptive  system,  shown  in 
Figure  3.2,  and  described  by 


e  =  H(e)w  .  (3.5) 

where  e(t)  e  Rm  is  the  error  signal  to  be  controlled,  w(t)  e  Rq  is  the 

«  k 

external  input  restricted  to  some  set  W,  and  e(t)  e  R  is  the  adaptive 
gain.  The  class  of  adaptive  controllers  considered  here  are  such  that  the 

k 

adaptive  gains  multiply  elements  of  internal  signals  z(t)  e  R  ,  referred  to 
as  the  regressor,  to  produce  the  adaptive  control  signals, 

fi  *  01  z1  ,  1  e  [1 ,m]  (3.6) 

where  and  z^  are  k^ -dimensional  subsets  of  the  elements  in  e  and  z, 
respectively.  Thus, 


m 

k  *  z  k.  (3.7) 

1*1  1 

Define  the  adaptive  gain  error. 


0 ( t)  :*  e(t)  -  8*  (3.8) 

where  8*  e  R  Is  the  tuned  gain  (3.4).  Also,  define  the  adaptive  control 
error  signals. 


Vj  :*  sj  z^  ,  1  =  1,  ....  m 


(3.9) 


An  equivalent  expression  Is, 


v  =  z'e 


(3.10a) 


where  the  time-varying  matrix  Z  is  defined  by 


Z  =  block  diag(zlf  Z2,  .  .  . ,  z^) 


(3.10b) 


To  describe  the  relations  among  the  signals  e,  z,  v,  and  w  we  Introduce 
the  interconnection  system  Hj  :  (w,v  )  ♦  (e.z)  ,  as  shown  In  Figure  3.3.  In 

particular,  let  Hj  e  R( s) (m+k) x(m+q)  ^  and  where  Hj  1s  defined  by. 


:=  H. 


Hew\  /’Hev 


H  /  \-H 
zw  '  zv 


(3.11) 


In  effect,  this  structure  serves  to  isolate  the  adaptive  control  error  v,  from 
the  rest  of  the  system.  When  the  adpative  control  Is  tuned,  9=0  and  v  *  0; 
consequently,  the  tuned  error  signal  (3.4)  is. 


e*  :s  =  Hp  " 


(3.12) 


We  can  also  define  a  tuned  regressor  signal, 


z„  :=  H  w 
*  zw 


(3.13) 


In  general,  all  the  subsystems  In  Hj  are  dependent  on  the  tuned  gains  e*  . 
The  interconnection  system  can  also  be  written  as, 
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i  •  7  •  Mi  i  ■  s  *  *  ■  r* 


-»  *v  - . 


t* 


e  *  e*  -  H  v 
*  ev 


2  *  2*  -  H2V  V 


(3.14a) 


(3.14b) 


with  v  given  by  (3.10).  To  complete  the  error  model  requires  describing  the 
adaptatlve  algorithm.  I.e.,  the  means  by  which  e(t)  Is  generated.  Me  will 
consider  two  typical  algorithms.  A  constant  gain  (gradient)  algorithm  [7]: 


e  =  r  Z  e 


(3.15) 


where  r  e  Rkxk,  r  *  r‘  >0  ,  and  a  similar  but  nonlinear  gain  algorithm: 


e  =  r(Ze  -  p(e)e) 


(3.16a) 


where  p  :  RK  ♦  R+  is  a  retardation  function,  whose  purpose  is  to  prevent 
e  from  growing  too  quickly  in  certain  situations.  Although  many  functions 
will  suffice  we  will  select  the  one  proposed  in  [24],  namely: 


(isi/c  -  1)  ,  lei  >  c  :*  maxie*i 


p ( 0 )  :  = 


(3.16b) 


1 01  <  c 


The  complete  adaptive  error  system.  Is  shown  in  Figure  3.4.  Note  that 

the  error  system  Is  composed  of  two  subsystems:  a  linear  subsystem  and  a 

non-linear  subsystem  eu  . 

N 


i 


i 


4 
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4.  CONDITIONS  FOR  GLOBAL  STABILITY 


The  theorems  stated  below  give  conditions  for  which  the  adaptive  error 
system  (Fig.  3.4)  Is  guaranteed  to  have  certain  stability  and  performance 
properties.  Proofs  are  given  In  Appendix  A.  Heurlstlcally,  however,  the 
basis  for  the  proofs  Is  application  of  the  Passivity  Theorem  ([19],  pg.  182). 
It  turns  out  that  the  map  e  +  v  Is  passive.  Thus,  If  Hey  Is  SPRm  ,  then 
the  map  e*  ♦  (e,v)  is  L2_stable  even  though  2  and/or  9  can  grow  without 
bounds.  Further  restrictions,  provided  below,  cause  9  and  z  to  be  bounded. 
(We  use  the  notation  "x  ♦  0  (exp.)"  to  mean  that  x(t)  ♦  0  (exponentially)  as 
t  ♦  •  •) 

Theorem  A:  Global  Stability 


For  the  adaptive  error  system  shown  in  Figure  3.4,  assume  that: 

(Al)  The  system  Is  well -posed  In  the  sense  that  all 
Inputs  w  e  W  produce  signals  e.v.z,  e  ,  and 


9  In  L  . 

•e 

(4.1a) 

( A2) 

u  ckxm 

zv  e  So 

(4.1b) 

(A3) 

hcv  c  spr; 

(4.1c) 

Under  these 

conditions: 

(1) 

If  (e*.  e*)  e  L*  Hl*  (-►  e*-^0) 

and  (2*,  2*)  e  L*  then  with 

algorithm  (3.15)  or  (3.16): 


( 1  -a ) 

<0.i)  c  l*. 

e  c 

L2  OLi  »  and  8  “►  0. 

(4. a) 

( 1-b) 

e  c 

c 

Lm,  and  e-e*  -►0. 

(4.2b) 

(1-c) 

.in.  ,m 
v  c  l9  n  L 

• 

.  V 

c  L™  ,  and  v  0  . 

(4.2c) 

(1-d)  (z,i)  e  t*  .  (z-z„,  £-£*)  €  LjnL* 


0  exp. 
(4. 2d) 


I 


m 


■ 

*  < 

a 


* 


,  and  z-z*  — ► 


( 1-e) 


If,  In  addition,  e*  *  0  (matched)  and  z*  e  PE  then 
(9,  9,  e-e*.  v,  z-z*)  —►0  exp. 

(4.2e) 


(ID  If  (e*.  e*)eL*  and  (z*,  z*)  e  L*  »  then  w1th  al90r1thm  (3.15): 


(1i-a)  z  e  Lk  (4.3) 

flp 

(11-b)  With  the  addition  of  either  algorithm  (3.16)  or  z  e  PE  It  follows 
that  the  elements  of  e,  9,  e,  e,  v,  v,  and  £  are  In  L  . 

(4.4) 


Theorem  IB:  Global  Stability 

Replace  (A3)  In  Theorem  1  by 

(A3)*  Hev  E  SPR0  (4*5) 

(1)  If  (e*.  e*)  e  (->  e*-^0)  ,  and  (z*.  2*)  c  L*  then  with 

algorithm  (3.15)  or  (3.16) 


( 1-a) 

(9,  9  )c  l!  ,  9  c  ,  e-^0 

(4.6a) 

(1-b) 

e  .  e  e  l"  ,  e  -  e*-^0 

(4.6b) 

(1-c) 

(v.v)  c  i; 

(4.6c) 

(1-d) 

(z,z)  e  l*  ,  (z-z*,  £-£*)  E  L*  . 

and  z-z+ — ►  0  . 

(4.6d) 

(1-e) 

If,  In  addition,  e*  ■  0  (matched)  and  z*e  PE  , 
then  (e,  v)— *>0  exp. 

(4.6e) 

(4.6e) 


(11)  If  (e*,  e*)  c  L™  and  (2*,  2*)  e  L*  .  then  with  algorithm  (3.15): 

(11-a)  2  c  l||  (4.7d) 

(11-b)  With  the  addition  of  either  2£PE  or  algorithm  (3.16),  the 
elements  of  e,  0,  e,  e,  v,  v  ,  and  2  are  In  L^. 

(4.7b) 


Corollary  1:  Performance  Bounds 

Suppose  2*  and  e*  satisfy  the  conditions  In  (1)  of  Theroems  1A  or  IB. 


(I)  Let  Hev  e  SPRj  ,  i.e.,  3  u.  Y  >  0  such  that  VweR  , 

o[Hev(jw)]  <  y  and  ^Hev(ju,)  +  Hev(-ju») 1  ]  >  u  Iffl  (4.8a) 

Then,  bounds  on  iei„  and  1 ei  can  be  obtained  from: 

1 e-e*i 2  <  Jjp  [ie*i2  +  (te*!*  +  2v  e(0)'  r’1e(0))1/2]  (4.8b) 

10 ,r~1ei—  <  e(0)'  r-1e(0)  +  2iei2ie-e*i2/Y  (4.8c) 

(II)  Let  Hgv  e  SPRq  ,  I.e.,  3W,  q,  k  >  0  such  that  V  u  e  R  , 

^{Hev(ju)  ♦Hev(-ju)']  >y  Hey(-ju)'  Hev(jM)  (4.9a) 

|{G  (ju)  ♦  Gev(-ju.)')  >  k  Im  (4.9b) 

Geyfs)  :e  (1  ♦  qs)  H^s)  (4.9c) 

Then,  bounds  on  ie_i  and  iei  can  be  obtained  from: 

L  m 

«e|2  <  ^r^{ie*+qe*i  ♦  (ie*+qe*i2  ♦  2k2we(0),r‘1e(0))1/2]  (4.9d) 

««,r"le»-  <  e(0)'r"1e(0)  +  ^  ie*+qe«i2iei2  (4.9c) 
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Discussion 


(1)  Theorems  1A  and  IB  give  conditions  under  which  the  adaptive  error 

system  is  globally  stable.  Essentially,  conditions  are  Imposed  on  the 

Interconnection  subsystems  In  Hj  .  In  particular,  e  SPRm  and 

H  e  S*xm  are  direct  requirements,  whereas  the  restrictions  on  the  tuned 
zv  0 

signals  e*  and  z*  ,  Indirectly  Impose  requirements  on  Hew  and  H  .  These 
latter  requirements  are  dependent  on  knowledge  about  w  e  W  .  For  example.  If 
w  Is  a  constant,  then  the  assumption  that  e*  0  (Theorem  lA-i)  requires 
that  the  tuned  feedback  system  Is  a  Type-I  robust  servomechanism,  l.e.,  the 

transfer  junction  Hew(0)  =  0  for  all  (u,y)  e  P  . 

(2)  Corollary  1  gives  explicit  bounds  on  signals  In  the  error  system. 

These  bounds  can  be  used  to  evaluate  the  adaptive  system  design.  Moreover, 

the  bounds  allow  a  coarse  determination  as  to  the  efficacy  of  adaptive  control 

vs.  robust  control.  By  comparing,  for  example,  the  adaptive  error  le^  from 

(4.8)  with  the  robust  error  ie  i.  from  (1.5),  It  Is  possible  to  obtain  a 

o  2 

quantifiable  measure  of  performance  degradation  during  adaptation. 

(3)  Although  Theorems  1A  and  IB  are  essentially  the  same,  there  are 
slight  difference  worth  noting.  These  differences  arise  because  In  1A, 

HeveSPR™  Hev(s)  is  proper  but  not  strictly  proper,  whereas  In  IB, 
HeveSPRj=->  Hgv(s)  Is  strictly  proper.  Thus,  comparing  part  ( 1 )  In  1A  and 
IB,  we  see  that  In  IB,  v,  v  e  lm  whereas  In  1A,  v  Is  addl tonally  In 

m 

l*  and  v-*-0  . 

(4)  The  use  of  persistent  excitation  or  gain  retardation  Is  seen  In  part 
(11)  of  theorems  1A  and  IB  to  provide  the  means  to  guaranty  bounded  signals. 
Other  schemes  based  on  signal  normalizations  or  dead-zones  can  provide  similar 
results,  e.g.  [32], [33].  The  effect  of  these  conditions  Is  to  provide  an 

L^-stablllty  which  Is  not  present  otherwise.  The  persistent  excitation 
condition  actually  supplies  exponential  stability,  which  is  stronger  than 
L^-stablllty,  as  provided,  for  example,  by  the  gain  retardation  (see  proof  In 
Appendix  A). 

(5)  The  persistent  excitation  requirements  In  parts  (1)  and  parts  (11) 
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are  different.  In  parts  (1),  z*ePE  ,  whereas  In  parts  (11),  ZePE  .  The 
different  assumptions  arise  because  In  parts  (1)  we  enforce  the  matched 
condition  e*=0  .  Hence,  z*ePE  =>  ZePE  .  This  follows  from  (1-d)  where 
z  -  z^  ♦  0  expoentlal ly.  Also,  with  e#  *  0  ,  a  bounded  disturbance  added  to 
the  reference  can  cause  z  e  PE  without  forcing,  e.  e  L  .  In  parts  (11), 

W  OD 

which  Is  more  realistic,  we  disallow  the  matched  condition,  and  hence, 
e*  e  L  .  Thus,  z  c  PE  Is  the  weakest  assumption  to  make.  However,  since  z 
Is  inside  the  adaptive  loop,  it  is  very  different  to  guarantee  z  e  PE  by 
injecting  external  signals.  Note  also  (In  both  parts(ii))  that  without 
retardation  or  PE  it  Is  possible  for  the  regressor  to  remain  bounded  even 
though  the  adaptive  parameters  may  grow  unbounded.  Similar  results  have  been 
reported  elsewhere,  e.g.  [24]. 


Robustness  to  Unmodeled  Dynamics 


Since  the  theorems  Impose  requriements  on  the  input/output  properties  of 
the  Interconnection  system.  It  follows  that  the  effect  of  model  error  on  these 


properties  determines  the  stability  robustness  of  the  adaptive  system.  For 

example,  both  theorems  reaulre  that  Hgv  e  $PRm  •  Suppose,  however,  that 

H  has  the  form, 
ev 


VHe, 


(4.10) 


where  flev  Is  the  projection  onto  Hgv  of  the  plant  uncertalny  operator  a  ; 

and  H  Is  the  nominal  transfer  function  when  there  is  no  uncertainty,  i.e., 
ev  _ 

when  a  =  0  .  Thus,  H  is  a  function  of  the  tuned  parametric  model  P*  and 

ev 

the  tuned  controller  gains  .  (See  Section  V  for  more  specific  formulae, 
e.g.  (5.5).) 


Conditions  to  Insure  that  H  .  e  SPR™  despite  uncertainty  In  H~v  Is 

ev  + 

provided  by  the  following: 


Lemma  4.1;  Let  Hev  be  given  by  (4.3).  Then  H@y  e  SPR+  If  the  following 
conditions  hold: 


(1)  He  SPR* 
ev  ▼ 


(4.11a) 
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(4.11b) 


(11)  flevc  B^(k)  where  V  u»  e  R  , 

Mu.)  <yX|»ev(Ju,)  +  Hev( -J« ) 1  ] /7[ Hev( J») ]  O.llc) 

Proof:  Define  ^( . ) :  Cmxm  -*•  R  by 

L  (A)  =  j  x_( A+A*) 

where  *  denotes  conjugate  transpose.  Then,  using  definition  (2.8)  with  (4.10) 
-  (4.11)  we  obtain 

MlHev(ju)]  s  ElHev(ju)  +  ^ev(j“^Hev( 

>  it!  Hpy  (<*«*>)]  -  >  0  . 


Hnece,  Heye  SPR*  • 
Comments 


(1)  In  order  to  apply  Lemma  4.1  It  Is  necessary  to  have  a  detailed 
description  of  how  the  plant  uncertainty  a  propagates  onto  the 
Interconnection  uncertainty  fley  .  This  type  of  uncertainty  propagation  was 
explored  In  depth  by  Safonov  [25]  and  more  sophisticated  expressions  then 
(4.4b)  are  available  to  describe  the  uncertain  operator  Hev  .  Section  5 
contains  more  detail  on  this  Issue. 


(2)  In  the  scalar  case  (4.11c)  becomes 


Mu.)  <  Re[Hev(ju)]/|Hev(J»)| 

■  COS  *  [HgyUw)] 


(4.12) 


Since  H  e  SPRm  by  assumption.  Mu)  Is  always  positive  for  u  c  R  ;  but 
ev 

because  of  the  cosine  function,  k(u)  <  1  •  In  Section  6  we  show  that  this 
limitation  on  the  effect  of  model  error  Is  easily  violated  by  even  the  most 
benign  type  of  unmodeled  dynamics  In  the  plant.  Methods  which  overcome  this 


limitation  are  discussed  In  Section  7.  The  requirement  that  k(u»)  <  1  also 
holds  for  any  multivariable  Hgv  e  SPRm 
decomposition. 


To  see  this  let 


H"  have  the  polar 
ev 


H  !  G  W  =  W  G 
ev  t  ev  ev  r 


(4.13) 


where  G„  ,  Gr  are  Hermltlan  and  V/ev  Is  unitary. 


_  l 

o(H  )  =  cT(G  )  *  o(G  ) 
ev  »  r 


Since 


it  follows  that 


Mw)  <  o[wevUu)]  <  1 


(4.14) 


In  the  case  of  scalar  systems,  the  condition  k(w)  <  1  can  be  interpreted  In 
terms  of  a  limitation  on  relative  degree  of  Hey(s)  .  A  necessary  condition 
for  H  e  SPR  is  that  the  relative  degree  of  H  (s)  does  not  exceed  one 
l.e.,  phase  limited  to  ±90°.  Rohrs,  et  al .  [12]  show  that  this  necessitates 
precise  knowledge  of  plant  order,  and  hence.  Is  not  a  feasible  requirement  In 
the  presence  of  an  unstructured  uncertainty  (2.12),  where  the  order  is 
unknown.  In  the  multivariable  case  it  is  awkward  to  talk  about  relative 
degree  or  phase,  however,  (4.14)  expresses  the  same  limitation. 


(3)  In  several  Instances,  e.g.,  [9] ,[26], [27],  it  has  been  reported  that 

the  SPR  condition  has  been  eliminated.  In  each  case,  however.  It  can  be 

verified  that  the  operator  Hey  =  positive  constant  ,  which  is  SPR.  But, 

these  studies  do  not  account  for  unmodeled  dynamics,  thus,  in  the  notation  of 

(4.10),  only  IT  *  positive  constant  .  Lemma  4.1  then  provides  the  means  to 
ev 

evaluate  the  effect  of  unmodeled  dynamic. 


5.  APPLICATION  TO  MODEL  REFERENCE  ADAPTIVE  CONTROL 


Consider  the  model  reference  adaptive  control  (MRAC)  system,  shown  In 
Figure  5.1,  consisting  of:  an  uncertain  scalar  plant  P  e  Rq(s)  ;  a 

reference  model  e  ;  and  filters  with  F  e  .  The  plant  is 

affected  by  a  disturbance  d  and  a  reference  command  r.  The  system  eouations 
are; 


e  =  y  -  yr  (5.1a) 

yr  =  Hr  r  (5.1b) 

y  =  d  ♦  Pu  (5.1c) 

u  =  -  0'2  =  -(efzj  +  e2z2)  (5. id) 

Zj  =  F  u,  z?  =  F(y-r)  ( 5 . le) 

Assume  that  the  adaptive  law  is  given  by  (3.15),  thus, 

e  *  r  z  e  (5. if) 


Let  the  plant  uncertainty  be  described  by{2.12),  i.e.. 


A  := 


c  Bs(6) 


(5.1g) 


where  P.  e  R  (s)  is  a  tuned  parametric  model  for  P.  Let  the  filter  dynamics 
o 

be  given  by 


F  (  s )  * 


(5.1h) 


where  L(s)  is  a  stable  monlc  polynomial  of  degree  t  .  Thus, 
flj(t),  e2(t)  e  Rl  and  so  e(t)  e  R2t  .  Using  the  definition  of  tuned  gain 
(3.2)  we  get. 


It  (O 


v  :*  s 1 z  from  (3.6) 


u 


-  9 'z  *  -(e*+a)'z 

■{0*121  +  0*2Z2)  "  v 
^*i  ^*2 

-  u  +  -p-  (r-y)  - 


t 


V 


Finally, 


u 


1 

*  Tnq7T  (r-y)  -  T7S7J7T  v  :=  c*,r'y) 


1 

T+a7[7l 


v 


(5.2) 


whPTA  & 

and  A*2  are  polynomials,  each  of  degree  t-1  ,  whose  coefficients  are 
the  elements  of  the  tuned  gains  a*j  and  8*2  »  respectively;  and  C*  denotes 
the  tuned  controller.  The  tuned  system  (  0*0  )  is  shown  In  Figure  5.2. 


In  terms  of  the  uncertain  plant  P,  the  adaptive  error  system  (Fig.  3.4) 
corresponding  to  this  MRAC  system,  has  tuned  signals: 


e*  =  (1  +  PC*)-1d  +  ((l+PC*)"1PC*-Hr]r 
F(l+PC*)’1C*(r-d) 

Z*S 

F(1+PCJ  !(d-r) 
and  Interconnections: 

Hev  *  (1+PC*)’1P(1+A*1/L)"1 


(5.3a) 


(5.3b) 


(5.3c) 


Fd+pc^rMi-H^j/ir1 

F(1+PC*)"1P(1+A*1/L)'1 


(5.3d) 


The  error  system  can  also  be  described  so  as  to  highlight  the  model  error 
a  .  The  following  definitions  are  convenient: 

T.  :•  (WP.C.r^.C.  1  -  S.  (5.4.) 
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*1 


K*  H 


(5.4b) 


ev|As0  *  u^r^a-vL,-1 

Thus,  the  error  system  (5.3)  can  be  also  be  expressed  as: 

e*  =  S*(l+aT#)_1d  +  (T<,(l+a)(l+aT*)"1-Hr)r  (5.5a) 


F  S^C*(l+AT*)’1(r-d) 

**  =  .j 

_F  S*(l+aT*)  i  d-r) 
Hev  •  K^d+aid+aT*)"1 


)5.5b) 


(5.5c) 


F  K^d+aT*)'1 
F  K*(1+a)(1+aT*) 


(5.5d) 


The  result  that  follows  In  Lemma  5.1  gives  conditions  under  which 

H  e  SPR„  and  H  e  S2*xl  ,  despite  model  error;  thus  conditions  (A1)-(A3) 
ev  o  zv  o 

of  Theorems  1A  and  2B  are  satisfied.  Additional  requirements  are  necessary  to 
establish  the  class  of  tuned  signals  e*  and  z*  as  given  by  (5.5a)  and  (5.5b), 
respectively.  These  requirements  are  discussed  following  Lemma  5.1. 


Lemma  5.1:  For  the  adaptive  system  (5.3)  or  (5.5)  H  e  SPR 
. — .  ev  o 

H  e  IF  the  following  conditions  are  all  satisfied: 

zv  o 


(i) 


P*(s) 


g( sn_1+  ...  +  Bn-1) 

rT  n-1 ^  7 

S  +aj  S  ♦  . . .  ♦  an 


9N.(s) 

u^sr 


and 


(5.6a) 


(11)  N*(s)  Is  a  stable  monlc  polynomial 


(5.6b) 


(111)  g  >  0  (5.6c) 

g  Kj(s) 

dv)  K*(s)  ■  -^-(sy  c  SPR0  where  Kj(s)  and  K2(s)  are  monlc  stable 
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polynomials. 

(v)  i  =  deg  l(s)  >  n  +  deg  Kj(s)  -  1 


( 5.6d) 
(5.6e) 


( vi ) 


a  e  B^( 6 )  Is  such  that 

6(h))  <  6(a»)  :®  n(u>)[n(u)  |T*(  jw)l  +  |S*(ju>)|]  ^ 

V  w  c  R, 

n(u)  :=  cos  *  [Kjju)] 

(5.6f) 


Proof:  See  Appendix  B. 


Discussion 


(1)  Condition  (1)-(v)  of  Lemma  5.1  are  restatements  of  known  results, 
but  normally  they  apply  to  the  actual  plant  P,  e.g.  [7].  In  Lemma  5.1, 
however,  these  conditions  apply  to  the  parametric  model  P*  —  not  to  the 
actual  plant.  As  such,  they  are  easier  to  satisfy,  since  the  parametric  model 
Is  somewhat  arbitrary.  This  flexibility  Is  penalized  by  an  increase  In  model 
error.  For  example.  If  the  actual  plant  has  a  relative  degree  of  2,  then 
choosing  a  parametric  model  of  relative  degree  1  —  as  required  by  condition 
(1)  —  Incrases  the  high  frequency  model  error. 

(2)  Condition  (vi)  Imposes  an  upper  bound  T  on  the  model  error 
associated  with  the  chosen  parametric  model.  This  condition  simultaneously 
Insures  that  Hgv  e  SPRQ  despite  model  error,  and  that  the  tuned  system  Is 
stable  (see  proof  In  Appendix  B). 

(3)  It  Is  easily  verified  that  T(u>)  <  1  ,  as  was  discussed  following 
Lemma  4.1.  In  fact,  even  the  "optimally  tight"  bound  (see  [25]  for  details  on 
this  calculation)  given  by, 

T  ‘  TfTTn”  H1’11  ♦  (ll+T|2  ♦  4n  Re(KT/|K|  )1/2]  (5.7) 


Is  also  restricted  to  be  less  than  1.  This  limitation  severely  restricts  the 
type  of  admlssable  model  error.  This  Issue  Is  pursued  In  Section  6. 


(4)  To  guarantee  global  stability  using  the  adaptive  law  (5. If), 
property  (1)  of  Theorem  1  requires  that  e*  ♦  0  and  z*.  i*  e  L^*  for  all  r 
and  d.  For  example,  let  r  and  d  be  any  bounded  signals  such  that 

r  *  constant  and  d  constant  as  t  ♦  •  .  Property  (1)  of  Theorem  1  is 
satisfied  if: 

6(0)  =  0  (5.8a) 

T*{0)  =  Hf(0)  =  1  (5.8b) 

Zero  model  error  at  DC  (5.8a)  is  certainly  to  be  expected  from  even  the  most 
crude  tuned  parametric  model . 

(5)  Let  r  be  bounded  such  that  r  ♦  constant  as  t  ♦  •  ,  but  let  d  be 
just  bounded,  i.e.,  d  e  L^  .  In  this  case  it  is  not  possible  to  guarantee 

e^  ♦  0  ,  but  we  can  guarantee  that  e*  e  L  •  To  obtain  global  stability  in 
this  case,  requires  the  introduction  of  the  retardation  term  (3.16)  into  the 
adaptive  law  (5.1f),  see  part  { i i )  of  Theorems  1A  or  IB. 

(6)  It  is  possible  to  obtain  versions  of  Lemma  5.1  for  adaptive  systems 

of  different  forms,  e.g.,  indirect  adaptive  [5],  Also,  the  use  of 

"multipliers",  e.g.  [4],  can  be  accounted  for  as  well.  The  multiplier 

* 

effectively  makes  use  of  the  availability  of  e  as  a  signal;  and  this  allows 
rel  deg  (P*)  *  2  rather  than  1  as  required  by  condition  (1)  of  Lenina  5.1. 


6.  LIMITATIONS  IMPOSEO  BY  THE  SPR  CONDITION 


The  fact  that  the  model  error  bound  given  In  condition  (vl)  of  Lemma  5.1 
can  not  exceed  one  has  unfortunate  consequences. 

Example  1 


Consider  a  plant  with  transfer  function. 


P(s)  *  P*(s) 


ab 

( s+aH  s+b) 


(6.1) 


where  P*  Is  the  parametlrc  model,  with  two  unmodeled  stable  poles  at  -a  and 
-b.  Suppose, also,  that  b  Is  much  greater  than  a,  and  that  a  Is  much  greater 
than  the  bandwidth  of  P*(s)  .  This  situation  seems  benign  —  and  most  likely 
a  certainty.  Comparing  (6.1)  with  (5.1g)  gives. 


«.)  .,[Ai^  i1/?  >  i 

9  9  9  9  * 

/*■*..*  W.T  ik4! 


for  all  frequencies  u  >  (ab/2)1^2  ,  thus,  condition  (vl)  of  Lemma  5.1  is 
violated,  and  global  stability  cannot  be  guaranteed.  The  following  example 
illustrates  this  point. 


Example  2 


Consider  the  example  MRAC  system  (Fig.  5.1)  studied  by  Rohrs  et  al .  [12], 
where: 


P(s) 


2  229 

***"  (s+15)Z  +  4 


HR(s)  *  sTT 

A  A 

u  ■  -  Ojy  +  e2  r 
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Let  r  =  constant  and  d  a  0.  Thus,  e*  ♦  0  exponentially  when  the  tuned  gains 
are  such  that  (5.8)  Is  satisfied,  l.e. , 

T*(0)  *  1+76“  *  Hr{0)  *  1 

Even  though  (e*, ,  e*-)  exist  to  satisfy  this,  H  (s)  Is  not  SPR,  and  so 
l  i  ev 

global  stability  Is  not  guaranteed.  Simulation  runs  with  r  *  .4  and  r  =  4.0 
are  shown  In  Figures  6.1  and  6.2,  respectively.  With  the  small  input  (Fig. 
6.1)  we  see  a  stable  response  which  tracks  the  reference  very  well.  With  the 
large  input  (Fig.  6.2)  the  response  Is  still  stable,  but  large  oscillations 
are  taking  place.  Larger  Inputs  will  eventually  drive  the  system  unstable, 
e.g.  [12]. 

In  this  example.  If  the  tuned  model  Is  taken  to  be  P*(s)  *  l/(s+l)  then 
It  Is  easily  verified  that  model  error  $(u)  Is  greater  than  one  at  some 
frequency. 


7.  SPR  COMPENSATION 


In  this  section  we  heurlstlcally  develop  a  means  to  obtain  global  robust 
adaptive  control.  Since  the  SPR  condition  Is  violated  whenever  model  error 
exceeds  one,  a  natural  scheme  Is  to  construct  an  SPR  compensator  which 
allelvates  the  problems  by  ''filtering"  the  plant  output;  thus,  avoiding  the 
trouble.  However,  direct  filtering  does  not  change  the  size  of  model  error. 
For  example,  with  the  plant  p  =  (1+a)P*  ,  let  yw  denote  the  output  of  the 
filtered  plant,  where 

yw  :=  Wy  =  Wd  ♦  (1+a)WP*u  (7.1) 

Thus,  model  error  Is  uneffected.  Even  filtering  Hpv  directly  by  W  offers  no 
help,  since  the  bound  (4.4c)  is  still  less  than  one,  i.e., 

|HeJ  <  Re(W  TTev)/|W  ffgVl  <  1  (7.2) 

for  any  stable  W.  What  we  seek  is  an  SPR  compensator  which  only  effects  the 
unmodeled  dynamics,  but  leaves  the  paramtric  model  intact. 

A  compensation  scheme,  which  offers  some  promise  as  an  SPR  compensator, 

is  that  proposed  in  [22],  as  shown  in  Figure  7.1.  To  see  the  desired  result 

suppose  that  P  *  (l+A)Pm  with  a  c  B^( 6 )  .  Then,  the  compensator  is 

equivalent  to  a  plant  which  maps  (u,d)  Into  y  where 

c 

y„  *  Wd  +  P  u  (7.2a) 

c  c 

ac  ;»  -jj — —  c  BS(W«)  (7.2b) 

nt 

Thus,  whenver  a(w)  >  1  ,  select  W(s)  such  that  |w( ju») !«(<*»)  <  1  •  The  filter 
W  acts  like  a  "frequency  switch"  whose  function  Is  to  Insure  condition  (vl)  of 
Lemma  5.1. 

There  are  two  ways  to  Implement  this  compensator  In  an  adaptive  system. 

The  first  way  Is  to  use  a  fixed  model  of  the  plant  for  Pm»  1*«*»  Pm  *  P  • 

The  second  way  Is  to  replace  P  with  an  adaptive  observer,  I.e.,  p  ■  p  . 


In  either  case,  to  obtain  the  benefit  of  the  SPR  compensator,  the  signal  to  be 
controlled  Is  the  compensator  output  yc  ,  not  the  plant  output  y.  Both  of 
these  compensators  will  now  be  examined. 

Fixed  SPR  Compensator 


Let  P  =  P  ,  a  fixed  model,  and  let  the  actual  plant  be  given  by  (2.17), 
m 

P  =  (1+a)P„  with  a  e  •  Then  fixed  compensator  plant  equivalent 

model  error  (7.2b)  Is: 


A 


C 


Bs(6j) 


(7.3a) 


where 

,Tr(ju»)-P*(juJ), 

6^(u>)  |W(  ju>)  1 6 (<d )  +  |l  •  W ( jo) )  |  •  I - p — TjoD - 1  (7.3b) 

This  scheme  is  motivated  by  the  fact  that  at  low  frequencies  the  tuned 
parametric  model  P*  is  close  to  P;  thus  6  Is  small  and  W  «  1  .  At  high 
frequencies  &  is  large  but  (P-  -  P*)/P*  is  small,  W  -  0  and  so  Is 

small.  Of  course  the  compensator  Is  limited  if  there  Is  large  model  error  at 
intermediate  frequencies. 

Example  2 


Example  1  Is  modified  to  Include  a  fixed  SPR  compsnator  with  U(s)  * 
l/(s+l)  and  Fts)  *  2/(s+l)  .  Simulation  results  with  the  large  step  command 
(r«4)  are  shown  In  Figure  7.2.  Comparing  these  to  Figure  6.2,  without 
compensation.  It  Is  readily  verified  that  the  Instability  tendencies  are 
eliminated.  Also,  direct  calculations  reveal  that  Hgv  e  SPRQ  ,  thus  global 
stability  Is  Insured. 

Adaptive  SPR  Compensation 

An  adaptive  SPR  compensator,  together  with  the  adaptive  controller.  Is 
shown  In  Figure  7.3.  The  adaptive  controller  Is  described  by. 


•  ‘9czc  •  zc  8  (Fcu  •  Fc(Vr)) 


F'(s) 


rc  zc  ec  *  ec  ■  yc  -yr 
n  -1 


■  (1/Lc(s),...,s  c  /Lc(s))  ,  n£  =  deg  l£(s) 


and  the  adaptive  observer  Is  described  by, 


(7.4a) 

(7.4b) 

(7.4c) 


=  o  ' 


eo  zo 


=  /c 


IK  u.  -f:  y) 


f;(s) 


rozoeo  •  eo  -  y  -  y 
n  -1 


•  (l/LJs) . s  0  /L.lsl)  .  n  • 


deg  Lq(s) 


( 7.4d) 
(7.4d) 
(7.4f ) 


where  L  (s)  and  L  (s)  are  both  monic  and  stable.  To  generate  the  error 
o  c 

system  Interconnection  operators  associated  with  this  system,  let  e  and 

c 

6^  denote  the  tuned  parameters  with  respective  gain  errors,  0  and  0  ;  and 

*0  c  0 

let  v  :=  e'z  and  v  :=  e'z  be  the  corresponding  adaptive  control  errors 
q  c  c  0  0  0 

(3.6).  By  analogy  with  the  procedure  used  in  Section  5  we  get. 


u 


C*(r-yc)  -  1+A*  /I  vc 


(7.5) 


where 


B*i 

l —  d  *  (1  -  —  4}P*u  +  vc 
0  o 


C* 


a*2/lc 


(7.7) 


(7.6) 


b*2/l0  gN* 
p*  "  FKfc  B  "D7 


(7.8) 


and  where  (A*lt  a*2)  are  polynomials  whose  coefficients  are  the  parameters  In 
a*  ;  (B*.  ,B*0)  are  polynomials  whose  coefficients  are  the  parameters  in 

;  and  N*,  p*  and  g  are  as  defined  by  (5.6a).  The  adaptive  error  model  Is 
given  below  In  terms  of  T*,  S*,  and  as  defined  In  (5.4).  In  additon. 
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define: 


D* 

R  :*  1  ♦  (W-l)  r- 

Lo 

The  tuned  signals  are: 

e*c  *  S*(l+aRT*)"1R  d  +  (T*(l+aR)(l+aRT*)*1-Hr)r 
e*o  =  O^^d+aRT*)"^  +  D*L”^  T*a(  1+aRT* ) -1  r 

FcA*2Lclp*1,c*(  1+aRT*1“1  ( p-Rd) 

Z*c  =  _j 

FcS*(l+aRT*)  MRd-r) 

FoA*2Lclp*lK*{1+ART*rl(r-Rd) 

Z*o  *  -1 

FJJI+aRT*)  A(d  -  (l+a)r) 


(7.9) 


(7.10a) 

(7.10b) 


(7.10c) 


(7.10d) 


The  interconnections  are: 


MI+aRHI+aRT*)"1 

-(1-W)S*(1+aRT*)"1 

l+( l-WjT^D^L”1 ( 1+aRT*) 

(7.11a) 


FcP;1K*(l+aRT*)’1 

FcK*{l+AR)(l+ART*rl 


FcA*2Lilp*lK*(1-W,(1+ART*rl 

-FcS*(l-W)(l+aRT*)“1 


(7.11b) 


FoP^MI+aRT*)"1 

-F0<*n+A)(l+ART^)’1 


FoA*2Lclp*1|C*(1‘W,(1+ART*)"1 

-FqTJI-WHI+aXI+aRTJ-1 


(7.11c) 


The  factor  (1+aRT*)"*  appears  In  all  the  terms  above.  The  transfer 

function  R  (7.9)  reduces  the  effect  of  unmodeled  dynamics;  however  not  exactly 

by  the  amount  anticipated,  vis  a  vis  (7.2).  This  Is  due  to  addltonal  model 

error  Introduced  by  the  adaptive  observer.  Nonetheless,  the  model  error 

attenuation  Is  greater  than  with  the  fixed  SPR  compensator.  In  particular,  at 

low  frequencies  a  -  0  and  at  high  frequencies  R  »  0  ,  since 

W  -  0  and  -  1  .  Without  further  testing  of  H  (7.11a)  It  Is  not 

o  j  ev 

possible  to  state  that  H  e  SPR  at  intermediate  frequencies.  Note, 

ev  o 

however,  that  the  nominal  value  of  H  is: 

ev 


TT 

ev 


-(l-W)S* 

1 


(7.12) 


which  Is  SPR0  provided  that  k*  e  SPR  and 

Re  KJja.)  >  |l(l-W(ju,))S*Utt)|2  ,  a  c  R  (7.13) 

Applying  (4.11)  to  (7.11a),  a  tedious  procedure,  would  give  an  upper  bound  on 

model  error  to  Insure  H  e  SPR  . 

ev  o 
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8.  CONCLUSIONS 


This  paper  has  presented  an  Input/output  view  of  multivariable  adaptive 
control  for  uncertain  linear  time  Invariant  plants.  The  essence  of  the 
results  are  captured  In  Theorems  1A  and  2B  which  provide  conditions  that 
guarantee  global  stability.  Corollary  1  also  give  specific  L£  and  bounds 
on  significant  signals  In  the  adaptive  control  system.  These  bounds,  for 
example,  can  be  used  to  guarantee  that  the  adaptive  system  performs  as  well  as 
a  robust  ( non-adaptlve)  system  using  the  same  structure,  but  with  fixed 
gains.  By  distinguishing  between  a  tuned  system  and  a  robust  system,  we 
establish  formulae  which  can  be  used  to  restrict  the  minimum  performance 
improvement  possible  with  the  same  control  structure. 

Although  the  stability  results  (Theorem  1A,  IB)  are  not  entirely  new  (see 
e.g.,  [7], [8]),  the  input/output  setting  provides  the  means  to  directly 
determine  the  system  robustness  properties  with  respect  to  model  error.  The 
type  of  model  error  examined  can  arise  from  a  variety  of  causes,  such  as 
unmodeled  dynamics  and  reduced  order  modeling.  It  is  very  difficult  to  treat 
this  type  of  "unstructured"  dynamic  model  error  by  using  Lyapunov  theory, 
since  the  system  order  may  not  be  known  —  in  fact,  it  may  be  Infinite. 
Although  Infinite  dimensional  (distributed)  systems  were  not  considered  here. 
Theorem  1  can  be  modified  to  Include  them,  e.g.,  [26]. 

The  structure  of  Theorems  1A  and  IB  require  that  a  particular  subsystem 
operator,  denoted  Hev  ,  Is  strictly  positive  real  (SPR).  This  requirement  Is 
not  unique  to  this  presentation  -  passivity  requirements,  in  one  form  or 
another,  dominate  proofs  of  global  stability  for  practically  all  adaptive 
control  systems.  Including  recursive  Identification  algorithms. 

Unfortunately,  although  Hgv  e  SPR  is  robust  to  model  error  (Lemma  4.1),  the 
bound  on  the  model  error  Is  too  small  to  be  of  practical  use.  Even  the  most 
benign  neglected  dynamics  violate  the  bound. 

Although  this  paper  Is  concerned  with  continuous- time  systems,  the 
theorems  carry  over  virtually  Intact  to  discrete-time  systems.  This  Is  a 
direct  consequence  of  the  portable  nature  of  the  Input/output  view.  However, 
there  Is  an  Important  Issue  unique  to  discrete-time  systems:  plant 
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uncertainty  Is  critical  to  where  performance  Is  actually  measured,  which  Is  In 
contl nuous-tlme,  not  at  the  sampled-data  points.  As  a  consequence.  It  may  be 
necessary  to  map  the  discrete  portions  of  the  adaptive  system  (most  likely  the 
controller)  Into  continuous-time,  l.e.,  the  Lg-IJalns  of  the  discrete-time 
operators  In  the  Interconnection  map,  which  are  associated  with  the  adaptive 
discrete-time  controller,  would  be  needed  rather  the  discrete-time  ^-gains  • 

Another  area  worth  pursuing  Is  the  adaptive  control  of  non-linear 
plants.  The  plant  uncertainty  description  (2.11)  does  not  exclude  non-linear 
plants.  Note  that  slowly  drifting  parameters  In  an  otherwise  perfectly  known 
LTI  plant  could  yield  the  same  uncertainty  description  as  a  non-linear  plant 
approximated  by  a  parametric  LTI  model.  All  that  Is  required  is  that  there 
exists  a  (possibly)  infinite  dimensional  LTI  system  which  matches  the 
input/output  behavior  of  the  plant  for  each  possible  Input/output  pair.  Of 
course,  if  the  plant  Is  truly  non-linear,  then  the  tuned  control  is  likely  to 
be  non-linear,  which  raises  some  very  interesting  issues  for  further  research. 

One  final  remark:  the  stability  results  presented  here,  as  well  as  other 
known  results,  provide  global  stability.  This  is  achieved  by  requiring 
Hev  e  SPR  .  a  condition  which  is  difficult  to  maintain  in  normal 
circumstances.  On  the  other  hand,  this  is  a  sufficient  conditon;  violation  of 
which  does  not  necessarily  lead  to  instability.  The  simple  example  presented 
here  in  Figure  6. 1-6. 2,  Illustrates  the  point.  Other  examples  of  this 
phenomena  abound,  e.g.,  [12].  It  would  appear  then,  that  a  more  valid 
approach  to  providing  a  system- theoretic  setting  for  adaptive  control  Is  to 
develop  local  stability  conditions,  which,  hopefully,  do  not  require  that 
Hev  6  SPR  *  Pre11m1nary  results  on  local  stability  support  this  hope,  e.g., 
[33],  [34]. 
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PROOF  OF  THEOREMS  1  AND  2 

Prel iminaries 

The  main  ingredient  In  the  proof  is  to  show  stability  by  means  of 
passivity.  Although  there  are  many  variations  on  this  theme,  a  general  result 
is  given  by  the  following. 

Theorem  A.l  ([21],  [35] 

Consider  the  feedback  system  of  Figure  A.l  below  with  causal  operators 
Gj  and  G2  • 


Figure  A.l  Feedback  System 


Suppose  there  exists  real  constants  e^,  ,  1*1,2  ,  such  that 

2  2 

<x,G1x>t  >  e j I XI  ^  4  5l,Glx,t2  +  a* »  Vt  >  v  x  c  l2C0,t]  (A.l) 

for  1*1,2.  Then  the  following  holds  V  t  >  0, 

2  2 

(e2+«1)iy1lt2  4  (e1+62)iy2lt2  <  I y  1 1 12 ^ 1 U1 1 12  *  >u2*t2^ 

2  2 

♦  «y2lt2(,Vt2  *  2|c J !•■  Uj1 12)  4  Uji-iUji ^2  ♦  I e 2 1  • 


4  la,  |  ♦  la, | 


(A. 2) 


Proofs  of  both  theorems  also  rely  on  well  known  results  for  systems 
H  c  SgXm  •  The  results  required  here  are  summarized  In  the  following. 

Theorem  A- 2  [see  [19],  Thm.  9,  pg.  59] 

Let  H  e  Snxm  ;  then: 
o 

(1)  If  u  e  L™  ,  then  y  =  Hu  e  ^  e  L^  .  y  Is  continuous,  and 

y(t)  ♦  0  as  t  ♦  •  . 

(ii)  If  u  e  Lm  .  then  y  =  Hu  e  Ln  .  $  e  Ln  .  and  y  Is  uniformly 
continuous. 

(Ill)  If  u  t  L™  and  u(t)  +  constant  c  e  Rm  as  t  ♦  •  ,  then 

y ( t )  ♦  H(0)c  exponentially  as  t  ♦  •  . 

In  order  to  simplify  notation  we  drop  the  superstrlct  on  Ln  which 

P 

indicates  vector  size. 

We  will  establish  Theorem  1A  first.  Some  of  the  steps  will  be  repeated 

for  IB.  Also,  without  loss  of  generality,  the  matrix  r  In  the  adaptation  law 

(3.15) ,(3.16)  Is  set  to  Identity.  Corollary  1  Is  established  as  a  by-product. 

Proof  of  Theorem  1A 


Part  (1) 


Identify  GJt  G2  In  Figure  A.l  with  e  v  and  HftV  respectively.  Also, 
let  Uj  «  e^,  u^  «  0,  Cj  ■  e,  •  v,  and  y?  «  H^v.  .  Using  adaptive 

law  (3.15)  we  obtain. 


<e,v>T  *  <e,Z'e>T  «  <Ze,e>T  ■  <e,  0>T 
■  y  ie(T)i2  -  ^  ie(0)i2 


(A. 4) 
(A. 5) 


1  ...  Z 


I 


Thus,  using  (A.l)  gives. 


£1  ‘  S1  *  °1  ‘  ‘  7  « 0 ( o ) I ^  (A. 7) 

Since  G_  *  H  e  SPR  by  assumption,  3u,  y  >  0  such  that  V  x  e  L,  , 
<x,Hevx>T  >  u 1 XB T2 *  ,HevX,T2  <  YlX,T2  *  Hence‘  from 


e2  *  M •  ~  ®2  *  ^  1A.8) 

Using  Lemma  A.l,  together  with  (A.4MA.8)  gives, 

,V,T2  <  I^,e*,T2  +  (,e*'T2  +  2wl9(0)|2)1/2]  (A. 9) 

ie-e^iT2  <  y«v«T2  (A. 10) 

|e <T) 1 2  <  1 0 (0 ) 1 2  +  2ieij2  *vly2  (A. 11) 


The  bounds  shown  In  (4.8)  follow  using  the  assumption  e*  e  L2  •  Hence, 
e,v  t  L,  and  e  e  L  . 

C  00 

• 

Having  established  that  v  e  L2  ,  Theorem  A-2  ^  z:  =  z-z*  e  L2nL_,  z  e 

z  ♦  0,  and  z  Is  continuous.  Since  z*.  z*  e  l  by  assumption.  It  follows 

that  z  e  L  and  i  e  L  (*♦  z  is  uniformly  continuous).  Using  v  *  Z'0  with 
00  00 

z,  0'  e  L  V  e  L  .  Using  e  *  e*-H  v  with  e*  e  L  and  H  e  S  (by 
assumption),  and  v  c  L  =>  e  e  L  .  Hence,  5  =  Ze  e  l  8  is  uniformly 

mm  00 

continuous  -0  v  *  Z'0  Is  uniformly  continuous  (since  z  1s)*4  v  +  0  since 
v  e  l2  Is  established.  Using  v  ♦  0^  e  -  e*  0  ,  and  since  e*  ♦  0  by 
assumption,  e  ♦  0  .  Furthermore,  v  ♦  z  *  0  exp.  and 
0  *  Ze  *  Ze  ♦  Z*e  ♦  0  ,  because  z  and  e  ♦  0  .  Using  v  «  Z'0  +  Z’§  with 
z,  I.  i  t  1.4  i  *  .  Hence,  e*  =  6*  -  Hey  h  ,  because  e*  e  by 

assumption.  Thus,  e  •  2e  ♦  Ze  e  L  •  This  establishes  properties  (1-a)- 
( 1-d) . 


To  show  (1-e)  consider  (3.15)  written  as: 


(A. 12) 


J 


o  *  -Z#  h  z*  e  ♦  w 
*  ev  * 

-(Z*HevZ‘  ♦  Z  Hey  Zi  ♦  Z  Hev  Z’)e 

Since  we  have  already  established  that  2  ♦  0  exp.  and  eel  .  It  follows 
that  w  ♦  0  exp.  Since  z#  e  PE  by  assumption  (provided  e*  *  0)  ,  w  (-►o  is 
stable  by  (2.15).  Hence,  e  -►  0  exp.  0,  v  ♦  0  exp.-^  e-e*  ♦  0  exp.  This 
completes  the  proof  of  part  (1)  with  adaptive  law  (3.15). 


To  show  that  (l-a)-(i-a)  hold  with  adaptive  law  (3.16)  requires  showing 
that  Gj:e  ->  v  Is  passive.  Consider  the  typical  time  Interval, 


I  = 


j  Ij  *  {t  e  CtQ.tj)!  ie(t)i  <  c} 

I  I2  *  {t  e  C 1 1 , t2 ) I  *9(t)i  >  c  >  maxie^i} 


(A. 13) 


Hence, 


Thus, 


<e,v>  »  <e,v>  +  <e,v> 

1  *1  l2 


(A. 14) 


<e,v>lj  *  <®*  0>i1  “  7  "  \  ■  9 ( tQ )■ 2 

•  *  2  A 

<e,v>.  «  <0  +  (1  -  lei/c)  0,  e>T 

t2  i2 

*  y  ie(t2)i2  -  j  ie( tj )■ 2  +  (1  -  ie»/c)2<0,  e>j 

! 

1  ..  .  2  1  . .  .  2 
>  f  leCtj,)!  -  j  i9(tj)i 


(A. 15) 
(A. 16) 
(A. 17) 
(A. 18) 


because  <e,  e>.  >  0  from, 

*2 
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•  u  * 


**  .■*  ^  •  •  *  . 


>  _  k  A.  k  _  *  _  ^.1 


exp 


•Y’-w-t  v— v;- 


■»  ■  -s  rr* 


m  » 's  •  v  -  v  —  *  -  v  *rw 


tw  r?  i  tj*  "  y 1  •  -jr  ytyj  v — wr 


I 

► 


■ 


■> 


«r 


I 

a 


i 

c’ 

V 

i 


e(t)'  e( t)  «  0( t> '[e(t)-e*] 

*  2  * 

■  ie(t)i  -  e(t)'e* 

•*  2  “ 

>  I0(t)l  -  l0(t)|C 

=  ie(t)i (ie(t)i-  c)  >0,  Vt  e  I2  .  (A. 19) 


Thus, 


<e,v>j  >  ^  ie(t2)i2  -  ^  ■ 6 C tQ ) i 2  (A. 20) 

Repeating  the  above  procedure  recurslvely.we  eventually  conclude  that 
1  2 

<e,v>T  >  -  y  i0(O)i  as  before  (A. 6),  and  hence,  G^e  l-^-v  Is  passive.  The 
results  In  (1)  now  repeat  for  adaptive  law  (3.16).  This  completes  the  proof 
of  part(i). 


Proof  of  Theorem  lA.Part  (11) 


Theorem  1A,  Part  (11)  Is  essentially  an  L  -stability  result.  The  method 
of  proof  requires  the  notion  of  "exponential  weighting"  which  Is  a  means  to 
obtain  L^-stabllity  of  a  system  from  the  Instability  of  an  exponentially 
weighted  version  of  the  system  (see  e.g.,  [19],  Chapter  9).  We  require  the 
following: 

Definition:  Given  a  real  number  a  define  the  exponential  weighting  operator 
by 


xa(t)  :•  eatx(t)  (A. 21) 

Consider  the  system  y  «  Gu.  An  exp.  weighted  version  of  this  system  Is 
denoted  by  y°  :«  Ga  ua  .  Note  that  If  G  Is  a  convolution  operator  with 
transfer  function  G(s)  then  G°  Is  also  a  convolution  operator  with  transfer 
function  G(s-o)  .  Thus,  the  corresponding  exponential  weighted  error  system 
corresponding  Is  described  by 


4 


4 


s.' 
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where  a  >  0  such  that 


0a  =  ae°  +  Zea  -  p(0)ea 


Ha  £  SPR™  and  Ha  c  Skxm  (A. 23) 

ev  +  zv  o 

Using  Theorem  A-l,  Identify  G^  with  e“  v°  and  G^  with  H^y.  Note  that  It  Is 
always  possible  to  find  some  a  >  0  such  that  (A. 23)  holds.  We  now  examine 
the  passivity  of  G^:  e°  ♦  v“  .  Thus, 

<e°,  va>T  *  <ea,  Z’oa>T  =  <  Zea,  ea>T 

=  <ea,  0a  -  a0a  +  p(0)0a>T 

*  j  e2aTi0(T)i2-  j  » e (0)i 2  +  <p{e)e°,  ea>-j-oieai22 
>  €^*^10 (T)i ^  -  2"  I0(O)|2  -  al8°|22  (A. 24) 


The  last  line  follows  from  (A. 19),  hence,  (A. 24)  holds  with  or  without  the 
retardation  term  In  the  adaptive  law.  At  this  point  there  are  two 
possibilities:  either  e  e  or  |e(t) |  ♦  •  as  t  ♦  •  .  If  8  c  La  then  3 
constant  c  <  •  such  that  iei  <  c  .  Then, 

O  m  O 

<e®,  v“>T  >  j  c2aT(ie(T)i2  c2)  -  ie(o)i2 

1  2aT  2  1  . o|n.  2  (A. 25) 

>  -  y  E  co  "  7 

If  le(t)|  ♦  •  as  t  ♦  -  then  It  Is  always  possible  to  select  an  arbitrarily 
large  T  such  that  ie(T)i«  I8IT-  •  Hence,  for  this  T,  (A. 24)  becomes, 

<e°,  v“>T  >  \  c2aT(i e(T)i 2  -  iei2J  -  £  ie(0)i2 


Thus, for  some  arbitrarily  large  T,  (A. 25)  and  (A. 26)  have  the  general  form. 


1  .e. , 


<ea,  va>T  >  -Cj  e2oT  -  C2  :*  -c(aT) 


where  c^,  are  non-negative  constants.  Hence, 


Cj  a  5j,  *  0,  aj  *  -c(aT) 


Since  =  H“  e  SPR+,  3  constants  v  ,  y  >  0  such  that 


<x.  «Sv  X>T  >  wIX,T2 
,Hev  X,T2  *  Y  |X,T2 


Then, 


e2  *  w  ’  62  =  a2  *  ^ 


Using  (A. 2),  we  get 


ivaiT2  < — ^ — [ I T2  +  (* e*» T2  +  Zv  ct«TM1/2] 
2u 


Since  e*  e  L  by  assumption, 

w  m 


aT/,  >-1/2. 


.e“«T2  <  ea'(2a)  le*^ 


Thus, 


a  el  £a 

1 T2  «  - 5T 


«T(2a)*l/2  .2  .  ,  -2.T  .,,1/2, 


[*«★*.  +  +  c  pC(oT))  ] 


Since  e  S«xm  ,  we  obtain 


zv  * 


|z(T)|  *  |  J  H2V(T-t)v(t )dr | 


|e‘oT  ^H5va-T)v0(T)dT| 


(A. 27) 


(A. 28) 


(A. 29) 


(A. 30) 


(A. 31) 


(A. 32) 


(A. 33) 


(A. 34) 


(A. 35) 


'  »  *  •  •  k  *  .  k  ‘k*-  , **  ,  -  i  * 


(A. 36) 


<  *  ,V<1,T2 


where  H“y(t)  Is  the  Impulse  response  matrix  associated  with  Ha 

Substituting  (A. 33)  and  (A. 27)  Into  (A. 36)  and  noting  that 
-2  T 

e  “  c(«T)  <  Cj  +  Cg  ,  we  obtain. 


I  z  ( T ) 


(2a) 


-1/2 


•Hzv<' 


)l 


[ 1  e* 


(,e*l_  +  4. 


ay(c J+C2)) 


1/2 1 


(A. 37) 


Since  the  right  hand  side  Is  Independent  of  T,  and  since  T  can  be  selected  to 
be  arbitrarily  large.  It  follows  that  z  e  L  .  Assuming  there  is  no 

oo 

retardation  or  persistent  excitation,  this  completes  the  proof  of  (11 -a )  to 
(ii-d). 


Assume  now  that  z  e  PE  ,  which  is  a  noncontradictory  assumption  since  we 
have  already  shown  that  z  e  L  .  Hence, 

9  =  -2  Hey  2'  e  +  Ze*  (A. 38) 

Since  z  e  PE,  H  e  SPR  and  z,  e*  e  t  ,  It  follows  from  (2.15)  that 
ev  +  «> 

(Ze*.  e(0))  Is  exp.  stable,  thus,  e,  0  £  L  .  The  remaining  results  In 

00 

(11 — e)  follow  Immediately. 


Suppose  now  that  the  adaptive  law  is  given  by  (3.16).  Then,  we  can 
write. 


0  *  2  e  -  p (9 )a  =  Z[e*-HevZ'{e-e*)]  -  p(e)e 
*  w  -  z  Heyz'  0  -  P(e)0 


(A. 39) 


where  w  :=  Z  e.  ♦  Z  H  Z'e*  £  l  ,  because  z,  e*  £  L  •  Consider  the 

ev  *  •  „  «  *  » 

candidate  Lyapunov  function  V:t  h-*ie(  t)i*  •  Hence, 

V  -  2  w'e  -  0*  Z  HevZ’0  -  p(e)V  (A. 40) 

Suppose  ie(t)i  ♦  •  as  t  ♦  •  .  Then  there  exists  a  time  T  >  0  such  that 
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(A. 41) 


VT  <  2IWI.  v‘/2.  •  I.f  Y.(He<)VT  -  (1  -  v{/2/c)\ 

Clearly,  there  exists  a  finite  constant  Cj  such  that  when  VT  >  Cj,  V  <  0  . 
Therefore,  9  can  not  grow  beyond  all  hounds,  and  hence,  0  E  L  •  So  then  Is 
e  and  §  ,  and  again  the  result  of  (1i-e)  follow.  This  completes  the  proof  of 
Theorem  1A.  Note  that  in  this  case  we  do  not  obtain  specific  bounds  on  e, 
because  the  proof  proceeds  by  contradltion. 

Proof  of  Theorem  IB 


Part  (1) 


Since  H  e  SPR  ,  there  exists  q  >  0  such  that  G  :=  (1  +  qs)H  c  SP 
ev  .  o  ev  ev 

and  furthermore,  G~*  e  S  .  As  a  result  we  can  write  (3.14a)  as, 

e  *  -Hev  y,  y  s  V  -  G‘J(e*  +  q  ej  (A. 42) 

Referring  to  Lemma  A-l,  let  Gj  :  vh^e,  G^  *  Hgv  ,  Uj  *  0  ,  and 

u-  *  -G’^e*  +  qe*)  .  Using  (A. 2)  together  with  (A. 42)  and  the  passivity 

c  ev  *  * 

properties  of  H  qives, 
ev  * 

,eiT2  *  Tyt ,U2IT2  +  (,U2*T2  *  2mI®(0)|2)1/2]  (A. 43) 

le (T) |  <  |e (0 ) |  +  2ieiT2  .  IU2IT2  (A.44) 

where  u  Is  defined  In  (4.9a).  Using  (4.9b)  gives, 

1 u2i t2  <  (l/k)ie*  +  qe*ij2  .  This  together  with  (A. 43),  (A.44)  and  the 

assumption  e*,  e*  e  l2  gives  the  bounds  shown  In  (4.9).  Hence, 

e  c  L„,  e  c  L  .  However,  we  can  not  conclude  that  v  e  L«  as  In  Theorem  1A, 
cm  <  Z 

part  (1).  From  (A. 42),  we  can  conclude  that  (1  +  qs)  v  e  l2  •  Since 

G  :■  (1  ♦  qs)H  e  S  ,  it  follows  from  lemma  A-2  that 

,zv  zv  ot 

z  :*  z-z*  c  l2  H  lB  ,  z  e  L2  and  z  ♦  0  .  Repeated  use  of  Lemma  A-2  and  the 
error  equations  (3.14)  gives  the  results  (1-a)  -  (1-d).  (1-e)  follows  from 

the  arguments  in  the  proof  of  Theorem  1A,  part  (1). 


Part  (I i) 


The  proof  Is  entirely  analgous  to  that  of  Theorem  1A,  part  (11),  where 
again  we  use  exponential  weighting. 
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APPENDIX  B 


PROOF  OF  LEMMA  5.1 

The  proof  utilizes  the  following  known  results: 

Deflnl tlon :  Let  J  denote  a  subset  of  S,  consisting  of  functions  In  S  whose 
Inverse  Is  also  In  S. 

Fact  [29]:  If  G  is  any  scalar  transfer  function  in  R(s),  then  G  has  a  cop rime 

factorization  In  S,  l.e.,  there  exists  N,  D,  A,  and  B  In  S  such  that 

G  =  N/D  and  AN  +  BD  a  1. 

Lemma  B-l:  Consider  the  tuned  adaptive  system  of  Figure  5.2.  Let 
P*  e  RQ(s)  and  c*  e  R0( s)  have  coprime  factorizations  In  S  given  by 

p*  *  Nn /D  and  C*  s  Nr/Dr  .  respectively.  Then,  the  elements  of  the 

transfer  matrix  from  (r.d)  Into  (e*.  z^.y,  u)  all  belong  to  S,  if: 


(1)  Q  :=  DO  +  NN  e  J  , 
pc  pc 


(from  [29]) 


(B.l ) 


(11)  6 (b>)  |T*(  ju) )  |  <  1,  V  «  e  R  , 


(from  [16]) 


where 


T*  :*  NpNc/Q  P*C*(1+P*C*V 


(B.2) 


Using  the  definition  of  Q  we  can  write  Hpv  and  Hjv  from  (5.5)  as. 


Hev  “  NpQ”1( 1+a)( 1+aT*)"1 


(B.3) 


FDpQ-^l+AT*)-1 


2V  FNpQ^d+aJd+AT*)'1 


(B.4) 


From  the  definition  of  K*  (5.4b),  we  also  obtain 


0 


(B.3) 


Proof  of  Lemma  5.1 

We  first  show  that  (1),  (11),  and  ( i v )  *>  Q  e  J  .  Let  P .  *  N  /D  be  a 

P  P 

coprime  factorization  of  P*  such  that  rel  deg  0p(s)  a  0.  Since  (1)  a>  rel  deg 
P*(s)  *  1,  It  follows  that  rel  deg  Np(s)  *  1.  Moreover,  (1v)  => 
rel  deg  K*(s)  *  1,  and  that  iq(s)  and  Kgts)  are  stable.  This,  together  with 
(11)  and  (B.3)  establishes  that  Q  e  J  . 

c  Sq  follows  Immediately  by  Inspection  of  (B.2),  since:  F  e  SQ  by 

assumption;  0  ,  N  e  S  ;  Q  e  J;  a  e  S  by  assumption  (vl);  and  finally  (vl) 

P  P  i 

*>  (11)  of  Lemma  B-l  *>  (1+aT*)"1  c  S  . 

Conditions  (iv)  and  (vl)  *>  H  e  SPR  .  This  follows  from  Lemma  4.1 

ev  o  , 

by  letting  TT  *  K .  and  letting  1  +  fl  «  ( l+a)( 1+aT*)~1  .  Thus,  (4.4a)  Is 
ev  *  ev  * 

satisfied  since  K.  e  SPR  from  (Iv).  Also,  from  (4.4b), 

*  o 

Mu)  *  |flev(ju)|  *  U(ju)S*(ju)[l-4(ja>)T*(ju)]‘1|  (B.4) 


6(u )  |S*(ja>)  I 


ld»M  l*Ud)l 


<  Mu»)  *  n(a>) 


(B.5) 


The  last  Inequality  comes  from  conditions  (vl)  and  the  definition  of 
F(u)  from  (4.4b). 


The  final  step  In  the  proof  of  Lemma  5.1  is  to  show  that  there  are  a 
sufficient  number  of  parameters  In  e*  to  Insure  a  solution  exists.  This  Is 
guaranteed  by  satisfaction  of  condition  (v).  To  see  this  combine  (B.3)  with 
the  definition  of  0  from  (B.l)  to  get 

Q  N  N  +0  D  -  NIC1  (B.6) 

c  p  pc  p  * 
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From  (5.2),  let  N  «  A../L  and  D  «1  4  A*.  be  a  coprime  factorization  of 

C  c  C  1/L 

C*,  and  let  Np  ■  g  N*/l  and  0p  *  1  +  D*/L  be  a  coprime  factorization  of  P*» 
where  P*  is  as  defined  In  (1).  With  K*  given  by  (1v),  (B.6)  becomes  the 
polynomial  equation, 

am  Ki0*  +  A*2KlN*  *  L{K2N*"Ki°*)  (B.7) 

Since  deg(K2Nw)  «  degUjD*)  and  K^,  N* ,  and  D*  are  all  monlc.  It  follows 

that  deg[L(K _N  -K  D  ]  *  deg(L)  +  deg(K  )  ♦  deg(0.)  -  1  .  Then,  using  known 
Z  *  1  1 

results  on  polynomial  equations,  e.g.  [30],  It  can  be  shown  that  (v)  Implies 
that  (B.7)  has  a  solution  (A*j,  A*2)  . 
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